Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,378 T→G 6.5% F620C (TTC→TGC) ‡ gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,3780TG 6.5% 272.7 / 6.9 92F620C (TTC→TGC) ‡gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (45/41);  new base G (2/4);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACC  >  NC_000913/869230‑869527
                                                                                                                                                    |                                                                                                                                                     
gAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTc                                                                                                                                                      >  1:1410652/1‑150 (MQ=255)
 aaCCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCa                                                                                                                                                     >  2:818076/1‑150 (MQ=255)
   ccGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGTACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATg                                                                                                                                                   <  1:1293872/150‑1 (MQ=255)
    cGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGc                                                                                                                                                  >  1:1630769/1‑150 (MQ=255)
    cGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGc                                                                                                                                                  >  1:1226341/1‑150 (MQ=255)
    cGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGc                                                                                                                                                  <  1:524744/150‑1 (MQ=255)
     gTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCg                                                                                                                                                 <  1:1349101/150‑1 (MQ=255)
      tCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGt                                                                                                                                                <  1:854129/150‑1 (MQ=255)
        ccGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAg                                                                                                                                              >  1:1547254/1‑150 (MQ=255)
          gACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGAt                                                                                                                                            >  1:463580/1‑150 (MQ=255)
          gACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGAt                                                                                                                                            <  2:1630769/150‑1 (MQ=255)
           aCAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATGAATGCGTattct                                                                                                                                           <  2:1272851/150‑5 (MQ=255)
           aCAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATGAATGCGTattct                                                                                                                                           <  1:13614/150‑5 (MQ=255)
           aCAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATGAATGCGTattct                                                                                                                                           <  2:1515581/150‑5 (MQ=255)
            cAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATaa                                                                                                                                          <  2:716528/150‑1 (MQ=255)
             aaCGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAAc                                                                                                                                         >  2:368986/1‑150 (MQ=255)
              aCGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATGAATGCGTattctgat                                                                                                                                        <  2:577717/150‑8 (MQ=255)
               cGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACAt                                                                                                                                       >  2:889926/1‑150 (MQ=255)
               cGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACAt                                                                                                                                       >  2:1384199/1‑150 (MQ=255)
                gACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACAtt                                                                                                                                      >  1:1326747/1‑150 (MQ=255)
                 aCCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTc                                                                                                                                     >  2:413000/1‑150 (MQ=255)
                 aCCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCCGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTc                                                                                                                                     >  1:1465819/1‑150 (MQ=255)
                    gCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATTCGTCGGGCCTGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAgg                                                                                                                                  <  1:1204281/150‑1 (MQ=255)
                       gCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCgg                                                                                                                               <  2:1114560/150‑1 (MQ=255)
                         gtgtACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGgag                                                                                                                             <  2:1410652/150‑1 (MQ=255)
                         gtgtACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGgag                                                                                                                             <  2:1621272/150‑1 (MQ=255)
                           gtACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCCGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGaa                                                                                                                           >  1:457631/1‑150 (MQ=255)
                           gtACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCCCAACCACAATCAGAATACGCATTCATGCGTAGATCACATTCAGGCGGAGaa                                                                                                                           >  1:488384/1‑150 (MQ=255)
                                  tgTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATg                                                                                                                    >  1:1434827/1‑150 (MQ=255)
                                   gTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATgg                                                                                                                   >  1:763107/1‑150 (MQ=255)
                                     cGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCTTTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCa                                                                                                                 <  1:56501/150‑1 (MQ=255)
                                      ggACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCaa                                                                                                                >  1:616964/1‑150 (MQ=255)
                                      ggACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGAAAACATTCAGGCGGAGAATAAAATGGCaa                                                                                                                >  2:1220098/1‑150 (MQ=255)
                                       gACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACAAACAGGCGGAGAATAAAATGGCAag                                                                                                               >  2:926102/1‑150 (MQ=255)
                                                 ccAGCAATATTCATCTGCGTGGAGAACAGGTGGCAGCCGTATCGTTGCAATGAGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTacac                                                                                                     >  2:963990/1‑150 (MQ=255)
                                                   aGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCg                                                                                                   >  1:1374197/1‑150 (MQ=255)
                                                            aaTCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGtt                                                                                          <  1:795134/149‑1 (MQ=255)
                                                              tCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGCTAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAg                                                                                        >  1:375311/1‑150 (MQ=255)
                                                              tCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAg                                                                                        <  1:1351472/150‑1 (MQ=255)
                                                                   gTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGgcg                                                                                   >  1:693956/1‑150 (MQ=255)
                                                                   gTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGgcg                                                                                   >  1:720598/1‑150 (MQ=255)
                                                                    tGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGgcgc                                                                                  >  2:1337621/1‑150 (MQ=255)
                                                                     ggCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCt                                                                                 <  1:571949/150‑1 (MQ=255)
                                                                            agGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGTTGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCAtt                                                                          <  1:1663685/150‑1 (MQ=255)
                                                                             ggTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTg                                                                         <  2:437951/150‑1 (MQ=255)
                                                                              gTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGc                                                                        <  1:755492/150‑1 (MQ=255)
                                                                              gTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGc                                                                        <  2:463580/150‑1 (MQ=255)
                                                                                    gCCGTGTCGTTGCAATGCGTCGTGCCGGGGCATTACGTCGCACAACCACAATCAGACTACGCCTTCCTGCGTAGATATCATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGc                                                                  <  1:273485/150‑1 (MQ=255)
                                                                                    gCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGc                                                                  <  2:616964/150‑1 (MQ=255)
                                                                                        tctcGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTg                                                              >  2:252600/1‑150 (MQ=255)
                                                                                        tctcGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCAAGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTg                                                              >  2:408636/1‑150 (MQ=255)
                                                                                           cGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCTCAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATg                                                           <  1:1178341/150‑1 (MQ=255)
                                                                                           cGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATg                                                           <  1:1221458/150‑1 (MQ=255)
                                                                                           cGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATg                                                           <  2:1465819/150‑1 (MQ=255)
                                                                                           cGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATg                                                           >  2:556090/1‑150 (MQ=255)
                                                                                           cGTTGCAATGCGTCGGGCCGGGACATTACGTCGCACAACCACACTCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGCGCTGTACACCGGAGTGGGGTAGTGGCGCTGGGCATTGCGACAGCGTTGATg                                                           >  1:333291/1‑150 (MQ=255)
                                                                                                 aaagCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCACTATACGCATTCATGCGTAGATAACATTCAGGCCGACACTAAAATTGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGGGTTGATGGCATCt                                                     <  2:1562435/147‑1 (MQ=255)
                                                                                                 aaTGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCt                                                     >  2:1612867/1‑150 (MQ=255)
                                                                                                     ggtggggCCGGGGCATTACGTCGAACAACCACAATCTGATTACGCATTCATGCTTATATAAGATTCAGGCTGAGAATAAAATGGCAAGAGCTGGACACCGTTGTTGGTTAGTGTCGATGGGCATTGCGACAGCGTTGATGGCATCTtgtg                                                 <  2:1142822/146‑1 (MQ=255)
                                                                                                        cGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACTCAATCATGCGTAGAAAACATTAAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTAAGTGGCGCTGTGCATTGCGACAGCGTTTATGGCATCTTGTGCAt                                              >  2:369970/1‑150 (MQ=255)
                                                                                                           gCCGGGGCATTAGGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGTCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGGGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCg                                           <  2:274418/150‑1 (MQ=255)
                                                                                                             cGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCt                                         <  2:1326747/150‑1 (MQ=255)
                                                                                                              ggggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTg                                        <  2:763107/150‑1 (MQ=255)
                                                                                                              ggggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTg                                        <  2:1547254/150‑1 (MQ=255)
                                                                                                              ggggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTg                                        <  2:1594562/150‑1 (MQ=255)
                                                                                                              ggggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGAATTCGCTg                                        >  2:436897/1‑150 (MQ=255)
                                                                                                                ggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGcc                                      <  2:917038/150‑1 (MQ=255)
                                                                                                                    ttACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAg                                  <  1:546495/150‑1 (MQ=255)
                                                                                                                       cGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGAtg                               >  2:993369/1‑150 (MQ=255)
                                                                                                                        gTCGCACAACCACAATCAGAATACGTATTCATGCGTAGATAACATTCAGGCGGAGATTAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGAtgt                              <  2:1023601/150‑1 (MQ=255)
                                                                                                                         tCGCACAACCACAATCAGAATACGCATTCAAGCGTAGATAACATTCAGGCGGAGAATAAAATGGAAAGAGCTGTACACCGTAGTGGGTTAGTGGCTCTCGGCATTGCGACAGCGTTCATGGCATCTTGTGCATTCGCTGCCAAAGAtgtg                             >  2:850516/1‑150 (MQ=255)
                                                                                                                          cGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATgtgg                            >  2:649145/1‑150 (MQ=255)
                                                                                                                            cacaACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATgtggtg                          <  2:1226341/150‑1 (MQ=255)
                                                                                                                            cacaACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATgtggtg                          >  1:1666049/1‑150 (MQ=255)
                                                                                                                             acaACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGCGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATgtggtgg                         <  1:864450/150‑1 (MQ=255)
                                                                                                                             acaACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATgtggtgg                         <  1:1371334/150‑1 (MQ=255)
                                                                                                                              caACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATgtggtggt                        >  1:1658487/1‑150 (MQ=255)
                                                                                                                                 ccACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGc                     <  1:1002347/150‑1 (MQ=255)
                                                                                                                                  cacaATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCg                    >  2:809986/1‑150 (MQ=255)
                                                                                                                                   acaATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGGGGTGGCgg                   <  2:463784/150‑1 (MQ=255)
                                                                                                                                      aTCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAg                >  2:1313790/1‑150 (MQ=255)
                                                                                                                                        cAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGa              <  1:413000/150‑1 (MQ=255)
                                                                                                                                         aGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGAt             >  1:788783/1‑150 (MQ=255)
                                                                                                                                          gAATACGCATTCATGCGTAGATAACATTCAGGCGTACAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATc            <  2:23039/150‑1 (MQ=255)
                                                                                                                                          gAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACCCCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTCCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGAta            >  1:574473/1‑149 (MQ=255)
                                                                                                                                          gAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATc            >  2:1625260/1‑150 (MQ=255)
                                                                                                                                          gAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATc            >  1:581543/1‑150 (MQ=255)
                                                                                                                                          atctACGCATGAATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATc            >  2:823/4‑150 (MQ=255)
                                                                                                                                            ctACGCATGAATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGa          >  1:578630/2‑150 (MQ=255)
                                                                                                                                             tACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGATGGGCATTTCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGAGGTGGCGGTAGGATCGaa         >  2:1578045/1‑150 (MQ=255)
                                                                                                                                               cGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATCTGGTGGTGGCGGTAGGAACGAAtt       >  2:1608199/1‑150 (MQ=255)
                                                                                                                                                  aTTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTca    <  1:1490704/150‑1 (MQ=255)
                                                                                                                                                   ttCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTcac   <  1:1309676/150‑1 (MQ=255)
                                                                                                                                                    tCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCTGTAGGATCGAATTTcacc  <  1:368986/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACC  >  NC_000913/869230‑869527

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: