Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 381,302 A→T 7.1% intergenic (‑42/‑49) yaiX ← / → insC1 pseudogene, interrupted by IS2A; hexapeptide transferase superfamily; putative transferase/IS2 repressor TnpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913381,3020AT 7.1% 171.4 / 2.7 77intergenic (‑42/‑49)yaiX/insC1pseudogene, interrupted by IS2A; hexapeptide transferase superfamily; putative transferase/IS2 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base T (4/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGA  >  NC_000913/381154‑381405
                                                                                                                                                    |                                                                                                       
ccAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCGTTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTaaa                                                                                                         <  2:535663/149‑1 (MQ=255)
 cAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                                                                        >  1:1548122/1‑149 (MQ=255)
 cAATCGTTGCTTCCGCTTCAATAACCGCAGTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGCCTGGCCCCCTGAATCTCCAGACTCCCAATATGACTTAAAt                                                                                                        >  2:1455387/1‑149 (MQ=255)
  aaTCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCGCCCTGAATCTACAGACAACCAATATCAGTTAAATa                                                                                                       >  1:131906/1‑149 (MQ=255)
   atggTTGCTTCCGCTTCAATAACCGCATTATTAATCTCGGTGGCAAAACCAATTTTTACGCCCTTGCTGATTTTTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATaa                                                                                                      <  2:1341264/146‑1 (MQ=255)
     cggtGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGt                                                                                                    <  2:1338139/146‑1 (MQ=255)
       ttGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACTCATATTTACTAGACTGGCCCCCTGATTCTCCAGACAACCAATAACACTTAAATAAGTGa                                                                                                  >  2:49331/1‑149 (MQ=255)
         gCTTCCTATTCAATAACCGCATTTTTAATCTCGGTGGCTAAACCAATTTTTTCGCCATTGCTGAGTATTGTGCCAGGACGAATAAACGCTTAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                                                                                <  2:1024388/149‑1 (MQ=37)
         gCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                                                                                >  2:1331122/1‑149 (MQ=255)
          cTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                                                                               >  1:657989/1‑149 (MQ=255)
          cTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTAACGCCAATCCTGATTAATGTGCCAGGACGAATAAACGCATAATTACTAGACTGACCACCATAATCACACGACAACCAATATCACTTAAATATGTGATAg                                                                                               >  1:926255/1‑149 (MQ=255)
           ttCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGACAGGACGAATAAACGCATAATTACTAGAATGTCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGt                                                                                              >  1:1559113/1‑149 (MQ=255)
            tCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAACCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAATTAAGTGATAGTc                                                                                             >  2:395702/1‑149 (MQ=255)
            tCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTc                                                                                             >  2:1108051/1‑149 (MQ=255)
            tCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTc                                                                                             >  1:6252/1‑149 (MQ=255)
             ccGCTTCAAAAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTTCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATattgg                                                                                            <  2:252635/149‑5 (MQ=38)
              cgcgTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCTTTGCTGATTATTGTGCCATGACGAATAAACGCATAATTAGTAGACTGGCCCCCTGAATCTCCAGAAAACCAATATCACTTAAATAAGTGATAGTCtt                                                                                           <  2:884573/149‑1 (MQ=37)
                   cAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACACCCAATATCACTTAAATAAGTGATAGTCTTAATAc                                                                                      <  2:1207267/149‑1 (MQ=255)
                     attaCCGCATTTTTAATCTCGGTGGAAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTTCTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                                                                                    <  2:1308341/146‑1 (MQ=255)
                     aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATat                                                                                               <  2:752928/138‑2 (MQ=255)
                     aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATat                                                                                               >  1:752928/1‑137 (MQ=255)
                     aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTTATAGTCTTAATACTa                                                                                    >  1:370881/1‑149 (MQ=255)
                     aTAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTCCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                                                                                    >  1:1484383/1‑149 (MQ=255)
                      ttaCCGCATTTTTAATCTCGGTTGCATAACCATTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAACTAAGTGATAGTCTTAATACTAg                                                                                   <  2:689498/147‑1 (MQ=255)
                        aCCGCATTTTTAATGTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                                                                                 <  2:1346265/149‑1 (MQ=255)
                          cGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                                                                               >  2:243172/1‑149 (MQ=255)
                          cGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                                                                               >  2:1384246/1‑149 (MQ=255)
                            cATTTTTAATCTCTGTGGCAAAACCATTTTTTACGCCATTTCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCAAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                                                                             <  2:392283/149‑1 (MQ=255)
                            cATTTTTAATCTCGGTGGCAACACCAATTTTTACGCCATTGCTGATTATTCTGCCAGGACTAATAAACGCATAATTACTAGACTGTCCCCATGATTCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                                                                             <  2:148402/149‑1 (MQ=25)
                            cATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                                                                             <  1:586715/149‑1 (MQ=255)
                             aTTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                                                                            >  2:1497987/1‑149 (MQ=255)
                               ttttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc                                                                          <  2:329388/149‑1 (MQ=255)
                                tttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACt                                                                         >  2:801464/1‑149 (MQ=255)
                                tttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACATAAATAAGTGATAGTCTTAATACTAGTTTTTAGACt                                                                         <  1:532792/149‑1 (MQ=255)
                                  tAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAg                                                                       <  1:1327/149‑1 (MQ=255)
                                     tctcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCATGACGAATAAACGCATAATTACTAGACTGGCCCCCCGAATCTCCAGACAACCAATATCACTTAAATAAGTGTTAGTCTTAATACTAGTTTTTAGACTAGTCa                                                                    <  1:532624/149‑1 (MQ=255)
                                       tcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATACACGCATAATTACTAGACTGGCCCCATGAATCTCCAGACCACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTATACTAGTCAtt                                                                  >  2:1173525/1‑149 (MQ=255)
                                        cGGGGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                                                                 <  2:355812/149‑1 (MQ=255)
                                         gttGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                <  1:1459171/147‑1 (MQ=255)
                                         ggTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                <  2:268193/149‑1 (MQ=255)
                                          gTGGCAAAACCATTTTTTACGCCATTGATGATTATTGTGCCAGGTCGAATAATCGCATAATTACTAGACTGGGCCGCTGAATATCCAGACAACCAATATCACTTAAATTAGTGATAGTCTTATTACTAGTTTTTAGACTAGTCATTGga                                                               <  2:926255/149‑1 (MQ=21)
                                          gTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGga                                                               <  1:1453692/149‑1 (MQ=255)
                                           tGGCAAAACCAATTTTTACGCCATTGCTGATTATTGGGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag                                                              <  2:1164640/149‑1 (MQ=255)
                                               aaaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACa                                                          >  2:1227140/1‑149 (MQ=255)
                                                aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATAACACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGACAAAAg                                                         >  2:1014982/1‑149 (MQ=255)
                                                 aaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAga                                                        <  1:1384246/149‑1 (MQ=255)
                                                  aCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAGTAGTTTTTAGACTAGTCATAGGAGAACAgat                                                       >  2:1458441/1‑149 (MQ=255)
                                                   ccAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAACCGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCCATATCACTTTATTAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGCACAgatg                                                      >  2:865057/1‑149 (MQ=21)
                                                    cAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATTAGTGATAGTGTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatga                                                     >  1:1498008/1‑149 (MQ=255)
                                                    cAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGATCAgatga                                                     >  2:1149510/1‑149 (MQ=255)
                                                    cAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatga                                                     >  2:692540/1‑149 (MQ=255)
                                                    cAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGACCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTGAGACTAGTCATTGGAGAACAgatga                                                     >  1:1169125/1‑149 (MQ=255)
                                                     aaTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatgat                                                    <  2:66189/149‑1 (MQ=255)
                                                       tttttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgattg                                                  <  1:582180/149‑1 (MQ=255)
                                                        ttttACGCGTTTGCGGATTATTGTGCCAGGACGATTAAACGCATAATTCCTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTGTTAATACTAGTTTTTAGACTAGTCATTGGCGAACAGAtgattgt                                                 <  2:600232/149‑2 (MQ=17)
                                                         tttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgattgat                                                <  2:6252/149‑1 (MQ=255)
                                                          ttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATg                                               >  1:724725/1‑149 (MQ=255)
                                                               ccATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTa                                          <  1:607377/149‑1 (MQ=255)
                                                                  ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAggg                                       <  2:533669/149‑1 (MQ=255)
                                                                    gCTGATTATAGTGCCAGGACGAATACACGCATAATTACTAGACTGGCCCACTTAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTACTACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGcc                                     <  1:1024850/149‑1 (MQ=17)
                                                                     cTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCg                                    >  2:200646/1‑149 (MQ=255)
                                                                     cTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGGCTTAGGGCCg                                    <  1:99240/149‑1 (MQ=255)
                                                                       gATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGga                                  >  1:1134958/1‑149 (MQ=255)
                                                                        attattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGCCAACCAATATCACTTAAATAAGTGATAGTCTTAATAGTAGTTTTTAGACGAGTCATTGGAAACAAGATGATAGATGTATTAGGGCCGgat                                 >  2:1389111/1‑148 (MQ=17)
                                                                            ttGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAAc                             >  1:142023/1‑149 (MQ=255)
                                                                                 ccAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGa                        >  1:1367848/1‑149 (MQ=255)
                                                                                 ccAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGa                        >  1:254107/1‑149 (MQ=255)
                                                                                     gACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGc                    <  2:1311444/149‑1 (MQ=255)
                                                                                        gAATAAACGCTTTATTACTAGAGTGGCCCACTGAATCTCCAGACAACCAATATCAGTTAAATAAGTTATAGTCTTAATCCTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCAGGAGAAACGCAGACGGCGTa                 <  2:260791/149‑1 (MQ=11)
                                                                                          aTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTAcc               >  2:1287564/1‑149 (MQ=40)
                                                                                          aTAAACGCATAATTACTAGACTGGCCCCCTGAACCTCCAGACAACCAATATCACTTAACAAAGGGATAGTCTTAATACTAGAGTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTAcc               <  1:300460/149‑1 (MQ=14)
                                                                                              aCGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAATTAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcacac           >  2:523609/1‑149 (MQ=38)
                                                                                              aCGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAGTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTATGGCCGTAGAAAAGCAGACGGCGTAacacac           >  2:335702/1‑149 (MQ=11)
                                                                                                gCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAg         >  1:287173/1‑149 (MQ=35)
                                                                                                  aTAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAAAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAgta       >  1:1032499/1‑147 (MQ=19)
                                                                                                  aTAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGa       >  1:280111/1‑149 (MQ=34)
                                                                                                      ttACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAg   <  1:917315/149‑1 (MQ=31)
                                                                                                       tACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTTATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGCCGGCGTACCACACAGGAAAAGa  >  2:1273192/1‑149 (MQ=11)
                                                                                                                                                    |                                                                                                       
CCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGA  >  NC_000913/381154‑381405

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