Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,438,800 T→C 7.5% D662G (GAC→GGC)  pfo ← pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,438,8000TC 7.5% 259.8 / 9.7 95D662G (GAC→GGC) pfopyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (48/38);  new base C (3/4);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTT  >  NC_000913/1438653‑1438944
                                                                                                                                                   |                                                                                                                                                
ggcggcACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCt                                                                                                                                                >  1:574639/1‑150 (MQ=255)
ggcggcACCACTTTTGCGCGAATAGCTGAGAGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCt                                                                                                                                                >  1:1096935/1‑150 (MQ=255)
    gcACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTggcg                                                                                                                                            <  1:655066/150‑1 (MQ=255)
        cacTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAg                                                                                                                                        >  1:1677112/1‑150 (MQ=255)
        cacTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAg                                                                                                                                        >  1:844052/1‑150 (MQ=255)
         acTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGCCAGACGGCacc                                                                                                                                       <  1:1503583/150‑3 (MQ=255)
             ttGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGCCAGACGGCacctggc                                                                                                                                   <  1:1266940/150‑7 (MQ=255)
             ttGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGCCAGACGGCacctggc                                                                                                                                   <  1:1266935/150‑7 (MQ=255)
             ttGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGGTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGCCAGACGGCacctggc                                                                                                                                   <  2:72972/150‑7 (MQ=255)
               gcgcGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGaa                                                                                                                                 >  1:765307/1‑150 (MQ=255)
                ggcctaTAGCTGAGTGAGGGCAAGCGGGAACGCAGTGGTTAGATTTGGTACAGAGTTCCTCATTGCAGATGGGGATGTCTTCGGTGATATTGCGTTTATCCCAGCTCTTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGaaa                                                                                                                                <  2:822070/145‑1 (MQ=255)
                cgcgAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGaaa                                                                                                                                >  2:1117684/1‑150 (MQ=255)
                cgcgAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTATCCAGATGGGGATCTCTTCGGCGATATTCCGTATTTCCCAGCGCGTAGTGCCCATCGGCCACGAGCCGTCTGGCGGCAGCGCCGaaa                                                                                                                                >  1:445746/1‑150 (MQ=255)
                  cgAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGGGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGAGGCAGCGCCGAAACg                                                                                                                              >  1:1072976/1‑150 (MQ=255)
                     aTAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACgggg                                                                                                                           <  2:361097/150‑1 (MQ=255)
                      tAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGa                                                                                                                          <  1:1297765/150‑1 (MQ=255)
                       aGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAg                                                                                                                         <  1:1554650/150‑1 (MQ=255)
                        gCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAgg                                                                                                                        <  2:753004/150‑1 (MQ=255)
                        gCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAgg                                                                                                                        >  2:122431/1‑150 (MQ=255)
                        gCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAgg                                                                                                                        <  1:10813/150‑1 (MQ=255)
                        gCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCGCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAgg                                                                                                                        <  1:399917/150‑1 (MQ=255)
                         cTGAGGGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGAGCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAggg                                                                                                                       <  1:839722/150‑1 (MQ=255)
                          tGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGc                                                                                                                      <  1:1589780/150‑1 (MQ=255)
                          tGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGc                                                                                                                      <  1:1481247/150‑1 (MQ=255)
                              tgcgggCAAGCGGCAACGCAGTGGGTCCATTGGGTACAGCGTTCCTCTATCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCa                                                                                                                  <  2:168812/147‑1 (MQ=255)
                                tgGGCAAGCGGCAACGCCGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGTGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCAcc                                                                                                                <  2:45952/150‑1 (MQ=255)
                                 gggCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCa                                                                                                               <  1:1049617/150‑1 (MQ=255)
                                    cAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGc                                                                                                            <  1:1229097/150‑1 (MQ=255)
                                    cAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGc                                                                                                            >  1:646422/1‑150 (MQ=255)
                                      aGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAACCGGGTAGGGCGTCACCAAGccc                                                                                                          >  2:854760/1‑150 (MQ=255)
                                       gCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCg                                                                                                         <  1:1117684/150‑1 (MQ=255)
                                        cGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCTCGTCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCgg                                                                                                        >  2:774493/1‑150 (MQ=255)
                                        cGGCAACGCAGTGGTTACATTGGGGACCGAGTTCCGCTATCCCGATGGGGATCTCTTCGGCGCTATTGCGTTTTTCCCAGCGCGTACTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCgg                                                                                                        <  1:163791/150‑1 (MQ=255)
                                         ggCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGAGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGc                                                                                                       >  1:672894/1‑150 (MQ=255)
                                              cGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCa                                                                                                  >  1:1455801/1‑150 (MQ=255)
                                                cAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATc                                                                                                >  2:639653/1‑150 (MQ=255)
                                                    ggTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAACCGGGGAGGGCGTCCCCAAGCCAGGCGcgcagcgcag                                                                                            >  2:1498658/1‑150 (MQ=255)
                                                     gTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGa                                                                                                     <  1:927798/140‑1 (MQ=255)
                                                     gTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGa                                                                                                     >  2:927798/1‑140 (MQ=255)
                                                       tACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCgg                                                                                         >  2:1178787/1‑150 (MQ=255)
                                                              ggTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGt                                                                                  <  2:574639/150‑1 (MQ=255)
                                                                tACAGAGTTCCTCTTACCAGATCGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCAGCAAGCCCGGCGAGCGTCGCAGCCGTTACGGtgt                                                                                >  1:92826/1‑148 (MQ=255)
                                                                  cAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTc                                                                              >  1:1455698/1‑150 (MQ=255)
                                                                  cAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTc                                                                              >  2:1557015/1‑150 (MQ=255)
                                                                   agcgTTCCTCTTTCCCGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCa                                                                             <  1:893570/147‑1 (MQ=255)
                                                                   agagTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGGGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCa                                                                             <  2:1544637/150‑1 (MQ=255)
                                                                     agTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACg                                                                           <  1:122431/150‑1 (MQ=255)
                                                                        tCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGaaa                                                                        >  1:1423844/1‑150 (MQ=255)
                                                                         ccTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGTGCCGAAACGGGGAGGGCGTCACCAAGCCCGGTGAGCATCGCAGCGGTTACGGTTTTCACGAAAt                                                                       >  2:1298032/1‑150 (MQ=255)
                                                                           tctTTCCAGATGGGGATCTCTTCGGCGATATTGCGCTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCa                                                                     <  1:1348272/150‑1 (MQ=255)
                                                                             tttCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGGTTTCACGAAATCAgg                                                                   >  1:1707519/1‑150 (MQ=255)
                                                                               tCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGAGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAgggg                                                                 <  1:1669911/150‑1 (MQ=255)
                                                                                 cAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGggcg                                                               >  2:901671/1‑150 (MQ=255)
                                                                                  aGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGggcgg                                                              >  1:430263/1‑150 (MQ=255)
                                                                                  aGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGggcgg                                                              <  2:1455801/150‑1 (MQ=255)
                                                                                       gggATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCCTCGCAGCGGGTACGGTTTTCCCGAAATCAGGGGCGGCCTCg                                                         >  1:174147/1‑150 (MQ=255)
                                                                                       gggATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCg                                                         <  2:1096935/150‑1 (MQ=255)
                                                                                        ggATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCgg                                                        >  2:1353885/1‑150 (MQ=255)
                                                                                        ggATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCgg                                                        >  2:680145/1‑150 (MQ=255)
                                                                                                tCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACt                                                >  1:548972/1‑150 (MQ=255)
                                                                                                  ggCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTgg                                              <  2:765307/150‑1 (MQ=255)
                                                                                                        atTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCg                                        >  2:441572/1‑150 (MQ=255)
                                                                                                         ttGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGa                                       >  2:359566/1‑150 (MQ=255)
                                                                                                          tGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGAt                                      >  2:702185/1‑150 (MQ=255)
                                                                                                           gCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGAtt                                     >  2:1223546/1‑150 (MQ=255)
                                                                                                            cGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGTAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGAATCGGAAACCACTGGCGGTCGATTg                                    >  2:532409/1‑150 (MQ=255)
                                                                                                              tttttCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGgc                                  <  1:901671/150‑1 (MQ=255)
                                                                                                                  tCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCtg                              <  1:20716/150‑1 (MQ=255)
                                                                                                                  tCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCtg                              <  1:724067/150‑1 (MQ=255)
                                                                                                                   cccAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCtgt                             >  2:996365/1‑150 (MQ=255)
                                                                                                                    ccAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCtgtg                            >  1:1460615/1‑150 (MQ=255)
                                                                                                                    ccAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCtgtg                            >  1:517007/1‑150 (MQ=255)
                                                                                                                     cAGCGCGTAGTGCCCATCGGCCAGGTGCCGGCTGGCGGCAGCGCCGAAACGGGGAGGGCGCCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTATTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGGCGATTGGCGCTGTGc                           <  1:761691/150‑1 (MQ=255)
                                                                                                                      aGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCg                          >  2:1019036/1‑150 (MQ=255)
                                                                                                                      aGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCg                          >  2:46629/1‑150 (MQ=255)
                                                                                                                      aGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGCTTTTCACGAAATCAGGGGCGGCATCCGAAACCACAGGCGGTCGATTGGCGCTTTGCg                          >  2:1093387/1‑150 (MQ=255)
                                                                                                                       gcgcGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCgg                         >  1:1176363/1‑150 (MQ=255)
                                                                                                                         gcgTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGAt                       <  2:11025/150‑1 (MQ=255)
                                                                                                                          cgTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGAtt                      <  2:646422/150‑1 (MQ=255)
                                                                                                                              gTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCAGCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTAcc                  <  1:734903/150‑1 (MQ=255)
                                                                                                                                gCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGGAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCgg                >  1:1715920/1‑150 (MQ=255)
                                                                                                                                gCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCgg                >  1:1296565/1‑150 (MQ=255)
                                                                                                                                 cccATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGt               >  2:1253830/1‑150 (MQ=255)
                                                                                                                                 cacATCGGGCAGGGGCCGTGTGGCGGCAGCTCCGAAACGGGGAGGTCGGCACCAAGCCAGTCGAGCATAGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGTAAACACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGt               <  2:1512114/148‑1 (MQ=255)
                                                                                                                                    aTCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGc            <  2:844052/150‑1 (MQ=255)
                                                                                                                                     tCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCa           >  2:904430/1‑150 (MQ=255)
                                                                                                                                      cGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCgg                         >  2:966871/1‑135 (MQ=255)
                                                                                                                                      cGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCgg                         <  1:966871/135‑1 (MQ=255)
                                                                                                                                        gCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACg        <  2:1296565/150‑1 (MQ=255)
                                                                                                                                         ggcgctGCCGCCAGACGGCAGCGCCGAAACGTGGCGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCCGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACgg       >  2:748171/6‑150 (MQ=255)
                                                                                                                                         ggcgctGCCGCCAGACGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACgg       >  2:1312508/6‑150 (MQ=255)
                                                                                                                                         ggcgctGCCGCCAGACGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACgg       >  1:731793/6‑150 (MQ=255)
                                                                                                                                            ggcgCCGGCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAAc    <  2:1237018/147‑1 (MQ=255)
                                                                                                                                             gTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAAct   >  1:556/1‑150 (MQ=255)
                                                                                                                                              tGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAActt  <  1:639653/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCCAGATGGGGATCTCTTCGGCGATATTGCGTTTTTCCCAGCGCGTAGTGCCCATCGGCCAGGTGCCGTCTGGCGGCAGCGCCGAAACGGGGAGGGCGTCACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTT  >  NC_000913/1438653‑1438944

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