Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,751,5370AC47.8% 82.2 / 15.1 67V111G (GTT→GGT) yiaTputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (24/11);  new base C (0/32);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 4.40e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGG  >  NC_000913/3751390‑3751684
                                                                                                                                                       |                                                                                                                                                    
aCGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAAccc                                                                                                                                                  >  2:541376/1‑151 (MQ=255)
  gCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAAccccc                                                                                                                                                >  1:946857/1‑151 (MQ=255)
   cccAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCg                                                                                                                                               <  1:608495/151‑1 (MQ=255)
   cccAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCg                                                                                                                                               >  2:709675/1‑151 (MQ=255)
      aGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGcc                                                                                                                                             >  2:26164/1‑150 (MQ=255)
        cgcAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCAATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCAt                                                                                                                                           >  1:396763/1‑150 (MQ=255)
               gtCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTAc                                                                                                                                   >  1:170323/1‑151 (MQ=255)
                      aCAGACGCCCTTTCTGCTTTTGGTGGAATCCCGAACGCCCCCCCCCCCGGCCGTTGCTGTTACCGAGACCGCCCACGGCTGCAGCGCTTTTCACGGTCCCCCACCGCGCGTGTTGACACCGGGCACCCCCCGCCATAGCCGTACTATCAc                                                                                                                             <  1:769161/150‑1 (MQ=255)
                      aCAGACGCCCTACCTGCATTTTGGGAAATACCGATGGCCCCCCCCCCCAGCCTTTGCTTTTATCCAGACCGTCCGCAGCGGCAGAGGCTTTCCCGGTTCCCCACCGCTCGGGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCAc                                                                                                                             <  1:523222/150‑1 (MQ=255)
                       cagagacccTATCTTCATTTTGGGTAATCCCGATATCCCCCCCCCCCGCCCTTTCCTCTTTCCCAGACCCTCCGCGCCGGCAGAGGCTTTCCGGTTCCCCCTCCCCCGTTGGTGACACCAGACACCCCCCGCCATAGCCGTACTATCTCg                                                                                                                            <  1:1054673/143‑1 (MQ=255)
                          acccccTACCGGCATTTTGTGAAATCCGGTAACCCCCCCCACCCAGCGTTTTCTTTTATCCAAATGGTCGGCAGCTCGAGAGCTTTCCCCGGTTCCCCAACGCTCGTGTGTATGCCGGCCAACCCCCGCCATAGCCGTACTATCACGCtt                                                                                                                         <  2:170901/146‑1 (MQ=255)
                           cgacCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTAt                                                                                                                       >  1:237512/1‑151 (MQ=255)
                             acCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATc                                                                                                                      >  1:5972/1‑150 (MQ=255)
                              ccTTTCTCCTTTTTGTGAAATCCCGAAGGCCCCCCCCCCCCCCCGTTGCTGTTACCAAAACCGCCCGCGGCGGCGGAGCTTTTCCCGTTCCCCCACCCCTGGTGGTAATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCa                                                                                                                     <  1:1076334/150‑1 (MQ=255)
                                  tCTGCATTTTGTGAAATCCCGATAGCCCCCCCCCCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGc                                                                                                                  <  2:415086/149‑1 (MQ=255)
                                    tGCTTTTTGTGAAATCCCGATAGCCCCCCCCCCCAGCCTTTGCTTTTACCCAGACCGCCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTAATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATc                                                                                                                      <  1:431252/143‑1 (MQ=255)
                                       aTTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGc                                                                                                            >  2:1206913/1‑150 (MQ=255)
                                        ttttGTGGAATACGGAAACCCCCCCCCCCCAGCCGTTGCTGTTACCCACACCGTCCGCAGCTGGAGAGCCTTTCACGTTTCCCCAACACTCGTGGGGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGc                                                                                                            <  2:15793/149‑1 (MQ=255)
                                         tttGTGAAATCCCGACAGCCCCCCCCCCCAGCCGTTCCTGTTATCCAAACCGCCCGCGCCGGCAGAGGTTTTCAGGGTCCCCCAACCCTCGTGGTAATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGt                                                                                                               <  1:996650/145‑1 (MQ=255)
                                              ggAATCCCGATACCCCCCCCCCCCGGCCGTTGCTGTTACCCAGAGCGTCCGCAGCTGCAGGGGCTTTCACGGTCCCCCAACGCTCGTGGTGATCCCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCAt                                                                                                     <  2:828554/150‑1 (MQ=255)
                                                aaTACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAAt                                                                                                  >  1:1075744/1‑151 (MQ=255)
                                                 aTACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATc                                                                                                 >  2:677266/1‑151 (MQ=255)
                                                 aTACCAATACCCCCCCCCCCCAGCCGTTGCTGTTACCCAGACCGTCCGCAGCTGCAGAGGCTTTCACGTTTCCCCACCCCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAAt                                                                                                  <  2:1014451/150‑1 (MQ=255)
                                                    ccgatgcccccccccccCAGCCGTCGCTGTTATCCAGAACGTCCGCAGCTGCGGTGCTTTTCCCGGTTCCCCACCGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCg                                                                                                <  2:508338/140‑1 (MQ=255)
                                                    ccGATAGCCCCCCCCCCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTCCCCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTc                                                                                              <  2:772263/151‑1 (MQ=255)
                                                    ccGATACCCCCCCCCCCCAGCCGTTGCTGTTACCCAGACCGCCCGCAGCTGCAGGGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTc                                                                                              <  1:233976/151‑1 (MQ=255)
                                                    ccGAGAGCTCCCCCCCCCACCCGTTGCTGTTACCCAGAACGCCCGCGGCTGCAGGGGCTTTCACGGTCCCCCACCGCTCGTGGTGATACCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTc                                                                                              <  1:476764/151‑1 (MQ=255)
                                                    cagATACCCCCCCCCCCCACCCGTCGCTGTTATCCAGATCGTCCGCACCTGCGGGGGCTTCCACGGTCCCCCACCGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGt                                                                                               <  2:1083719/148‑1 (MQ=255)
                                                       atacccccccccccCACCCGTTGCTGTTACCCAGAACGTCCGCGGCGGCAGAGGTTTCCACGTTCCCCCACCGCTCGTGGTGATCCCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCAtt                                                                                           <  1:1137341/145‑1 (MQ=255)
                                                             ccccccccCAGCCGTTGCTGTTATCCAGACCGTCCGCAGCTGCAGAGCCTTTCCCGGTCCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                     <  2:350345/151‑1 (MQ=255)
                                                              ccccACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                     >  2:380293/1‑150 (MQ=255)
                                                                 cccccaCCCGTTCCTGTTACCAAGACCGGCCCCAGCGGCAGGGGCTTTCGCGGTCCCCAACCGCTGGTGTTGACCCCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAggg                                                                                 <  1:376255/149‑1 (MQ=255)
                                                                   cccaGCCGTTGTTGTTATCCAGACCTTCCGCAGCTGCAGAGGCTTTCACGGTCCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGt                                                                                <  2:552111/150‑1 (MQ=255)
                                                                        cctttGCTGTTATCCAGAACGTCCGCAGCTGCAGGGGCTTTCACGTTTCCCCAACGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAAtt                                                                          <  2:1142754/148‑1 (MQ=255)
                                                                        ccGTGGCTGTTGTCAAAATCGTCCGCGGGGGCGGAGGCTTTCCGGGTCCCCCAACGCTGGTGGGGATACCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAAtt                                                                          <  1:241452/151‑1 (MQ=255)
                                                                           ttGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTcc                                                                        >  2:1006868/1‑150 (MQ=255)
                                                                           ttGCTGTTACCCAGACCGCCCGCAGCTGCAGAGGCTTTCACGGTCCCCCAACGCTCGTGGTGATCCCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTc                                                                         <  1:162085/149‑1 (MQ=255)
                                                                            tGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTcc                                                                        >  2:929231/1‑149 (MQ=255)
                                                                             tcttttATTCAGAACCGTCCCGAGCGGAGGAGTCTTTCACGTTCCCCCACGTCCCGTGGTAATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTc                                                                         <  1:351228/144‑1 (MQ=255)
                                                                               tGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTGTCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCgg                                                                   >  1:37002/1‑151 (MQ=255)
                                                                                      ccaaacctcccccAGCGGCAGAGGTTTTCCCGTTCCCCCAACCCTCGTGGTGACCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCa                                                             <  1:588120/139‑1 (MQ=255)
                                                                                       aGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAg                                                            <  2:780641/150‑1 (MQ=255)
                                                                                       aGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAgg                                                           >  1:976319/1‑151 (MQ=255)
                                                                                            cGTCCGCAGCTGCGGAGGCTTTCGCGTTCCCCCACCGCTCGTGGTGATCCCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGAt                                                         <  1:8610/149‑1 (MQ=255)
                                                                                              tCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCgg                                                     >  1:755002/1‑151 (MQ=255)
                                                                                               ccccggCTGCGGAGCTTTTCACGGTCCCCCAACGCTCGTGGTGATACCGGCAACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCg                                                      <  2:286841/144‑1 (MQ=255)
                                                                                                cGCAGCTGCGGAGGTTTCCACGGTCCCCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTc                                                   <  1:470293/151‑1 (MQ=255)
                                                                                                cGCAGCTCCGGAGGCTTTCCCGTTCCCCCAACGCCCGTGTGGATACGAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATg                                                        <  2:993035/146‑1 (MQ=255)
                                                                                                     cTGCAGAGTCTTTCACGTTCCCCCACCGCTCGTGGTGATACCAGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAg                                                 <  1:26163/148‑1 (MQ=255)
                                                                                                        cgggggCTTCCACGGTCCCCCAACGCTCGTGGTGATCCCGGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCa                                                  <  2:698546/140‑1 (MQ=255)
                                                                                                          ggggCTTCCCCGCTCCCCCAGCGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTGTCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGtt                                         <  2:652423/149‑1 (MQ=255)
                                                                                                               tttCCCGTTCCCCCACCGCTCGTGGGGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAtt                                     <  1:576080/150‑1 (MQ=255)
                                                                                                                tctcccgTTCCCCAATCGCGCGTGGAGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTc                                        <  2:1068548/141‑1 (MQ=255)
                                                                                                                  cgcGGTTCCCCAACGCTGGTGGTGATACCAGACACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAtttt                                   <  1:541350/147‑1 (MQ=255)
                                                                                                                    cGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTc                               >  1:208664/1‑151 (MQ=255)
                                                                                                                     ggTCCCCCAACGCGGGGGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCAttt                                    <  2:970526/145‑1 (MQ=255)
                                                                                                                       tcccccACCGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCa                              <  2:170327/147‑1 (MQ=255)
                                                                                                                          cccACCGCTCGGGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTg                         <  1:84358/151‑1 (MQ=255)
                                                                                                                           ccAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGc                        <  2:367522/151‑1 (MQ=255)
                                                                                                                                 gCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGAc                     >  2:271010/1‑149 (MQ=255)
                                                                                                                                 gCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGAc                     >  2:1203475/1‑149 (MQ=255)
                                                                                                                                         ggtgATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGaa                 <  2:557990/147‑1 (MQ=255)
                                                                                                                                             aTACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCt         >  2:4420/1‑151 (MQ=255)
                                                                                                                                               aCCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTa        >  1:678504/1‑150 (MQ=255)
                                                                                                                                                ccAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAAc      >  2:565804/1‑151 (MQ=255)
                                                                                                                                                   ggCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGt    <  2:237518/150‑1 (MQ=255)
                                                                                                                                                    gCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTgg  >  2:415004/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTCCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGG  >  NC_000913/3751390‑3751684

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