Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,308,664 C→T G269D (GGC→GAC)  alsC ← D‑allose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,308,6640CT99.2% 517.7 / ‑4.9 132G269D (GGC→GAC) alsCD‑allose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (68/63);  total (68/64)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 136 total aligned reads displayed.

GACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGGAATTTCA  >  NC_000913/4308516‑4308811
                                                                                                                                                     |                                                                                                                                                   
         gATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAatgatg                                                                                                                                          >  1:514043/1‑151 (MQ=255)
          aTAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTTTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATgg                                                                                                                                         <  1:1210552/151‑1 (MQ=255)
           taattaatCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATgg                                                                                                                                         >  2:1355070/1‑150 (MQ=255)
            aattaacCCGCCCATCACCCCCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTTTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAatga                                                                                                                                            <  1:963199/146‑1 (MQ=255)
             attaatCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGgcgc                                                                                                                                      <  3:412962/151‑1 (MQ=255)
               taatCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTg                                                                                                                                    >  1:1073932/1‑151 (MQ=255)
                 atCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGc                                                                                                                                  >  2:1165297/1‑151 (MQ=255)
                  tCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGc                                                                                                                                  >  4:39945/1‑150 (MQ=255)
                   ccGCCCATCACCACCAGTTGGTAATAGGTTTGTCCCTGCAAAATATTCAGACCGTTGTTGAGGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCaa                                                                                                                                <  2:832859/151‑1 (MQ=255)
                   ccGCCCACCCCCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGTTGGTGCCGATGATCACCCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTTGCat                                                                                                                                <  4:140695/151‑2 (MQ=255)
                        cATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGc                                                                                                                            >  1:845651/1‑150 (MQ=255)
                         aTCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCAt                                                                                                                          >  3:124127/1‑151 (MQ=255)
                          tCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCAt                                                                                                                          >  1:1002635/1‑150 (MQ=255)
                           caccacCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCAt                                                                                                                          >  3:591114/1‑149 (MQ=255)
                           caccacCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCAt                                                                                                                          <  1:242844/149‑1 (MQ=255)
                            accaccAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAgg                                                                                                                       >  3:128868/1‑151 (MQ=255)
                            accaccAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAgg                                                                                                                       >  1:511806/1‑151 (MQ=255)
                             ccaccaGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAgg                                                                                                                       <  2:127571/150‑1 (MQ=255)
                              caccaGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGtt                                                                                                                     >  2:1375179/1‑151 (MQ=255)
                                ccaGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGttt                                                                                                                    <  4:12236/150‑1 (MQ=255)
                                   gtgggTAATAGGTTTGAACCTGAAAAATATCCAGACCGTTGTTGATGGGGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCa                                                                                                                  <  1:897093/146‑1 (MQ=255)
                                         aaTAGGTTTGTCCCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATAAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCaaa                                                                                                                <  1:507096/145‑1 (MQ=255)
                                         aaTAGGTTTGTACCTGAAAAATTTTCAGACCTTTTTTATGGTTGCGGATAATCAACCCGCCAATCTCCACAGAGAAAATGCGCCCCTTCCCGCCGAAGAAACTGCTGCCGCCACTGATGGCGCTGGCAATGGCATAGGTTTCaaatc                                                                                                              <  4:492175/147‑3 (MQ=255)
                                           tAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATa                                                                                                         <  2:1668038/150‑1 (MQ=255)
                                            aGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATAcc                                                                                                       <  2:810457/151‑1 (MQ=255)
                                             ggTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGACCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCg                                                                                                      <  2:1171678/151‑1 (MQ=255)
                                               tttGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCgg                                                                                                     <  2:957454/150‑1 (MQ=255)
                                                 tGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCa                                                                                                   >  2:1614513/1‑150 (MQ=255)
                                                  gTACCTGCAAAATATTCAGACCGTTGTTCAGGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCaa                                                                                                  <  2:1292458/150‑1 (MQ=255)
                                                    aCCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGc                                                                                                >  2:1163012/1‑150 (MQ=255)
                                                    aCCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGGTTCAAAACCCATACCGGCAAGCg                                                                                               >  1:915253/1‑151 (MQ=255)
                                                      cTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGt                                                                                             >  2:1576125/1‑151 (MQ=255)
                                                       tGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGtt                                                                                            >  1:580045/1‑151 (MQ=255)
                                                          aaaaTATTCAGACCGTTTTTGATGTTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCt                                                                                          <  1:1591656/150‑1 (MQ=255)
                                                          aaaaTATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTg                                                                                         >  3:67806/1‑151 (MQ=255)
                                                            aaTATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCg                                                                                       >  2:793664/1‑151 (MQ=255)
                                                               atTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGc                                                                                     >  2:213069/1‑150 (MQ=255)
                                                                ttCAGACCGTTGTTGGTGGTGCCGATGATCACCCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCa                                                                                    <  1:494004/150‑1 (MQ=255)
                                                                ttCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCa                                                                                    <  4:107955/150‑1 (MQ=255)
                                                                  cAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCAcc                                                                                  <  3:15460/150‑1 (MQ=255)
                                                                   aGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCg                                                                                 >  2:1416333/1‑150 (MQ=255)
                                                                    gACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCg                                                                                 >  1:749864/1‑149 (MQ=255)
                                                                    gACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCg                                                                                 >  2:1548394/1‑149 (MQ=255)
                                                                    gACCGTTGTTGATGGGGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGa                                                                                <  1:1405543/150‑1 (MQ=255)
                                                                      ccGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTc                                                                             >  2:912684/1‑151 (MQ=255)
                                                                      ccGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTc                                                                             <  4:388069/151‑1 (MQ=255)
                                                                      ccGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTc                                                                             <  3:275291/151‑1 (MQ=255)
                                                                           gttcATGGTGCGAAGGATCAACCCGCCAAACGCCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGc                                                                         <  1:1572463/150‑1 (MQ=255)
                                                                             tgATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAg                                                                       <  2:1285545/150‑1 (MQ=255)
                                                                             tgATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGt                                                                      >  1:210931/1‑151 (MQ=255)
                                                                                 ggTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTg                                                                    >  2:1639370/1‑149 (MQ=255)
                                                                                  gTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTgaga                                                                 <  3:84105/151‑1 (MQ=255)
                                                                                       cGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAgacgac                                                             >  2:1250116/1‑151 (MQ=255)
                                                                                        gatgatCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGacgacg                                                            >  2:122937/1‑151 (MQ=255)
                                                                                          tgatCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGcc                                                          >  2:141950/1‑151 (MQ=255)
                                                                                          tgatCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGc                                                           <  3:39946/150‑1 (MQ=255)
                                                                                           gatCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCt                                                         <  1:1165162/151‑1 (MQ=255)
                                                                                             tCAGCCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGc                                                           <  1:325262/147‑1 (MQ=255)
                                                                                              cAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCAGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGcc                                                      >  1:1630958/1‑151 (MQ=255)
                                                                                               aaCCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGcct                                                     <  1:680517/151‑2 (MQ=255)
                                                                                               aaCCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGcc                                                      >  2:685167/1‑150 (MQ=255)
                                                                                                 cccGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAg                                                    <  1:1002534/150‑1 (MQ=255)
                                                                                                 cccGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGa                                                   <  3:230466/151‑1 (MQ=255)
                                                                                                  ccGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGAc                                                  >  2:457156/1‑151 (MQ=255)
                                                                                                   cGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGAcc                                                 <  2:944288/151‑1 (MQ=255)
                                                                                                      cAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGc                                              >  1:1073883/1‑151 (MQ=255)
                                                                                                       aaTCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGc                                              >  3:294326/1‑150 (MQ=255)
                                                                                                         tCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGcac                                            >  2:249615/1‑150 (MQ=255)
                                                                                                         tCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGcac                                            >  2:1279794/1‑150 (MQ=255)
                                                                                                         tCACCACAGAGAAAATGCGACCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGAGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGcac                                            >  1:1002684/1‑150 (MQ=255)
                                                                                                          caccacAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGcaca                                           >  1:891645/1‑150 (MQ=255)
                                                                                                            ccacAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACaaa                                         <  2:1364006/150‑1 (MQ=255)
                                                                                                            ccacAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAac                                        <  2:883446/151‑1 (MQ=255)
                                                                                                            ccacAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAac                                        >  2:1102223/1‑151 (MQ=255)
                                                                                                             cacaGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACaaa                                         >  4:124181/1‑149 (MQ=255)
                                                                                                               caGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACAcc                                     >  2:574385/1‑151 (MQ=255)
                                                                                                                agagAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCg                                    <  1:688654/151‑1 (MQ=255)
                                                                                                                 gagaAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCg                                    <  1:1040360/150‑1 (MQ=255)
                                                                                                                     aaaTGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGgaga                                <  2:646719/150‑1 (MQ=255)
                                                                                                                     aaaTGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAgat                               <  1:307560/151‑1 (MQ=255)
                                                                                                                        tGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAgatgat                            >  1:813644/1‑151 (MQ=255)
                                                                                                                        tGCGCCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAgatgat                            >  4:346577/1‑151 (MQ=255)
                                                                                                                          tgcCCCTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCATTGGCATAGGTTTCAAAACCCATACCGGAAACCGGTTCGGCGGCACCGAGTCGTCCGGTTGGGACACCCCCGGCCAGACCGGCACAAACACCGACgatgaaa                           >  1:1266089/2‑148 (MQ=255)
                                                                                                                            ccccGTGCCGCCGAAGAAACTGGTGTCGCCTATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGGGTCGTGCAGTTGATACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATaaa                         <  1:1691176/150‑1 (MQ=255)
                                                                                                                             cccTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAac                        >  1:1097264/1‑150 (MQ=255)
                                                                                                                              ccTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAacac                      >  2:852673/1‑151 (MQ=255)
                                                                                                                              ccTTGCCGCCGAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAacac                      >  2:1525202/1‑151 (MQ=255)
                                                                                                                                   ccgccgAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACaccacc                  >  3:422675/1‑150 (MQ=255)
                                                                                                                                   ccgccgAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACAccacca                 <  3:83840/151‑1 (MQ=255)
                                                                                                                                   ccgccgAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACAccacca                 <  2:1737359/151‑1 (MQ=255)
                                                                                                                                    cgccgAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAg                <  4:520264/151‑1 (MQ=255)
                                                                                                                                    cgccgAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAg                <  1:29596/151‑1 (MQ=255)
                                                                                                                                       tgaagaaACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGAt             >  1:1102562/2‑151 (MQ=255)
                                                                                                                                       cgAAGAAACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGa              >  2:614812/1‑150 (MQ=255)
                                                                                                                                        gaagaaACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGAtg            <  2:500868/151‑1 (MQ=255)
                                                                                                                                         aagaaACTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGAt             >  1:535626/1‑149 (MQ=255)
                                                                                                                                            aaaCTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTgg         <  1:1735849/150‑1 (MQ=255)
                                                                                                                                               cTGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGGAAtt     >  2:1480247/1‑151 (MQ=255)
                                                                                                                                                tGGTGTCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGGAAttt    >  2:1674266/1‑151 (MQ=255)
                                                                                                                                                  gtgtCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGGAATTTCa  <  1:1697035/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GACAGCCGCGATAATTAATCCGCCCATCACCACCAGTTGGTAATAGGTTTGTACCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGGAATTTCA  >  NC_000913/4308516‑4308811

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