Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,379 C→A 6.5% F620L (TTC→TTA) ‡ gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,3790CA 6.5% 274.2 / 7.3 93F620L (TTC→TTA) ‡gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (44/43);  new base A (2/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCA  >  NC_000913/869230‑869528
                                                                                                                                                     |                                                                                                                                                     
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:556090/1‑150 (MQ=255)
                                                                                           cGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATg                                                            <  2:1465819/150‑1 (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:1594562/150‑1 (MQ=255)
                                                                                                              ggggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTg                                         <  2:763107/150‑1 (MQ=255)
                                                                                                              ggggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTg                                         <  2:1547254/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)
                                                                                                                                                     cATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGGGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTcacca  <  2:1574838/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGATCTTCCCAGCAATATTCATCTGCGTGGCGAAGAGGTGGCAGCCGTCTCGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCA  >  NC_000913/869230‑869528

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: