Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6050TC41.0% 112.1 / 30.7 78G187G (GGA→GGGycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (36/10);  new base C (0/32);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.90e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGAT  >  NC_000913/1174457‑1174754
                                                                                                                                                     |                                                                                                                                                     
acaacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTcc                                                                                                                                                     <  2:328921/151‑1 (MQ=255)
acaacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTcc                                                                                                                                                     >  1:546807/1‑151 (MQ=255)
acaacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTcc                                                                                                                                                     >  2:1194540/1‑151 (MQ=255)
    ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccac                                                                                                                                                 <  2:1260595/151‑1 (MQ=255)
     cAGGCACCGACCGCCTGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                                                                                                                                <  2:592885/151‑1 (MQ=255)
      aGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                                                                                                                                >  2:455776/1‑150 (MQ=255)
        gCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACg                                                                                                                                             <  1:1407699/151‑1 (MQ=255)
          accgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGc                                                                                                                                           >  2:1404440/1‑151 (MQ=255)
            cgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGt                                                                                                                                         <  1:221608/151‑1 (MQ=255)
                cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTggg                                                                                                                                      >  1:1032122/1‑150 (MQ=255)
                  ccAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAAc                                                                                                                                   >  2:457809/1‑151 (MQ=255)
                  ccAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAAc                                                                                                                                   >  1:1422603/1‑151 (MQ=255)
                   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:230487/1‑151 (MQ=255)
                                                                                                                                       aTCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTgg               >  1:953381/1‑151 (MQ=255)
                                                                                                                                              ttttCATCCCCCCCCCGGGGGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg          <  1:965721/147‑1 (MQ=255)
                                                                                                                                              tGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAagctggagctgg         >  2:1267964/1‑150 (MQ=255)
                                                                                                                                                     tCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGat  >  1:1060620/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGAT  >  NC_000913/1174457‑1174754

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: