Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA45.3% 76.1 / 49.5 75intergenic (‑330/+60)ybgD/gltAputative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (30/11);  new base A (0/34);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 5.47e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCA  >  NC_000913/752977‑753273
                                                                                                                                                       |                                                                                                                                                    
aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                                                     >  2:213164/1‑150 (MQ=255)
aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaa                                                                                                                                                    <  1:781794/151‑1 (MQ=255)
 ttttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaa                                                                                                                                                    >  2:758834/1‑150 (MQ=255)
     cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAAt                                                                                                                                                <  1:669735/150‑1 (MQ=255)
      aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                               >  2:313451/1‑150 (MQ=255)
      aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                               >  1:456756/1‑150 (MQ=255)
      aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCa                                                                                                                                              <  2:267467/151‑1 (MQ=255)
       tAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCa                                                                                                                                              >  2:774988/1‑150 (MQ=255)
         gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                                           >  1:327806/1‑151 (MQ=255)
                     cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAc                                                                                                                               >  2:831332/1‑151 (MQ=255)
                         tACCATGAAGTTTTTAATTGAAAGTATTGGGATGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcg                                                                                                                           >  2:1035623/1‑151 (MQ=255)
                               gAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTa                                                                                                                      >  2:1025848/1‑150 (MQ=255)
                                aattttttAATTGAAAGTTTTGGGTTGCTGATAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTATTTTTACAATTTTTTGTTTTAAGGAATATAAAAATCACCCCGCCATATGAACGGCGGGTTaaa                                                                                                                    <  2:771079/148‑1 (MQ=255)
                                aaGTTTTAAATTGAAAGAATGGGTTGGCTAATTGTTTGGGTTGTTCTTTTTTTTTAAAATATCTAAATAGGTCTTTAAATGGATTTTATTTTGAAAATTTTTGGTGTTTAGGCAAAAAAAAAACAACCCGCCATATGAACGGCGGGTTaag                                                                                                                    <  1:402045/151‑2 (MQ=255)
                                   tttttAATTAAAAGTATTGGGTTGCTGATATTTTGAGCTGTTCTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTaaaa                                                                                                                   <  2:84850/149‑1 (MQ=255)
                                    ttttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTaaaa                                                                                                                   >  1:1033388/1‑148 (MQ=255)
                                    tttaaatttaaatGTATTGGTTTAGTTATAGTTTGCGTTGTTTTATTCTTTTAAAAAATTTATATGGGTATTTAAAGGTTTTTTTTTTTACCATTTTTTTTTTTTTGGGCTTATAAAAAACACCCCGCCATATGAACGGCGGGTTAAAat                                                                                                                  <  2:357760/150‑1 (MQ=255)
                                       tAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                              >  1:131730/1‑150 (MQ=255)
                                                   aTTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTAAAATTAga                                                                                                    >  1:186001/1‑148 (MQ=255)
                                                     tggtttcttGATAATTTGGGTTTTTTTTTTTTTTTAAAATATTTTTATAGTTCTGTTAATGGTTTTTTTTTTTACAATTTTTTGTTTTTGGGCATATAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                                 <  2:427646/142‑1 (MQ=255)
                                                     tGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                                 >  1:118117/1‑150 (MQ=255)
                                                      gggTTGTTGATATTTTGAGCTTTTCTTTTTTTTTAAAATATCTATATAGGTTTTTAAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                                 <  2:971857/149‑1 (MQ=255)
                                                      gggTTGCTAATAATTGGGGTTTTTTTATTTTTTTTAAATATAAAAATAGGTCTTTTATTGATTTTTTTTTTTACATTTTTTTGTTTTTAGACAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATc                                                                                                <  1:403730/150‑1 (MQ=255)
                                                         ttGCTGATAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTTTTTTAATGGATTTTTTTTTTACATTTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                                              <  2:739375/149‑1 (MQ=255)
                                                               aTAATTTGAGCTTTTTTATTTTTTTAAAATACTTATATGGTTTTTTAAATGATTTTTATTTTAATATTTTTTTTTGTTAAGGCAAATAAAAATAACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                                        <  2:930654/149‑1 (MQ=255)
                                                                tAATTTGAGCTTTTCTTTTTTTTTAAAATATTTTTATGGGTCTGTTAAGGTATTTTATTTTTAAAATTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTaa                                                                                       <  2:301790/149‑1 (MQ=255)
                                                                 attttGAGCTTTTTTATTTTTTTAAAATATCTAAATAGGTCTGTTAATGAATTTTTTTTTTACATTTTTTTGTTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTaa                                                                                       <  2:649721/146‑1 (MQ=255)
                                                                 aaTTTGAGTTGTTCTTTTTTTTTAAAATATCTATATAGTTCTGTTAATGTATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                                     <  2:198460/150‑1 (MQ=255)
                                                                   tttGTGCTTTTTTTTTTTTTTTAAATATTTAAATGGTTTTGTTAATGAATTTTATTTTTAAATTTTTTTGTTTTAAGAAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                                          <  1:28513/143‑1 (MQ=255)
                                                                   tttGAGTTTTTTTTTTTTTTTAAAATTTTTATATAGTTTTTTTAATGGATTTTATTTTTACATTTTTTTGTTTATAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAc                                                                                      <  1:930371/147‑1 (MQ=255)
                                                                    ttggtttttttttttCTTTTAAAATTTTTATATATTTCTTTTAATGTTTTTTTTTTTTACAATTTTTTGTTTGTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                                     <  1:680405/134‑1 (MQ=255)
                                                                      gagtttttCTATTTTTTTTAAATATTTATATAGTTTTTTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGAATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCtt                                                                                  <  2:906777/142‑1 (MQ=255)
                                                                          ttttttttttttttaAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTAAAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGata                                                                              <  2:862263/138‑1 (MQ=255)
                                                                          tGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTTGCAAACAACAATTAACACTTGATCTCg                                                                           >  1:280681/1‑151 (MQ=255)
                                                                            tttttttttttttAAATATCTATATAGTTTTTTTAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCg                                                                           <  1:911044/141‑1 (MQ=255)
                                                                            ttttttttttttaAAATATTTATATGGGTTTTTTAATGGATTTTTTTTTTACAATTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                          <  2:665742/142‑1 (MQ=255)
                                                                            ttCTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCg                                                                           <  2:43424/149‑1 (MQ=255)
                                                                             tCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtt                                                                        >  2:1025076/1‑151 (MQ=255)
                                                                                 ttCTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAGTCAACCCGCCATATGAACGGCGGGGTTAAATTTTTTCAACTTAAAAATAAAACAATAACGCTTGAAATCTCttt                                                                       >  1:605357/1‑148 (MQ=255)
                                                                                  tttttttAAATATTTATATAGTTTTTTAAATGTTTTTTTTTTTTACATTTTTTTTTTTTTAGAAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaa                                                                    <  1:831288/148‑1 (MQ=255)
                                                                                  tCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaa                                                                    >  2:1306/1‑150 (MQ=255)
                                                                                    tttttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAACAATCAACCATTAACGCTTGATATCTCTTTTaaa                                                                   >  1:66696/1‑149 (MQ=255)
                                                                                       taaaaTATTTATATAGGTCTTTTAATGTATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTc                                                                <  1:807044/147‑1 (MQ=255)
                                                                                          attatCTATATAGGTTTGTTAATGTATTTTTTTTTTACAATTTTTTGTTTTTAGATATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcg                                                               <  1:831546/145‑1 (MQ=255)
                                                                                            tatCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                         >  2:83319/1‑151 (MQ=255)
                                                                                                 atataGGTTTTTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                         <  1:427330/146‑1 (MQ=255)
                                                                                                   ataGGTTTGTTAATGTTTTTTATTTTTAAATTTTTTTGTTTTTAGACTAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                                    <  1:112593/149‑1 (MQ=255)
                                                                                                    taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTATAAATTAACAATCAACCATTAAAGCCTGATATCGCTTTTAAAGTCGCGtttttta                                                      >  1:118125/1‑144 (MQ=255)
                                                                                                     aGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                                 >  1:925556/1‑150 (MQ=255)
                                                                                                      ggTTTGTTAATGGATTTTATTTTTAAATTTTTTTTTGTTTAGGCAAATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                                        <  2:27957/142‑1 (MQ=255)
                                                                                                       gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAACCATTAACGCTTGGTATCGCTTTTTAAGTCGCGTTTTTCATATCCTGt                                              >  1:908381/1‑151 (MQ=255)
                                                                                                       gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTCAAATTTTTTCAACTTAGCCATCAACCATTAACGCTTGGTATCGCTTTTAAAGTCGCGTTTTTCATATCCt                                                >  2:873310/1‑149 (MQ=255)
                                                                                                        tttGTTAATGGTTTTTTTTTTTACAATTTTTTGTGTTTAGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                              <  2:626371/148‑1 (MQ=255)
                                                                                                          ttttAATGGTTTTTATTTTTACAATTTTTTTTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtat                                            <  2:78944/148‑1 (MQ=255)
                                                                                                          tGTTAAGGGTTTTTTTTTTTAAATTTTTTTGTGTTTAGGAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                             <  1:1073724/149‑1 (MQ=255)
                                                                                                              aaTGGTTTTTATTTTTAAAATTTTTTGTGTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                        <  1:774945/150‑1 (MQ=255)
                                                                                                              aaTGGATTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                        <  2:999155/150‑1 (MQ=255)
                                                                                                              aaTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGt                                              >  1:861523/1‑144 (MQ=255)
                                                                                                               aTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAACCCTTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGc                                       >  1:404023/1‑150 (MQ=255)
                                                                                                                 ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCt                                      >  2:1010333/1‑149 (MQ=255)
                                                                                                                    ttttttttttAAATTTTTTTGTTTTTAGGTATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg                                  <  2:728427/150‑1 (MQ=255)
                                                                                                                     ttttttttttTAATTTTTTGTTTTTTAGAAAAAAAAAAACAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGa                                    <  1:742773/148‑1 (MQ=255)
                                                                                                                      ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                               >  2:521783/1‑151 (MQ=255)
                                                                                                                          tttaACATTTTTTTGTTTTTGGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACg                            <  2:517147/150‑1 (MQ=255)
                                                                                                                               aaaaTTTTTTTTGTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa                        <  1:518194/148‑1 (MQ=255)
                                                                                                                                 atttttttttttttAGTAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCaa                       <  2:544624/147‑1 (MQ=255)
                                                                                                                                   ttttttttGTTAAGGAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATc                     <  1:25563/149‑1 (MQ=255)
                                                                                                                                     ttttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt                   >  1:189343/1‑149 (MQ=255)
                                                                                                                                      tttGTGTTTAGTAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa                 <  1:520171/150‑1 (MQ=255)
                                                                                                                                       ttGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCAt                <  2:458414/150‑1 (MQ=255)
                                                                                                                                           tttttAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa               <  1:368246/146‑1 (MQ=255)
                                                                                                                                           tgtTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc             >  1:47265/1‑150 (MQ=255)
                                                                                                                                             tttAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt           <  1:758790/150‑1 (MQ=255)
                                                                                                                                             tttAGGAAAAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCg            <  1:313905/149‑1 (MQ=255)
                                                                                                                                                  gCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGGTTTTCATATCCTGTATACAGCTGa                                    >  2:344056/1‑120 (MQ=255)
                                                                                                                                                   cATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACtgt     >  2:520810/1‑150 (MQ=255)
                                                                                                                                                      ataAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgc                                 <  1:344036/119‑1 (MQ=255)
                                                                                                                                                      ataAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCa  <  1:724996/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCA  >  NC_000913/752977‑753273

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