Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,436,956 2 bp→TT coding (202‑203/741 nt) cysQ → 3'(2'),5'‑bisphosphate nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,436,9560CT97.2% 235.1 / NA 71P68F (CCC→TTC) cysQ3'(2'),5'‑bisphosphate nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (42/27);  minor base G (0/2);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,436,9570CT98.6% 225.5 / NA 71P68F (CCC→TTC) cysQ3'(2'),5'‑bisphosphate nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (41/28);  minor base G (0/1);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATGCGCGCC  >  NC_000913/4436713‑4437199
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
aaaaCAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTg                                                                                                                                                                                                                                              <  2:81209/251‑1 (MQ=255)
 aaaCAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTgg                                                                                                                                                                                                                                             >  1:140721/1‑251 (MQ=255)
  aaCAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTggg                                                                                                                                                                                                                                            <  1:348981/251‑1 (MQ=255)
      aTAAGTTAACACCGCTAACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGt                                                                                                                                                                                                                                        >  1:114595/1‑251 (MQ=255)
       tAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGCTGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCCGCGGATATTGCCGCTCACCCCGTTATCAGGGACGGTTTACGTACGCTGACACCGGATGTTACGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAGGTc                                                                                                                                                                                                                                       >  2:5805/1‑251 (MQ=255)
             aaCACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCa                                                                                                                                                                                                                                  <  1:52983/250‑1 (MQ=255)
                  cGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGGCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTg                                                                                                                                                                                                                            >  1:1385/1‑251 (MQ=255)
                   gCTCAAAGAGACGAGGAGGAGCAATGTATGATCAAGTATAACAGCTTGCACCGAATGCAGGCGATGCCATTCAGCAGGTCAACGACGGGACGAAACCGATGGACGTCGTCAGCAAAACGGACAATTCTCCGGCAACGGCAGCGGATATTGCAGATCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGCCCTTTCTCAAGCAGCTCCTTTCGGTTGGGCAGTCCGTCAGCACTgg                                                                                                                                                                                                                           >  2:465391/1‑251 (MQ=255)
                                   tGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTActggctgg                                                                                                                                                                                                           >  2:461050/1‑251 (MQ=255)
                                               aGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATTCTTTCGGTTGGGAAGTCCGTCAGCACTGGCCGCGTTACTGGCTGGTAGACCCGCTgg                                                                                                                                                                                               >  2:264912/1‑251 (MQ=255)
                                                gATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGCTGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTgg                                                                                                                                                                                               >  1:263430/1‑250 (MQ=255)
                                                 aTCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGAAACGGCAGCGGATATTGCCGCTCAAACCGTTATCATGGACGGTTTACGTACGCTGACACCGAATGATCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGGCCGTCAGCACTGGCAGCGTGACTGGCTGGTAGAACCGCTGGAt                                                                                                                                                                                             >  2:251990/1‑251 (MQ=255)
                                                  tCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGCAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATg                                                                                                                                                                                            >  1:64748/1‑251 (MQ=255)
                                                    aaGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATggg                                                                                                                                                                                          >  2:37937/1‑250 (MQ=255)
                                                     aGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGGTGACCCGATGGATGGt                                                                                                                                                                                          >  2:127167/1‑250 (MQ=255)
                                                        aTGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACTCCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACt                                                                                                                                                                                       <  2:167484/250‑1 (MQ=255)
                                                          gtcAGCTTGCACCGATGGCAGTCTATTCCATTATGGAAGGCTACGACGGTGCGAAACCGATAGTCGACGTTCGCAAAGCGGATAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTAATCATGGACGGTTTACTTTCGCTGGCACCGGATGTTCTGGTCTTTTCTCAAGAATATCCTTTCGTTTGGGAATTGCGTCAGCACTGGCGGCGTTACGTTCTCGTAGACCCGCTGGATTGTACTaa                                                                                                                                                                                     <  2:54138/248‑1 (MQ=255)
                                                            cAGCTTGCACGGAATGCAGGCGATGCCATGATGCATGTCTACGACGGGACGAAACCGATGGACGTCGGCAGCAAAGCGGACAATTCTCCGGTAACTGCAGCGGATATTGCCGCTCACACCGTTATCAGGGACGGTTTACGTACGCTGACACCTGATGTTCCGGTCGTTGCTGAAGAATATCGTTTCGGTTGGGAAGACCGTCAGCACTGGCAGCGTTACTGGCTGTTAGACGCGCTGGATGGTGCTAAaga                                                                                                                                                                                  <  2:432189/251‑1 (MQ=255)
                                                             aGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGTTGGTAGACCCGCTGGATGGTACTAAagag                                                                                                                                                                                 <  2:449224/251‑1 (MQ=255)
                                                               cTTGCACGGAATGCAGGCGAAGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGtt                                                                                                                                                                               >  1:148501/1‑251 (MQ=255)
                                                                 tGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGtt                                                                                                                                                                               <  1:111254/249‑1 (MQ=255)
                                                                  gCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTttat                                                                                                                                                                            <  2:404652/251‑1 (MQ=255)
                                                                  gCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTtta                                                                                                                                                                             <  1:37937/250‑1 (MQ=255)
                                                                  gCACGGAATGCAGGCGATGCCATTATGAAGGTCTACGACGGGACGAAACCGATGGACGTCGTCATCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCAATGGAAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTttat                                                                                                                                                                            <  2:134229/251‑1 (MQ=255)
                                                                   cACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTttatt                                                                                                                                                                           <  1:457790/251‑1 (MQ=255)
                                                                      ggAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACACTTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGCCGGTTTACGTACGCTGACACCGGCTGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTACAGAGTTTATTaa                                                                                                                                                                         >  2:388542/1‑250 (MQ=255)
                                                                       gAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTaa                                                                                                                                                                         >  1:31082/1‑249 (MQ=255)
                                                                       gAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTaa                                                                                                                                                                         <  2:31082/249‑1 (MQ=255)
                                                                           gCAGGCGATGCCATTATGCAGGTCTACGACGGGCCGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTa                                                                                                                                                                    >  1:127219/1‑250 (MQ=255)
                                                                                 gATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCg                                                                                                                                                              >  1:246846/1‑250 (MQ=255)
                                                                                    gCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAAtt                                                                                                                                                          <  1:321947/251‑1 (MQ=255)
                                                                                     ccATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCCGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTc                                                                                                                                                         >  1:273453/1‑251 (MQ=255)
                                                                                                       acgGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACTGCAGCGGATATTGCCGCTCACTCCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGAGCCTGTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCtt                                                                                                                                       <  2:399028/251‑2 (MQ=255)
                                                                                                         gggACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCAt                                                                                                                                     >  1:361696/1‑251 (MQ=255)
                                                                                                            aCGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGCTATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTCGACCCGCTGGATGGGACTAAAGCGTTTATTAAACGTAATGGCGAATTTACCGTTACCATTGAGCACATTg                                                                                                                                   >  2:39933/1‑250 (MQ=255)
                                                                                                              gAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCCCATTGAcc                                                                                                                                >  2:374938/1‑251 (MQ=255)
                                                                                                                 aCCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATg                                                                                                                             <  1:325817/251‑1 (MQ=255)
                                                                                                                 aCCGATGGACGTCGTAAGCAAAGCGAACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTATGCCCGCTGGATGGTACTAAAAGGGTTATTTAACGTAATGGCCAATTCCCCGTTAACATTGCGCTCATTGCCCa                                                                                                                               >  2:384547/1‑249 (MQ=255)
                                                                                                                   cGATGGACGTCATCAGCAAAGCTGACAGTTCGCCGGTAACGGCAGCGGAAATTGCCGCTCAATCCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGACAGCATTGGCATCGTTAATGGCTGGTAGTGCCGCTGGGTGGTAATAAAGAGTTTATTTAACGTGAGGGCGGATTCACCGTTGACATTGCGCACATTGTCCATggt                                                                                                                           <  2:271690/251‑2 (MQ=255)
                                                                                                                         aCGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCAt                                                                                                                                     <  2:121562/235‑1 (MQ=255)
                                                                                                                         aCGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCAt                                                                                                                                     >  1:121562/1‑235 (MQ=255)
                                                                                                                           gtcgtcCGCAACGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCAAACCGTTATCATGGACGGTTTACGTACGCTGACCCCGGATGTTCCGGTACTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGAt                                                                                                                   >  2:191940/1‑251 (MQ=255)
                                                                                                                                cAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGAtttt                                                                                                                >  1:152777/1‑249 (MQ=255)
                                                                                                                                    aaaGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACCCCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCgt                                                                                                          >  1:416633/1‑251 (MQ=255)
                                                                                                                                      aGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCgtgg                                                                                                        <  1:284413/251‑1 (MQ=255)
                                                                                                                                                                 ggATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGt                                                                             >  1:430918/1‑251 (MQ=255)
                                                                                                                                                                               acacCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGa                                                                >  1:252782/1‑250 (MQ=255)
                                                                                                                                                                                 acCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAgg                                                             <  1:293398/251‑1 (MQ=255)
                                                                                                                                                                                     ttATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCaa                                                          >  1:158671/1‑250 (MQ=255)
                                                                                                                                                                                       aTCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTGCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAAGGAACGTAATGTACAGCGCGGCAGAAGGCAAAGc                                                       <  2:361696/251‑1 (MQ=255)
                                                                                                                                                                                        tCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAgg                                                             <  1:281641/244‑1 (MQ=255)
                                                                                                                                                                                        tCATGGACGGGTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGACCGTAATGTACAGCGCGGCCGAAgg                                                             >  2:281641/1‑244 (MQ=255)
                                                                                                                                                                                           tGGACGGTTTACGTACGCTGACACCGGATCTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGACGTCCGTCAGCACTGGCCGCGTTACTGGCTGGCAGACCCGCTGGATGGTACTAAAGCGTTTATTAAACGTAATGGCGAATTCACCGTTAAAATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTATTGAACGTAATGTACAGCGCGGCAGAAggaacagtgtgg                                                   >  2:39235/1‑241 (MQ=255)
                                                                                                                                                                                             gatggtatacgaacgCTGACTCCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTCAGGATGCGAAGACCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGCTGGTCCTAAAGAGTTTATTAAACGTACTGGCGAATTCACCGTTAACATTACGCTCATTGACCATGGCAAACCGAATTTAGGCGTGGTGTAATCGACGGAAACGAAAGAAATGTACAGAGGGGGAGAAGGccaacagtgga                                                  >  2:49268/8‑240 (MQ=255)
                                                                                                                                                                                              aCGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGaa                                                 >  1:204431/1‑250 (MQ=255)
                                                                                                                                                                                                ggTTTACGTACGCTGACACCGGATGTTCAGGTCCTTTCTGAAGCAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGGACCGCGCGGCAGAAGGCAAAGCGTGGaa                                                 >  2:166827/1‑248 (MQ=255)
                                                                                                                                                                                                 gTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTGCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAaagaa                                             <  1:264912/251‑1 (MQ=255)
                                                                                                                                                                                                      cgtacgCTGACACCGGATGTTCCGGTACTTTCTGAAGAAGATCCTTTCGGTTGGGCAGGCCGTCAGCACTGGCAGCGTTACTGGCTGGGAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGACGAATTCACCGATAACATTTGGCTCATTGAACATGGCAAACCGATTTGAGGCGTGGTGTATGAGCCGGGAATGAACGTAATGGAAAGCGCGGCAGAAGGGAAAGCGTGGAAAGAagaggg                                        >  2:267310/1‑249 (MQ=255)
                                                                                                                                                                                                       gtacgCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGGACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGCCCATGGCAAACCGATTTTAGGCGTGGTGTCTGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGAGGAAAGAAGAGTgg                                       >  2:247/1‑250 (MQ=255)
                                                                                                                                                                                                          cgCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCggg                                    >  1:217643/1‑250 (MQ=255)
                                                                                                                                                                                                              gACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTGTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTgc                                <  1:196841/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccGGATGTTCCGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCaa                            <  2:343240/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGGTCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGCTTc                    >  2:355395/1‑249 (MQ=255)
                                                                                                                                                                                                                               tCCTTTCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTc               >  1:94138/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCTGAAGAAGATCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCcgcg           <  1:374938/250‑1 (MQ=255)
                                                                                                                                                                                                                                         agaagaTCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATgcgc     >  1:229732/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gaagaTCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATgcgcg    <  2:152777/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gaagaTCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATgcgcg    <  2:217643/251‑1 (MQ=255)
                                                                                                                                                                                                                                            agaTCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATGCGcgcc  >  1:87048/1‑251 (MQ=255)
                                                                                                                                                                                                                                            agaTCCTTTCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGGAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGAGTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATGCGcgcc  <  1:230453/251‑1 (MQ=255)
                                                                                                                                                                                                                                            agaTCCTTTCGGATGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACACGCTGGATGGTACTAAAGAGTTTATTAAACGCCATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGCAGGCAAAGCGGGGAAAGAAGAGTGCGGGGTGCGCAAGCAGATTAAAGTCCGCGATGCGcgcc  >  2:291261/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
AAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTGTGCGCAAGCAGATTCAGGTCCGCGATGCGCGCC  >  NC_000913/4436713‑4437199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑