Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,510,364 C→T 10.4% V205I (GTT→ATT) ‡ yhfX ← putative PLP‑binding protein YhfX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,510,3640CT10.4% 159.0 / 7.0 58V205I (GTT→ATT) ‡yhfXputative PLP‑binding protein YhfX
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base T (4/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGGCATGAGTCACACCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAATC  >  NC_000913/3510215‑3510508
                                                                                                                                                     |                                                                                                                                                
ttCGGCATGAGTCACACCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAAc                                                                                                                                                  <  2:538994/150‑1 (MQ=255)
ttCGGCATGAGTCACACCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAAc                                                                                                                                                  <  2:1024163/150‑1 (MQ=255)
            cacaCCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGAAAAGTTTTGCCGACacc                                                                                                                                      <  1:977974/150‑3 (MQ=255)
              cacCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCt                                                                                                                                    <  2:1296969/150‑1 (MQ=255)
              cacCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCAGAAAGcca                                                                                                                                    >  2:174213/1‑149 (MQ=255)
               acCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGAAAAGTTTTGCCGACaccgaa                                                                                                                                   <  1:455188/150‑6 (MQ=255)
                 cGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCAt                                                                                                                                 <  1:1224503/150‑1 (MQ=255)
                  gTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATc                                                                                                                                >  1:480490/1‑150 (MQ=255)
                    aTTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATccc                                                                                                                              <  1:508107/150‑1 (MQ=255)
                       gcgcCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCaaa                                                                                                                           >  2:1342081/1‑150 (MQ=255)
                       gcgcCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCaaa                                                                                                                           <  1:967452/150‑1 (MQ=255)
                         gcCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGc                                                                                                                         <  1:242638/150‑1 (MQ=255)
                           cAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAg                                                                                                                       >  1:596887/1‑150 (MQ=255)
                            aGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCATGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAgg                                                                                                                      <  1:439115/150‑1 (MQ=255)
                              tAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCGGGCa                                                                                                                    <  1:240978/150‑1 (MQ=255)
                                   gCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGa                                                                                                               >  1:1298170/1‑150 (MQ=255)
                                     aGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACCGCCTCATCCCAAAGCAGGCAAGGGaaa                                                                                                             >  1:1121184/1‑150 (MQ=255)
                                                  tGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCgg                                                                                                >  1:228620/1‑150 (MQ=255)
                                                    gTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCt                                                                                              <  2:1298170/150‑1 (MQ=255)
                                                     tCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCCGTGTCGCCAAAACTATTAAAACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTa                                                                                             >  2:1218848/1‑150 (MQ=255)
                                                            aGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAg                                                                                      <  1:861636/150‑1 (MQ=255)
                                                                      aGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGa                                                                            >  2:1148383/1‑150 (MQ=255)
                                                                            tCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTgg                                                                      <  1:954937/150‑1 (MQ=255)
                                                                             cAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGa                                                                     <  1:900300/150‑1 (MQ=255)
                                                                             cAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTACGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGa                                                                     >  2:642040/1‑150 (MQ=255)
                                                                               aGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGAtt                                                                   <  1:1342081/150‑1 (MQ=255)
                                                                                 tGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTc                                                                 <  1:897160/150‑1 (MQ=255)
                                                                                    aaTACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGc                                                              <  2:272218/150‑1 (MQ=255)
                                                                                     aTACCAGATTTCGCCAGTTGATCCCGTGTCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCg                                                             <  2:315903/150‑1 (MQ=255)
                                                                                     aTACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTAGGTTTCTGCAAAACTATTAAAACATACTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCg                                                             >  2:966317/1‑150 (MQ=255)
                                                                                       aCCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGAc                                                           <  1:323648/150‑1 (MQ=255)
                                                                                        ccAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATAACGTTTCTCCACAAATAATAAAAAAAAATAATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACa                                                          >  2:408625/1‑150 (MQ=255)
                                                                                                cGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCgg                                                  >  1:1330294/1‑150 (MQ=255)
                                                                                                  ccAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCa                                                <  1:994337/150‑1 (MQ=255)
                                                                                                      ttGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAacac                                            <  1:1216408/150‑1 (MQ=255)
                                                                                                       tGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACAcc                                           >  1:25234/1‑150 (MQ=255)
                                                                                                        gATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCt                                          >  1:743707/1‑150 (MQ=255)
                                                                                                        gATCCCGTGCCTGTATCAGCGTGTGAAGATTAGGTGTCGGCAAAACTTTTCAGAAAGCCTCATCCCAAAGCAGGTAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCt                                          >  2:550153/1‑150 (MQ=255)
                                                                                                        gATCCCGTGCCTGTATCAGCGTGTGAAGATTACGTGTCTGCAAAACTATTCCGAAAGCATCATCCCAAAGCAGGTAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCt                                          >  2:554934/1‑150 (MQ=255)
                                                                                                          tCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTg                                        <  2:645991/150‑1 (MQ=255)
                                                                                                            ccGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAg                                      >  1:672758/1‑150 (MQ=255)
                                                                                                              gTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAggg                                    <  1:66456/150‑1 (MQ=255)
                                                                                                              gTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAggg                                    <  1:1275162/150‑1 (MQ=255)
                                                                                                                      aTCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAACAGGCAAGGGAAATGGGTAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGc                           >  1:888658/1‑150 (MQ=255)
                                                                                                                         aGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGctc                         <  1:1148383/150‑1 (MQ=255)
                                                                                                                         aGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGctc                         <  2:596887/150‑1 (MQ=255)
                                                                                                                             tgtgAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTGGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCt                     >  1:1290671/1‑150 (MQ=255)
                                                                                                                             tgtgAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCt                     <  1:800483/150‑1 (MQ=255)
                                                                                                                               tgAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTgg                   >  1:266585/1‑150 (MQ=255)
                                                                                                                                gAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGc                  >  1:330923/1‑150 (MQ=255)
                                                                                                                                    aTTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCgg              <  2:672758/150‑1 (MQ=255)
                                                                                                                                     gaggctgtCGGAAAAGTTTTGCCGACAGTCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGa             >  2:937716/6‑150 (MQ=255)
                                                                                                                                       cGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATa           >  2:285142/1‑150 (MQ=255)
                                                                                                                                        gctgtCGGAAAAGTTTTGCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATaa          >  1:186165/3‑150 (MQ=255)
                                                                                                                                        gctgtCGGAAAAGTTTTGCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATaa          >  2:829643/3‑150 (MQ=255)
                                                                                                                                        gctgtCGGAAAAGTTTTGCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAAACTCCGGCAACACCTTCAGGGCAAAAACGCTCGCCGGGCCCGGAGaa          >  2:1091796/3‑150 (MQ=255)
                                                                                                                                         gtgtCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATaaa         <  1:642040/150‑1 (MQ=255)
                                                                                                                                         gtgtCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATaaa         <  1:608542/150‑1 (MQ=255)
                                                                                                                                               gCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAat   >  1:1149877/1‑150 (MQ=255)
                                                                                                                                                cAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAatc  >  2:981308/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TTCGGCATGAGTCACACCGTATTGCGCCAGTAATGGCAGCGAAGTGCAGCTGGTCGCTGAAGGCGCGTTCAGTTGCTCAAGTGCAATACCAGATTTCGCCAGTTGATCCCGTGCCTGTATCAGCGTGTGAAGATTCGGTGTCGGCAAAACTTTTCCGACAGCCTCATCCCAAAGCAGGCAAGGGAAATGGGTAAGTCCGGCTAAATGCAGCCCTGGCAGATTCTGGATTTCGGCGACAATCTCCGGCAACACCTTGAGGGCAAAACCGCTCTCCTGGCCCGGATAAAGAAAATC  >  NC_000913/3510215‑3510508

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