Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6260GA46.0% 67.8 / 41.2 74F180F (TTC→TTTycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (30/10);  new base A (1/33);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 6.40e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACC  >  NC_000913/1174476‑1174775
                                                                                                                                                       |                                                                                                                                                     
cAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                                                                                                                                                       >  2:593012/1‑151 (MQ=255)
 aGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGa                                                                                                                                                      <  2:279033/151‑1 (MQ=255)
   aaCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaa                                                                                                                                                     >  1:974537/1‑150 (MQ=255)
      ggTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTa                                                                                                                                                 <  1:50906/151‑1 (MQ=255)
                  agCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                                                                     >  2:584230/1‑151 (MQ=255)
                     cATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                                                                     >  2:975449/1‑148 (MQ=255)
                            aatgggcggcGAACCGCCCCTTTTTTAACACAGGACTTTTTTGTTGCCCAAACCCAGGCGCAAGGGTTGAAAAAAGGAAAATGGCGGACCCCACTTTTAATCCCCCCCCCGGGCGGGGAAAAAAAATTAACCGCCACACTcagca                                                                                                                                 <  2:1365966/141‑1 (MQ=255)
                             gATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaa                                                                                                                           >  2:715384/1‑150 (MQ=255)
                               tGGGGGCGCAACGGCCCCTCTTTAACACAGGTAATTATTGTTGCGCAAAACCAGGCGCAAGGGCTGTAATAAGGAAAATGGCGAATCACCCTTTTCATCCCCCCCCACGGCGGGGGAACAAAATAAACCGCCACACTCAGCAGTaca                                                                                                                            <  2:1201074/147‑1 (MQ=255)
                                ggcggcGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                                        >  1:373627/1‑150 (MQ=255)
                                 gcggcgCAACGGCACCTTTTTAACACAGGTCTTTATTGTTGCGCCAAACCAGGCGCCAAGGCTGTAATAAGGCAAATGGCGGATCACACTGTTCATTCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                                        <  1:1065869/149‑1 (MQ=255)
                                   ggcgCAACGGCCCCTCTTAAACACAGGTCATTATTTTTGCGCAAACCCATGCCCCAAGGCTGTAATAAGGCAATTTGGGGATCACTCTTTTCATCCCCCCCCCCGGCGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAgg                                                                                                                     <  1:1434764/150‑1 (MQ=255)
                                     cgcAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGAAGTACAAACAAGGCg                                                                                                                   >  2:1280871/1‑150 (MQ=255)
                                     cgaaaCGCCCCCTCTTAAACACAGGTCATTTTTGTTCCCCAAACCCATGCCCCAGGGTTGAAATAAGGCAATTTGCGAATCACCCTGTTCTTCCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGt                                                                                                                               <  2:621473/135‑1 (MQ=255)
                                      gcAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGa                                                                                                                  >  1:1363916/1‑150 (MQ=255)
                                        aaCGCCCCCTCTTTAACACAGGTTATTATTGTTGCGCAAACCCATGCGCCAAGGCTGTAAAAAGGCAAATTGCGGATCCCACTGTTCATCCCCCCCCCCGGGGTGGGAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGaa                                                                                                                 <  1:282737/149‑1 (MQ=255)
                                          cgcccccTCTTTAACACAGGTCATTATTGTTGCGCAAACCCAGGCGCAAAGGCTGAAATAAGGCAAATTGCGGAACCCACTTTTCTTCCCCCCCCGGGGGGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                                                   <  2:102657/140‑1 (MQ=255)
                                              aCCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                                          >  1:177911/1‑150 (MQ=255)
                                                ctctTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgc                                                                                                       >  1:672343/1‑151 (MQ=255)
                                                ctctTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTAAAAACAAGGCGAATAGTGGCAg                                                                                                        >  2:990904/1‑150 (MQ=255)
                                                   tttAACAAAGTTCATTTTTGTTGCGCAAACCCATGCCCAAAGGCTGTAAAAAGGAAAATTGCGGACCACACTGTTCATCCCCCCCCGGGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgc                                                                                                       <  2:1298364/148‑1 (MQ=255)
                                                      aaCACAGGTAATTTTTGTTCCACCAACCCATCCGCCAGGGCTGAAAAAAGAAAATTTGCGGACCCCCCTTTTCTTCCCCCACCAGGGGGGGGAAACAAAATTACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAg                                                                                                  <  2:496996/150‑1 (MQ=255)
                                                         acaGGTATTTTTTGTTGCGCAAAACCATGCGCAAGGGCTAAAAAAAGGAAAATTGCGGATCACCCTTTTCATCCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                                               <  1:169605/150‑1 (MQ=255)
                                                            ggTAATTATTGTGGCCCAAACCCATGCGCAAAGGCTGTAAAAAGGAAAATTGCGGATCACACTTTTCTTCCCCCACCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                                            <  1:717159/150‑1 (MQ=255)
                                                             gTCATTATTGTTGCGCAAACCCGTGCGCAAAGGTTGTAAAAAGGAAATTTGCGGATCACACTGTTCTTCCCCCACCACGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACgg                                                                                              <  2:577413/147‑1 (MQ=255)
                                                              tCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                                          >  2:353611/1‑150 (MQ=255)
                                                               cttttttttggcgcCAAACATGGCGCCAAGGTTTTAATAAGCAAATTTGCGCACCTCACTGTTCATCCCCCCCCACGCCGGGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCa                                                                                                   <  1:127130/130‑1 (MQ=255)
                                                                 ttttggtggcgcCACACCCGGCCCCAAGGCTTTAAAAACCAATTTGCTGAACCCCATTTTTCCTCCCCCACCGCGGCGAGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                                                        <  2:262720/141‑1 (MQ=255)
                                                                 ttattTTCGCGCCAACCCTTGCGCCAAGCCGAAAATAAGGCAAATTGCGGACCCCACTTTTCTTCCCCCCCCGCGGGGGGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaaa                                                                                       <  2:669378/150‑1 (MQ=255)
                                                                  tattGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                                                     >  1:1078120/1‑151 (MQ=255)
                                                                  tattGTTGCGCCAAACCATGCGCAAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCTCACCACGGCGCGGGCACACAATTAACCGCCACACTCTGCAGTACAAACAAGTCGAATAGTGGCAGCGCCAGGCGGCTaaaaaat                                                                                     >  2:1392115/1‑149 (MQ=255)
                                                                  tattGTTGCGCCAAACCAGGCGCCAAGGCTGTAAAAAGGAAAATGGGGGATCACACTGTTCTTCCCCCCCCCCGGCGGGGGAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaaa                                                                                       <  2:1247811/149‑1 (MQ=255)
                                                                     tgttgCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTTTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtttt                                                                                   >  1:263243/1‑150 (MQ=255)
                                                                         gcgcAAAACAAGGCGCCAAGGTTGTAAAAAGGAAATTTGCGGACCCCACTTTTCATCCCCCACCCCGGGGTGGAAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                                               <  1:899215/150‑1 (MQ=255)
                                                                          cgcCAAACCATGCGCAAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCTTCCCCCACCAGGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                                               <  1:1019123/149‑1 (MQ=255)
                                                                           tcttccgatccccccAAGGCTGTAATAAGGCAAATGGCGGATCACCCTTTTCATCCCCCCCCCCGGGGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATa                                                                              <  1:1460942/136‑1 (MQ=255)
                                                                            ccAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                                            >  1:1383508/1‑150 (MQ=255)
                                                                            ccAAACCATGCGCAAGGGCTGTAAAAAGGAAAATTGCGGATCACACTTTTCTTCCCCCCCCCCGGCGTGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAc                                                                             <  2:698366/149‑1 (MQ=255)
                                                                              acccctgacgcAAAGGGTGAAAAAAGCAAAGTGGCGGACCACACTTTCCTTCCCCCACCGCGGTGTGGAAACAAAAATAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                                         <  2:70550/143‑1 (MQ=255)
                                                                              aaaCAATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACCCTTTTCTTCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAc                                                                          <  1:559586/150‑1 (MQ=255)
                                                                                aCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGCGGCAGCGACAGACGGCTAAAAATTTTACCTACCACg                                                                         >  2:297025/1‑149 (MQ=255)
                                                                                 ccaggccacAAGGCTGTAAAAAGCAATATTGAGCTTCACCATTTTCTTCCCCCACCGCGGCGGGGGAAAAAAATAATCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                                                      <  2:624516/143‑1 (MQ=255)
                                                                                 ccATGCGCCAAGGCTGTAAAAAGGAAAATTGCGGATCACACTGTTCATCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                                                      <  1:90703/151‑1 (MQ=255)
                                                                                    ggcgcCAGGGCTGTAAAAAGCAAAATTGCGGATCCCCCTTTTCATTCCCCCCCGGGGGGGGGAAACAAAATTACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATa                                                                    <  1:1180090/149‑1 (MQ=255)
                                                                                      cgcCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATa                                                                              >  2:1460863/1‑138 (MQ=255)
                                                                                        ccAGGGCTTTAATAAGGCAAATTGCGGATCACACTTTTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACa                                                                <  1:332187/150‑1 (MQ=255)
                                                                                        ccAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGCGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACa                                                                >  2:354541/1‑150 (MQ=255)
                                                                                          agggCTGAAAAAAGGCAAATTGCGGATCCCTCTTTTCTTCCCCCACCGGGGCGGGGAAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAtt                                                             <  2:389706/149‑1 (MQ=255)
                                                                                          aaGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAtt                                                             >  1:1470991/1‑151 (MQ=255)
                                                                                            ggttgaaaaaaGGAAAATGGCGGTCCCCCCTTTTCTTCCCCCCCCCCGGGGTGGAAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                                           <  2:194091/142‑1 (MQ=255)
                                                                                                gTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGc                                                        >  2:616136/1‑150 (MQ=255)
                                                                                                 taataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                                      >  2:1374727/1‑151 (MQ=255)
                                                                                                  aaaagggaaaaTGGCGAATCACCCTGTTCATCCCCCCCCCCGGCGTGGGAACGAAATGAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCa                                                       <  1:1003269/141‑1 (MQ=255)
                                                                                                   ataaGTCAAATTGGGAATCACTCTGTTCATCCCCCACCACGGCGGGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                    <  2:1350879/151‑1 (MQ=255)
                                                                                                      aGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                    >  1:1441292/1‑148 (MQ=255)
                                                                                                          aaaTTGCGGATCACACTGTTCATTCCCCACCACGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                                                <  1:128598/148‑1 (MQ=255)
                                                                                                           aattgggggATCACACTGTTCATCCCCCCCCCCGGGGGGGAAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                                             <  1:537668/144‑1 (MQ=255)
                                                                                                            aTTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACca                                            >  2:321784/1‑150 (MQ=255)
                                                                                                              tGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                                          >  2:425217/1‑150 (MQ=255)
                                                                                                               gCGGATCACCCTTTTCTTCCCCCCCCGGGGGGGGGAAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                                          <  2:858533/149‑1 (MQ=255)
                                                                                                                  ggATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGc                                       >  1:1464913/1‑150 (MQ=255)
                                                                                                                  ggATCACACTGTTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                                          <  1:226824/147‑1 (MQ=255)
                                                                                                                    aTCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTgg                                    >  1:953381/1‑151 (MQ=255)
                                                                                                                    aTCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTgg                                    >  1:230487/1‑151 (MQ=255)
                                                                                                                           ttttCATCCCCCCCCCGGGGGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg                               <  1:965721/147‑1 (MQ=255)
                                                                                                                           tGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAagctggagctgg                              >  2:1267964/1‑150 (MQ=255)
                                                                                                                             ttcttcccccccCACGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGc                             <  1:1181481/143‑1 (MQ=255)
                                                                                                                                  tCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGat                       >  1:1060620/1‑150 (MQ=255)
                                                                                                                                  cccccccCACGGGGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt                          <  2:1344345/146‑1 (MQ=255)
                                                                                                                                       cccaggcggggaaaaCAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATcc                  <  2:1177930/137‑1 (MQ=255)
                                                                                                                                         cgcggcggGAGAACGAATTAACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGt               <  2:517665/149‑1 (MQ=255)
                                                                                                                                         cacGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGCGATCCTccc                 >  1:360883/1‑149 (MQ=255)
                                                                                                                                          acGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATTCCCTg              >  1:1063206/1‑151 (MQ=255)
                                                                                                                                               gTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGa         <  2:318890/151‑1 (MQ=255)
                                                                                                                                                    aaCGAATTAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTc     <  2:598703/150‑1 (MQ=255)
                                                                                                                                                       gAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGaa        >  2:819652/1‑144 (MQ=255)
                                                                                                                                                       gAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCAcc  >  2:613632/1‑150 (MQ=255)
                                                                                                                                                       gAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCAcc  >  1:932271/1‑150 (MQ=255)
                                                                                                                                                       gAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGGGCAGGAACTCAcc  >  1:878576/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACC  >  NC_000913/1174476‑1174775

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