Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,045,994 G→T 7.5% intergenic (‑93/+174) ygfF ← / ← gcvP putative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,045,9940GT 7.5% 319.1 / 12.8 118intergenic (‑93/+174)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base G (42/66);  new base T (4/5);  total (46/72)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 119 total aligned reads displayed.

TAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCG  >  NC_000913/3045845‑3046140
                                                                                                                                                     |                                                                                                                                                  
tAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCg                                                                                                                                                    <  1:797740/150‑1 (MQ=255)
     cAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCtt                                                                                                                                               <  1:683832/150‑1 (MQ=255)
      aGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCttt                                                                                                                                              >  1:1055312/1‑150 (MQ=255)
      aGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCttt                                                                                                                                              >  1:44417/1‑150 (MQ=255)
        ttGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACACTTTTCCCCCTCGCCCTTTCa                                                                                                                                            <  2:1432976/150‑1 (MQ=255)
         tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTAAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAg                                                                                                                                           >  1:127871/1‑150 (MQ=255)
          gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAgg                                                                                                                                          <  1:966580/150‑1 (MQ=255)
          gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAgg                                                                                                                                          >  1:704093/1‑150 (MQ=255)
           cccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACATATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAggg                                                                                                                                         >  1:372817/1‑150 (MQ=255)
            ccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTAAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGga                                                                                                                                        <  2:1079732/150‑1 (MQ=255)
              gcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGgaga                                                                                                                                      <  2:586327/150‑1 (MQ=255)
              gcccaATGCCTCGACTGCCACCAGGCACTAGTGCTATAGCCATTTGTTTTTACGAGAAAAATTAATTCGTTGAGGTATCTTTGAGATTGCTTACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGgaga                                                                                                                                      <  1:1036818/150‑1 (MQ=255)
               cccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGagag                                                                                                                                     >  1:103711/1‑150 (MQ=255)
               cccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGagag                                                                                                                                     <  1:1422297/150‑1 (MQ=255)
                 cgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCTCCCTGAAAGGGCGAggg                                                                                                                                   <  1:262217/150‑1 (MQ=255)
                 cgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAggg                                                                                                                                   >  2:679534/1‑150 (MQ=255)
                  gATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCTCCCTGAAAGGGCGAgggg                                                                                                                                  <  2:592575/150‑2 (MQ=255)
                  gATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGc                                                                                                                                  <  1:207386/150‑1 (MQ=255)
                   aTGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCTCCCTGAAAGGGCGAggggg                                                                                                                                 <  2:1571910/150‑3 (MQ=255)
                    tGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCg                                                                                                                                <  2:1566662/150‑1 (MQ=255)
                       cgcgACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCgggg                                                                                                                             <  1:579004/150‑1 (MQ=255)
                       cgcgACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCgggg                                                                                                                             <  2:1617860/150‑1 (MQ=255)
                         cgACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCTCGCTTTTCCACCTCGCCCTTTCCGGGAGAGGGCCGGGGTg                                                                                                                           >  1:269777/1‑150 (MQ=255)
                           aCTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAg                                                                                                                         >  1:1399403/1‑150 (MQ=255)
                           aCTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAg                                                                                                                         <  2:1163078/150‑1 (MQ=255)
                             tGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTACCGAGAAAAATTAATTCTTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTACCCCTCGCCCTTTCAGTGAGAACGCCGGGGTGAggg                                                                                                                       >  1:1248121/1‑150 (MQ=255)
                              gCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGt                                                                                                                      >  1:827220/1‑150 (MQ=255)
                                caccaGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTaa                                                                                                                    <  1:1537190/150‑1 (MQ=255)
                                 accaGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTaaa                                                                                                                   >  2:626842/1‑150 (MQ=255)
                                 aacaGTCACTAGTGCTATATCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGTCTTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTaaa                                                                                                                   <  2:819678/148‑1 (MQ=255)
                                  ccaGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAat                                                                                                                  <  1:368094/150‑1 (MQ=255)
                                    aGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAatat                                                                                                                >  2:36482/1‑150 (MQ=255)
                                         caAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcgcc                                                                                                           >  1:156331/1‑149 (MQ=255)
                                         caAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACATATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcgcg                                                                                                           <  2:372817/150‑1 (MQ=255)
                                           aGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGcc                                                                                                         <  2:1683650/150‑1 (MQ=255)
                                            gTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCa                                                                                                        >  2:1425792/1‑150 (MQ=255)
                                             tGCTTTAGCCATTTGTTTTTCCTAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTGGGCACGCTGTTCCCCCTCGCCCTTTCAGGGAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAg                                                                                                       <  1:1305974/150‑1 (MQ=255)
                                                 ataGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGcc                                                                                                   >  2:1456596/1‑150 (MQ=255)
                                                  taGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGTTAGGGTAAATATTCGCGCCAGTGCCg                                                                                                  <  1:387496/150‑1 (MQ=255)
                                                  taGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGTGTAAATATTCGCGCCAGTGCCg                                                                                                  <  1:1429308/150‑1 (MQ=255)
                                                     ccATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGcc                                                                                               <  1:567399/150‑1 (MQ=255)
                                                     ccAATTGTTTTTCCGAGAAAAATTAATTCGGTGACGTATCTTTGAGATTGCTAACAACCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGcc                                                                                               <  2:306947/150‑1 (MQ=255)
                                                       aTTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTg                                                                                             >  2:870394/1‑150 (MQ=255)
                                                          tGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTTACTAACTGGCGCATGCTTCATTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTc                                                                                          <  1:1258636/150‑1 (MQ=255)
                                                           gTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTcc                                                                                         >  1:287711/1‑150 (MQ=255)
                                                            tttttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTccc                                                                                        <  1:1630737/150‑1 (MQ=255)
                                                              tttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAGCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                                      <  1:1169911/150‑1 (MQ=255)
                                                                 ccGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAc                                                                                   <  1:616067/150‑1 (MQ=255)
                                                                  cGAGAAAAATTATTTCGTTGACGTATCTTTGAGTTTGCTAACAATCTGGCTCAATCTTCTATGTTGGCACGCTGTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                  <  2:720157/150‑1 (MQ=255)
                                                                  cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                  >  2:561321/1‑150 (MQ=255)
                                                                  cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                  >  1:318772/1‑150 (MQ=255)
                                                                  cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGATCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCTCGCAAGTGCCGGCCTGTTCCCCACAcc                                                                                  >  2:877460/1‑150 (MQ=255)
                                                                  cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTCCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                  <  2:1399403/150‑1 (MQ=255)
                                                                   gagaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAccc                                                                                 <  2:804130/150‑1 (MQ=255)
                                                                   gagaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAccc                                                                                 <  1:1513701/150‑1 (MQ=255)
                                                                     gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCCCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAATTATTCGCGCCCGTGCCGGCCTGTTCCCCTCACCCTa                                                                               >  2:158256/1‑150 (MQ=255)
                                                                     gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTa                                                                               >  1:1658611/1‑150 (MQ=255)
                                                                           ttaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctc                                                                         <  1:601485/150‑1 (MQ=255)
                                                                              attCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                      <  1:381814/150‑1 (MQ=255)
                                                                              attCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                      <  2:514693/150‑1 (MQ=255)
                                                                                 cGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGTCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa                                                                   >  2:232148/1‑150 (MQ=255)
                                                                                   ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                 <  1:1564530/150‑1 (MQ=255)
                                                                                    tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                <  2:210625/150‑1 (MQ=255)
                                                                                        gTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACACTAACCCTCTACCCAAAGGGGc                                                            >  1:842958/1‑150 (MQ=255)
                                                                                         tATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                                           <  2:791479/150‑1 (MQ=255)
                                                                                          aTCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCATTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                          <  1:1327991/150‑1 (MQ=255)
                                                                                             tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAggg                                                       <  1:626842/150‑1 (MQ=37)
                                                                                              ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgggg                                                      <  2:827220/150‑1 (MQ=255)
                                                                                               tGAGATTGCTAACAAACTCGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                     <  1:1581155/150‑1 (MQ=255)
                                                                                                 agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgggttag                                                   <  1:980665/150‑5 (MQ=37)
                                                                                                    ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTc                                                <  2:1571291/150‑1 (MQ=255)
                                                                                                    ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTc                                                <  2:826801/150‑1 (MQ=255)
                                                                                                     tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGCCTGTCt                                               >  2:1448137/1‑150 (MQ=255)
                                                                                                      gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTg                                              >  1:1183645/1‑150 (MQ=255)
                                                                                                      gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGGGAGGGTAAATATTCGCGCCAGTGCCGGCATGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTg                                              <  2:858942/150‑1 (MQ=255)
                                                                                                       cTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGa                                             <  1:1456596/150‑1 (MQ=255)
                                                                                                        tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTTTCTGAg                                            >  2:851496/1‑150 (MQ=255)
                                                                                                         aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                           >  1:633865/1‑150 (MQ=255)
                                                                                                            aaaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGGCTGTCTGAGCACt                                        >  1:218198/1‑150 (MQ=255)
                                                                                                             accTGGCGCAATCTTCAATGTTGGCACGCTGTTCCCCCTCGCCCTTTCCGGGTGAGGGCCGGGGTCTGGGTAAATATTCGCGCCATTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACtt                                       <  1:1272724/148‑1 (MQ=25)
                                                                                                                 ggttCAATCTTCAAGGTGGGCAGGCTTTTCCCCCTGGCCCTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTa                                  <  2:198395/147‑1 (MQ=25)
                                                                                                                    tCAATCTTCAATGTTGGCACGGTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGGCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACt                                <  2:704093/150‑1 (MQ=255)
                                                                                                                         cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGt                                   >  2:252929/1‑142 (MQ=255)
                                                                                                                         cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGt                                   <  1:252929/142‑1 (MQ=255)
                                                                                                                          ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCa                          >  2:1195379/1‑150 (MQ=255)
                                                                                                                             aaTGTTGGCACGCTTTTCCCCCTCACCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCCGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCt                       <  1:1375787/150‑1 (MQ=255)
                                                                                                                              aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGGGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                      <  1:764522/150‑1 (MQ=255)
                                                                                                                              aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                      <  2:103711/150‑1 (MQ=255)
                                                                                                                              aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                      <  2:588361/150‑1 (MQ=255)
                                                                                                                              aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                      <  2:719102/150‑1 (MQ=255)
                                                                                                                                gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCTTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAc                    <  1:1669821/150‑1 (MQ=255)
                                                                                                                                    gCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaa                >  2:926482/1‑150 (MQ=255)
                                                                                                                                       cGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAg             >  2:473256/1‑150 (MQ=255)
                                                                                                                                           cccggcccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGCGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc         >  2:1389466/6‑150 (MQ=25)
                                                                                                                                            ccggcccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGcc        >  2:1167512/5‑150 (MQ=25)
                                                                                                                                                cccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg    >  2:1547790/1‑150 (MQ=25)
                                                                                                                                                cccTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAAAATTCGCGCCAGAGCCGGCCTGTTCCCCTCACCCTAACCCTCACCCCACAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg    <  2:1055312/150‑1 (MQ=37)
                                                                                                                                                 ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc   >  1:1103176/1‑150 (MQ=255)
                                                                                                                                                  cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGTGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg  <  2:44417/150‑1 (MQ=255)
                                                                                                                                                  cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg  >  2:1161476/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TAATGCAGTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCG  >  NC_000913/3045845‑3046140

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: