Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091340,5150GT42.0% 58.8 / 33.8 69Q473K (CAA→AAA) caiTputative transporter
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base T (29/0);  total (49/20)
Fisher's exact test for biased strand distribution p-value = 1.43e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAA  >  NC_000913/40372‑40664
                                                                                                                                               |                                                                                                                                                      
taaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCgg                                                                                                                                                 >  2:1115487/1‑150 (MQ=255)
 aaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGCTTGTATCgggt                                                                                                                                               >  2:11829/1‑149 (MQ=255)
   attaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGttt                                                                                                                                              >  2:306897/1‑150 (MQ=255)
       aaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTGATCGGtttta                                                                                                                                            >  2:1057156/1‑146 (MQ=255)
       aaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGTTTCAgg                                                                                                                                          >  1:1039390/1‑150 (MQ=255)
            cATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGTTTTAGgccgcc                                                                                                                                     >  2:682057/1‑150 (MQ=255)
             aTCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTATCGGGTTTAGggcgcc                                                                                                                                     >  1:509744/1‑149 (MQ=255)
             aTCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGgccgcc                                                                                                                                     >  2:1019176/1‑149 (MQ=255)
              tCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGTTTCTGgccccc                                                                                                                                     >  2:194491/1‑148 (MQ=255)
                 ttCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGa                                                                                                                                   <  2:1059794/147‑1 (MQ=255)
                  tCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATTGGTTTCAGGCCGCCGCGCGcc                                                                                                                              >  1:517448/1‑151 (MQ=255)
                   cAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGGTTCAGGCCCCCGGGCGcc                                                                                                                              >  2:569654/1‑150 (MQ=255)
                   cAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGGTTCAGGGCGCCGAgcgc                                                                                                                               >  1:377476/1‑149 (MQ=255)
                         ctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCCCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTATCTGGTTTTACCccgccgccgcgcgcag                                                                                                                          >  1:493889/1‑139 (MQ=255)
                            gATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACa                                                                                                                     <  1:274680/150‑1 (MQ=255)
                             atatTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAATGGGTTCCGGCCGCCGAGCGGCAGCCGAACaa                                                                                                                    >  2:801268/1‑150 (MQ=255)
                               atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTTAGGGCGCCCGGCCCCCGCAGaaacaaa                                                                                                                  >  1:646954/1‑145 (MQ=255)
                                 tttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGTTTCATGCCGCCCAGCGCCAGCAGAACaa                                                                                                                    >  2:141509/1‑146 (MQ=255)
                                 tttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATAGGGTTTAAGCCGCCGCGCGCCAGCAGAACAAt                                                                                                                   >  2:123254/1‑147 (MQ=255)
                                   tGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGa                                                                                                              <  2:935349/150‑1 (MQ=255)
                                   tGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGGTTCAGGCCGCCGGGCGCCAGCAGAACAATACCGa                                                                                                              >  2:680810/1‑150 (MQ=255)
                                    gggggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTTAGGCCCCCGAGCGCCAGCCGAACAATAACGAt                                                                                                             >  2:777186/1‑150 (MQ=255)
                                       ggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGGTTTAGGCCGCCGAGCGCCAGCAGAACAATACCGATaa                                                                                                           >  2:433470/1‑149 (MQ=255)
                                        gTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGc                                                                                                        <  2:731127/151‑1 (MQ=255)
                                            attaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAACGGGTTTTTTCCGCCGACGGCCTGCAGGAAAATAACGAAAATGCCa                                                                                                      >  2:388332/1‑149 (MQ=255)
                                             ttaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGGTTTAATAGGGTTTAGGCCGCCGAGCGCCAGCGGAACAAAACCCGTAATGCCaa                                                                                                     >  1:179560/1‑149 (MQ=255)
                                                 tCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTGATGGGGTTCTGGGCGCCGAGGGCCCGCAGAACAAAACCCATAAATCCAAccca                                                                                                 >  1:107568/1‑147 (MQ=255)
                                                   tttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAACGGTTTTAGGGCGCCGCGCGCCACCAGAACAAAAACGATAAAGCCAAACAGaa                                                                                               >  1:719816/1‑149 (MQ=255)
                                                    ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGa                                                                                                <  2:340093/147‑1 (MQ=255)
                                                     tCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTg                                                                                            <  2:398529/150‑1 (MQ=255)
                                                      ccAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGAc                                                                                          <  2:199037/151‑1 (MQ=255)
                                                        aGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAACGGTTTCAAGCCCCCGAGCGCCAGCAGAACAATAACGATAAAGCCAACCAGAATTGAcc                                                                                         >  2:109026/1‑150 (MQ=255)
                                                        aGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAACGGGTTTAGGCCGCCGAGAGCCAGCCCAACAATACCCATTATGCCCACCACAAATGagc                                                                                         >  1:203126/1‑148 (MQ=255)
                                                        aGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGCTTGAATCCGTTTTAGGCCCCCGTGGGCCAGACGAAAAAAACCCAACCTGACAACCCGAATTGatc                                                                                         >  1:776752/1‑148 (MQ=255)
                                                           tCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAAc                                                                                      <  2:328518/150‑1 (MQ=255)
                                                              gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAAcc                                                                                     <  2:883998/148‑1 (MQ=255)
                                                               tttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGGTTTAGGCCGCCGAGCGCCCGCAGAACAATAACGATAATGCCAACCAGAAATGACCAAcc                                                                                     >  2:680690/1‑147 (MQ=255)
                                                                ttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGCATCGGTTTCAGGGCGGCGCAGGCCGGCAGAAAAATACCGAAAAATTCAACCAAAATTGAACAAcc                                                                                     >  1:84395/1‑146 (MQ=255)
                                                                 tCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATAc                                                                                <  2:863040/150‑1 (MQ=255)
                                                                  cgcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGGGCGCCCGCAGAACAATAACGATAATGCCCACCAGAAATGACCAACCGAt                                                                                  >  2:1045192/1‑147 (MQ=255)
                                                                       cTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATTGGGTTTAGGCCCCCGAGGGCCAGCAGTACAAAAACGATAATAGCCACCAGAAAAGACCCACCCATACGc                                                                              >  2:887498/1‑146 (MQ=255)
                                                                       cTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGGTTTAGGCCGCCGAGCGCCCGCAGAAAAATACCGATAATGCCCACCAGAAATTACCCACCGAAACGCaac                                                                           >  2:282236/1‑147 (MQ=255)
                                                                        tttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGGTTCAGGCCGCCGAGCGCCAGCAGAACAAAACCGATAAAGCCAACCAGAATTGACCAACCCAAACGCAcc                                                                           >  2:327622/1‑148 (MQ=255)
                                                                              aaaGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTATCGGGTTTAGGCCGCCGCACGCCCGACGAAAAAAAACGATAATGCCCACCAGAAATGAACAACCAAAACGCaac                                                                           >  1:369743/1‑140 (MQ=255)
                                                                              aaaGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGTTTTAGGCCGCCGAGCGCCCGCAGAAAAAAAACGATAATGCCAACCAAAATTTACCAACCGAAACGCACCAGCAGAgg                                                                   >  2:926242/1‑150 (MQ=255)
                                                                                aGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGTATTGGTTTCTCGCCGCCCAGCGCCCGCAGAACAAAAACCCTAATGCCAACCCGAAATTACCCAACGAAAAGCACCAGCAGAggggg                                                                >  1:969629/1‑151 (MQ=255)
                                                                                   agagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGGTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAggggggt                                                              >  2:908288/1‑148 (MQ=255)
                                                                                       cGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATAGGTTTCAAGCCGCCGGGCCCCCACCGAAAAATAACCATAATTCCCACCACAAATGACCCACCGACAACCCCCAGCAGAGGTGGgtctt                                                           >  1:244209/1‑149 (MQ=255)
                                                                                         tCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGGTTCAGGCCGCCCAGCGCCCGCAGAAAAAAAACGATAATGCCAACCAGAAATGACCAACCGATACGCACCAGCAGAGGTGGGTCTTTAcc                                                       >  2:805572/1‑151 (MQ=255)
                                                                                            ccATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAATCGGTTTTAAGCCGCCGGGCGCCAGCAGAAAAAAAACCATAAAGCCAACCAGAATTGACCAACCCATAAGCCCCAGCAGAGGGGGGTCTTCACCaac                                                    >  1:532037/1‑149 (MQ=255)
                                                                                                 aTGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTGATGGGGTTCAGGGGCCCGCGCGGCAGCAGGACAATACCCGTAAAGCCAACCAAAAATAACCAACCGATAACCACCAGCAGAGGGGGGTCTTTACCAACGCGt                                                >  1:735663/1‑150 (MQ=255)
                                                                                                   gTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTAc                                              <  1:847538/150‑1 (MQ=255)
                                                                                                       aCGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAACGGGTTTAAGCCGCCGCGCGCCCGCCGAAAAAAACCGATAAAACCAACCAGAAATTACCCACCGATACCCACCAGCAGAGGGGGGTTTTCACCATCGCGTACtt                                            >  1:933867/1‑148 (MQ=255)
                                                                                                         gaagaaCAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCgg                                       <  2:75361/151‑1 (MQ=255)
                                                                                                           agaaCAGCGGGCATCCTCCGGCGATAATGGCGGGTTTAATCGGGTTCAGGCCGCCGAGCGCCAGAAGAAAAATACCGATAATTCCAACCACAATTGTCCAACCGAAACCCACCACCAGAGGGGGGTCTTCACCCACGCGTAATTcgcg                                        >  2:257366/1‑148 (MQ=255)
                                                                                                             aaCAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTAACGGTTTCAGGCCGCCGAGCGCCAGCAGAACAAAACCGATAATGCCAACCAAAATTGAACAACCCATACGCACCAACAGAGGGGGGTCTTCACCATCCCGTACTTCTCgg                                       >  1:715446/1‑147 (MQ=255)
                                                                                                              aCAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGt                                  <  2:387965/151‑1 (MQ=255)
                                                                                                                    ggCATCCTCCGGCGATAATGGCGGTTTTAATCGGGTTCAAGCCGCCGAGCGCCCCCAGAACAATACCGATAATGCCAACCAGAATTGACCCACCGATAAACACCCGCAGAGGGGGGTCTTTACCATCGCCTACTTCGCGGCAAATgg                                >  2:773249/1‑147 (MQ=255)
                                                                                                                    ggCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACAt                            <  1:884875/151‑1 (MQ=255)
                                                                                                                     gCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATc                           <  2:254673/151‑1 (MQ=255)
                                                                                                                      cATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCg                          <  1:130629/151‑1 (MQ=255)
                                                                                                                         cctccGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGcc                        <  1:1019027/150‑1 (MQ=255)
                                                                                                                                 gATAATGGCGGTTTGTAACGGTTTCAGGCCGCCGCGCGCCAGCAGAAAAAAAACGGTAATGCCCACCAGAATTTACCAACCGATAAGCACCCACAGAGGTGGGTCTTTACCCACGCGTACTTCGCGGCAAGTGGACATCGCCCggg                    >  1:132802/1‑146 (MQ=255)
                                                                                                                                     aTGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGGGGACATCGCCAGGGGATAagag            >  2:249912/1‑150 (MQ=255)
                                                                                                                                        gCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAgg        >  1:6016/1‑151 (MQ=255)
                                                                                                                                         cGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGc       <  2:444725/151‑1 (MQ=255)
                                                                                                                                            tttGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGtt    >  2:756796/1‑151 (MQ=255)
                                                                                                                                             ttGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGtt    <  2:403172/150‑1 (MQ=255)
                                                                                                                                             ttGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGtt    <  2:399615/150‑1 (MQ=255)
                                                                                                                                               gAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTaa  <  1:680697/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
TAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAA  >  NC_000913/40372‑40664

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