Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,928,868 C→A 5.1% A246S (GCG→TCG)  viaA ← protein ViaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,928,8680CA 5.1% 525.9 / 8.8 152A246S (GCG→TCG) viaAprotein ViaA
Reads supporting (aligned to +/- strand):  ref base C (64/80);  new base A (5/3);  total (69/83)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

AGGCGATAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACT  >  NC_000913/3928720‑3929014
                                                                                                                                                    |                                                                                                                                                  
aGGCGATAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGcc                                                                                                                                                   <  1:900309/150‑1 (MQ=255)
aGGCGATAGGAGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGcc                                                                                                                                                   <  1:1193995/150‑1 (MQ=255)
     aTAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTc                                                                                                                                              <  1:1090149/150‑1 (MQ=255)
      tAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAATCTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcg                                                                                                                                             <  2:1324029/150‑1 (MQ=255)
        ggTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcgcg                                                                                                                                           >  1:553830/1‑150 (MQ=255)
        ggTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcgcg                                                                                                                                           >  1:1072563/1‑150 (MQ=255)
        ggTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTcgcg                                                                                                                                           >  2:1265593/1‑150 (MQ=255)
        ggTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGGCGACGGTTcctg                                                                                                                                           <  1:996818/150‑4 (MQ=255)
          tGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCa                                                                                                                                         <  1:1172760/150‑1 (MQ=255)
           gAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCa                                                                                                                                         <  1:672087/149‑1 (MQ=255)
           gAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCa                                                                                                                                         >  2:672087/1‑149 (MQ=255)
           gAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCAc                                                                                                                                        >  2:529551/1‑150 (MQ=255)
           gAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGGCGACGGTTcctgagc                                                                                                                                        <  2:846389/150‑7 (MQ=255)
            aGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCAcc                                                                                                                                       <  1:1212648/150‑1 (MQ=255)
            aGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCAcc                                                                                                                                       <  1:89166/150‑1 (MQ=255)
               aaCTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATg                                                                                                                                    <  2:1072563/150‑1 (MQ=255)
                aCTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATgg                                                                                                                                   <  2:1164549/150‑1 (MQ=255)
                   gTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGc                                                                                                                                <  1:992443/150‑1 (MQ=255)
                    tttttCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCg                                                                                                                               <  1:681869/150‑1 (MQ=255)
                      tttCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATgg                                                                                                                                   >  2:108849/1‑144 (MQ=255)
                         ccaccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCgcgaac                                                                                                                          <  2:150757/150‑6 (MQ=255)
                         ccaccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAgg                                                                                                                          <  1:936085/150‑1 (MQ=255)
                         ccaccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAgg                                                                                                                          >  1:1090069/1‑150 (MQ=255)
                         ccaccaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAgg                                                                                                                          >  1:416714/1‑150 (MQ=255)
                           accaGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGtt                                                                                                                        <  1:186402/150‑1 (MQ=255)
                             caGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTc                                                                                                                      <  2:517762/150‑1 (MQ=255)
                               gCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCa                                                                                                                    <  1:529551/150‑1 (MQ=255)
                                ccGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCAt                                                                                                                   <  2:1356616/150‑1 (MQ=255)
                                 cGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATc                                                                                                                  >  2:18539/1‑150 (MQ=255)
                                 cGACGCTAAAACTCATCCTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTAGTTGAAGACCATCAACCAGCTCAGCTACCTTCGCCGGTTCGCGCAACATCGTGCGGAACGTTTCCAta                                                                                                                  >  2:1381978/1‑149 (MQ=255)
                                       tAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCa                                                                                                            >  2:76964/1‑150 (MQ=255)
                                       tAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCa                                                                                                            <  2:1235327/150‑1 (MQ=255)
                                                   tCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGt                                                                                                >  1:1083759/1‑150 (MQ=255)
                                                   tCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGt                                                                                                >  1:422104/1‑150 (MQ=255)
                                                   tCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGt                                                                                                >  1:1176731/1‑150 (MQ=255)
                                                       gTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTg                                                                                            >  1:776390/1‑150 (MQ=255)
                                                            gTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTg                                                                                       <  2:329236/150‑1 (MQ=255)
                                                             ttATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTgg                                                                                      >  1:145221/1‑150 (MQ=255)
                                                              tATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGc                                                                                     <  2:145221/150‑1 (MQ=255)
                                                              tATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGc                                                                                     <  2:691172/150‑1 (MQ=255)
                                                              aaTCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGc                                                                                     <  2:1315268/149‑1 (MQ=255)
                                                                tCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCtt                                                                                   >  1:883733/1‑150 (MQ=255)
                                                                 cccTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCAGGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGCTATTGATTTGGCTTc                                                                                  >  1:958025/1‑150 (MQ=255)
                                                                    tAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCg                                                                               >  1:1189156/1‑150 (MQ=255)
                                                                     aGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTg                                                                                                                                 <  1:539043/99‑1 (MQ=255)
                                                                      gtgtCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCgag                                                                             <  2:704018/150‑1 (MQ=255)
                                                                      ctgccGCCAGAACTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCgag                                                                             >  1:368522/5‑150 (MQ=255)
                                                                           gCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGc                                                                        >  1:332330/1‑150 (MQ=255)
                                                                           gCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGc                                                                        >  1:664018/1‑150 (MQ=255)
                                                                            ccAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGcc                                                                       >  2:728116/1‑150 (MQ=255)
                                                                            ccAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGcc                                                                       <  1:754129/150‑1 (MQ=255)
                                                                              aGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGcccc                                                                     <  2:1237205/150‑1 (MQ=255)
                                                                                ttCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGc                                                                                                                                >  1:155560/1‑89 (MQ=255)
                                                                                  cTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgct                                                                 >  2:919805/1‑150 (MQ=255)
                                                                                    ggcggcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTTCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgctgc                                                               >  2:437508/1‑150 (MQ=255)
                                                                                    ggcggcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgctgc                                                               <  1:1395474/150‑1 (MQ=255)
                                                                                     gcggcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgctgct                                                              <  1:1265593/150‑1 (MQ=255)
                                                                                      cggcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGctc                                                             <  1:1050633/150‑1 (MQ=255)
                                                                                        gcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATGAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTg                                                           <  2:842841/150‑1 (MQ=255)
                                                                                        gcAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTg                                                           >  2:965990/1‑150 (MQ=255)
                                                                                          aGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATc                                                                                                                  >  2:67207/1‑93 (MQ=255)
                                                                                          aGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGcc                                                         >  1:617972/1‑150 (MQ=255)
                                                                                           ggAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAg                                                                   >  2:737608/1‑139 (MQ=255)
                                                                                           ggAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAg                                                                   <  1:737608/139‑1 (MQ=255)
                                                                                           ggAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCa                                                        <  1:961068/150‑1 (MQ=255)
                                                                                              gaCGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGc                                                     <  1:562261/150‑1 (MQ=255)
                                                                                                cGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGt                                                   >  2:948846/1‑150 (MQ=255)
                                                                                                  tAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGttt                                                 <  1:1057915/150‑1 (MQ=255)
                                                                                                   aaaaTATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTc                                                >  2:293225/1‑150 (MQ=255)
                                                                                                    aaaTATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCa                                               >  1:218674/1‑150 (MQ=255)
                                                                                                     aaTATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAg                                              <  2:1090069/150‑1 (MQ=255)
                                                                                                        atcatcGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTc                                           <  2:416714/150‑1 (MQ=255)
                                                                                                          catcGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGTCGTTTCAGTTCCg                                         >  1:798301/1‑150 (MQ=255)
                                                                                                            tcGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGc                                       <  2:965231/150‑1 (MQ=255)
                                                                                                               cTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGt                                    >  2:582077/1‑150 (MQ=255)
                                                                                                                 ttgttgCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGGGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCTTTTCTGTTGCGGCTGTTc                                  <  2:1038783/150‑1 (MQ=255)
                                                                                                                     tgCAGACCATCATCCTGCTCAGTAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCGATGTGCGCATCGTTGCGCGGTATTTATTTGGCTTCGCGAGATCGCCCCAGGTGCTCTGCCAGGCGTTTCATTTCCGGCTTTTCGTTa                              <  2:64346/150‑1 (MQ=255)
                                                                                                                     tgCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGtt                                   >  2:338047/1‑145 (MQ=255)
                                                                                                                     tgCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTa                              >  2:860072/1‑150 (MQ=255)
                                                                                                                      gCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTaa                             >  1:312409/1‑150 (MQ=255)
                                                                                                                       cAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAg                            >  2:871357/1‑150 (MQ=255)
                                                                                                                       cAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAg                            <  2:368522/150‑1 (MQ=255)
                                                                                                                        aGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGa                           <  2:1191318/150‑1 (MQ=255)
                                                                                                                          aCCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGaaa                         <  2:553830/150‑1 (MQ=255)
                                                                                                                           ccATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAAt                        >  1:502485/1‑150 (MQ=255)
                                                                                                                               cAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCAc                    <  1:130658/150‑1 (MQ=255)
                                                                                                                                aaCCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCAcc                   <  2:885329/150‑1 (MQ=255)
                                                                                                                                  ccTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGt                 <  2:124816/150‑1 (MQ=255)
                                                                                                                                  ccTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGt                 >  1:381827/1‑150 (MQ=255)
                                                                                                                                    tGCTCAGGAACCGTCGCCTGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGACCGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAt               <  2:139536/150‑1 (MQ=255)
                                                                                                                                    ggtgcgcgAACCGGCGACGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTATGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAt               >  2:1049561/8‑150 (MQ=255)
                                                                                                                                     gCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAtt              >  2:960564/1‑150 (MQ=255)
                                                                                                                                      tgcgcgAACCGGCGACGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAttt             >  1:1213376/6‑150 (MQ=255)
                                                                                                                                        cgcgAACCGGCGACGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTca           >  2:986668/4‑150 (MQ=255)
                                                                                                                                        cAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGTGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTca           >  2:364378/1‑150 (MQ=255)
                                                                                                                                         aGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTcac          <  1:948846/150‑1 (MQ=255)
                                                                                                                                              ccGGCGACGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCa     <  2:1213376/150‑1 (MQ=255)
                                                                                                                                                gTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAAc   <  2:422104/150‑1 (MQ=255)
                                                                                                                                                 tCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACt  <  2:1178139/150‑1 (MQ=255)
                                                                                                                                                 tCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACt  <  2:617972/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGGCGATAGGTGAGCAACTGTTTTTCCACCAGCCGACGGTAAAACTCATACTCCAGTTCCGTTATCCCTAGTGTCGCCAGTTCTGGCGGCAGGAGACGTAAAATATCATCGCTTTGTTGCAGACCATCAACCTGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACT  >  NC_000913/3928720‑3929014

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