Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA42.2% 88.7 / 30.5 64intergenic (‑124/+79)ves/spycold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (31/6);  new base A (0/27);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.35e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAG  >  NC_000913/1824914‑1825207
                                                                                                                                                     |                                                                                                                                                  
tttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagta                                                                                                                                                  >  1:570191/1‑148 (MQ=255)
 ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                                                                 >  1:848360/1‑150 (MQ=255)
   aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                              >  1:11556/1‑151 (MQ=255)
   aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                              <  1:710536/151‑1 (MQ=255)
   aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                              >  1:298447/1‑151 (MQ=255)
    cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                             <  1:339121/151‑1 (MQ=255)
              aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGt                                                                                                                                   >  2:592260/1‑151 (MQ=255)
                  ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                               >  1:27184/1‑151 (MQ=255)
                   tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                              >  2:291531/1‑151 (MQ=255)
                       tcctcccccccGCCCCTTTTTTTTCATTTTTATCGTTCAGGTAACGGTGGTTGTTTAGCGGCACCTGGGGTGATTTACGTTAGGGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAcc                                                                                                                              <  2:604167/139‑1 (MQ=255)
                         cTCGCCCCCGCCCCGTTTTTTTAATTTTTATCGTTAGGGTAACCGTGATTGTTTAGCGGCACCTCGGGTGATTTCCTTTGTTGTGATTCAGGTTTTTTTTACAACTGTTGGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCa                                                                                                                           <  2:324322/148‑1 (MQ=255)
                                  ggCCCGTTCTTTTCATTGTTTTCGTTCAGGTACCGGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                <  1:909679/150‑1 (MQ=255)
                                      cGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                           >  1:807536/1‑151 (MQ=255)
                                          cTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                        >  2:568410/1‑150 (MQ=255)
                                             ttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                       >  2:762654/1‑148 (MQ=255)
                                             ttCATTGTTATCGTTAAGGTAACCGTGGTTGGTTATCGCAACCGGGGGTGATTTATTTTGGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAg                                                                                                     <  2:166200/150‑1 (MQ=255)
                                              tCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAc                                                                                                   >  2:757377/1‑151 (MQ=255)
                                                   gTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                              >  1:860021/1‑151 (MQ=255)
                                                     tATGGTTCAGTTAACCGTGGTTGGTTATGGGAACCTGGGGTGATTTACTTTGGTGGTAATCGGGTTTTTTTTAAACCTGTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                               <  2:430762/148‑1 (MQ=255)
                                                        cgttagggaaaCGGTGGATGGTTATCGGAACCTGGGGTGTTTTACGTTAGTGGTAATCAGGTTTTTTTTACACCTGTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                            <  1:716968/139‑1 (MQ=255)
                                                          ttCAGGTACCCGTGGATGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACATCTGTTGGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                            <  2:29219/146‑1 (MQ=255)
                                                          ttCAGGTACCCGTGGATGGTTATCGGCAACTCGGGTGATTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                                        <  2:127229/150‑1 (MQ=255)
                                                          ttCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCtt                                                                                         >  1:353935/1‑149 (MQ=255)
                                                             aGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt                                                                                    >  1:798347/1‑151 (MQ=255)
                                                              ggtaccgggggttGGTTACCGGCACCTGGGGTAATTTCCGTTGGTGGTGATCGGTTTTTTTTTACACTTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt                                                                                    <  2:571493/138‑1 (MQ=255)
                                                                 accggTGGATGGTTATCGGCAACTGGGGTTATTTACGTTGGTGGGGACCGGTTTTTTTTAACACTTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                                <  1:512444/147‑1 (MQ=255)
                                                                 aacgggggATGTTTTGCGGCAACTGGGGTGATTTGCGTTAGTGGTGATCGGTTTTTTTTTACATCTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                                <  1:396953/145‑1 (MQ=255)
                                                                 aaCCTTGGTTGGTTTCCGGCACCTGGGGTTATTTACTTTAGTGTTGACCAGGTTTTTTTTACAACTGTTTAAGAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGt                                                                                     <  1:22658/146‑1 (MQ=255)
                                                                     ggtgactttactCGGAAACGGGGGTAATTTCTTTTAGTGGTGATCTGGTTTTTTTAACATTTTTTTGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt                                                                                 <  2:835596/135‑1 (MQ=255)
                                                                      tGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGAAGTGGACAAGACCGGCGGACTTACGTTTTTTgg                                                                              >  1:895295/1‑148 (MQ=255)
                                                                           ggTTATCGGCAACTGGGGTGATTTACGTTGGGGGTGACCGGGCTTTTTTTACAACTTTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg                                                                       <  1:807495/150‑1 (MQ=255)
                                                                             ttATCGGAAACTGGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTTTTTGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGa                                                                          <  1:478870/145‑1 (MQ=255)
                                                                                 cGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                                 >  1:428779/1‑150 (MQ=255)
                                                                                  ggCAACTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                                <  1:77628/150‑1 (MQ=255)
                                                                                  ggCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                                >  1:494558/1‑150 (MQ=255)
                                                                                  gccAACTGGGTTGATTTAGGTTGGTGGTGATCGGGTTTTTTTTACAACTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTca                                                               <  2:450954/149‑1 (MQ=255)
                                                                                     acctcgtgttttttACTTTGGTGGTGATCAGGTTTTTTTTACATTTTTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                             <  2:392177/139‑1 (MQ=255)
                                                                                     aaCTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                              <  1:250321/149‑1 (MQ=255)
                                                                                     aaCTCGGGTGATTTACGTTAGGGGTGATCAGGTTTTTTTTACAACTTTTAGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                                                            <  1:738696/151‑1 (MQ=255)
                                                                                      aCTCGGGTGATTTACGTTAGGGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAAATTATTcagc                                                             >  2:923620/1‑149 (MQ=255)
                                                                                         cGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGtt                                                         >  2:696948/1‑150 (MQ=255)
                                                                                          ggttGATTTATGTTAGGGGTGATCAGTTTTTTTTTACAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGtt                                                         <  2:574521/146‑1 (MQ=255)
                                                                                            gTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                       >  1:132064/1‑149 (MQ=255)
                                                                                             tGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                                    >  2:181077/1‑151 (MQ=255)
                                                                                                 ttACTTTGGTGGTGACCAGGTTTTTTTTAAACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt                                                 <  1:141825/150‑1 (MQ=255)
                                                                                                  tACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAATTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt                                               >  2:791857/1‑151 (MQ=255)
                                                                                                       tAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt                                          >  1:521191/1‑151 (MQ=255)
                                                                                                           ggtgATCAGTCTTTTTTTACAACTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                      <  1:873893/151‑1 (MQ=255)
                                                                                                             ggtaacgagttttttttAACACTGTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                      <  2:221662/141‑1 (MQ=255)
                                                                                                              tgaccagtttttttttACACCTGTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                        <  2:721843/138‑1 (MQ=255)
                                                                                                               aTCAGTTTTTTTTAACATCTTTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc                                    <  2:272163/149‑1 (MQ=255)
                                                                                                               aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCg                                   >  1:628288/1‑150 (MQ=255)
                                                                                                               aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGc                                  >  2:35407/1‑151 (MQ=255)
                                                                                                                    gtttttttttACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGc                                  <  2:192774/144‑1 (MQ=255)
                                                                                                                    ggttttttttACAACTGTTTGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                               <  1:899862/147‑1 (MQ=255)
                                                                                                                          tttaACACCTGTTGGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGAc                             <  2:941251/146‑1 (MQ=255)
                                                                                                                            ttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                        >  1:268483/1‑149 (MQ=255)
                                                                                                                               cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                      >  2:49614/1‑148 (MQ=255)
                                                                                                                                aaCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAg                   <  1:762672/150‑1 (MQ=255)
                                                                                                                                  cTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg                >  1:530281/1‑151 (MQ=255)
                                                                                                                                   tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc               >  1:928985/1‑151 (MQ=255)
                                                                                                                                       ttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc            >  2:829203/1‑151 (MQ=255)
                                                                                                                                       tggagaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt             <  1:961289/148‑1 (MQ=255)
                                                                                                                                        tGAGAAGAGAAAGGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc            >  1:439987/1‑150 (MQ=255)
                                                                                                                                             agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa       >  1:755715/1‑150 (MQ=255)
                                                                                                                                             agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa       >  1:372757/1‑150 (MQ=255)
                                                                                                                                             aaagaaaagaaaCGCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCGGTCAGACGCTTCTCaaaa      <  2:312178/149‑1 (MQ=255)
                                                                                                                                               aaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt    <  1:49168/149‑1 (MQ=255)
                                                                                                                                                aaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa   <  2:361130/150‑1 (MQ=255)
                                                                                                                                                 aaaagaaagCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa      <  1:316373/147‑1 (MQ=255)
                                                                                                                                                 aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGCTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt    >  2:420987/1‑149 (MQ=255)
                                                                                                                                                 aaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg  <  1:610406/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAG  >  NC_000913/1824914‑1825207

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