Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913730,7440AC49.1% 70.9 / 12.9 59T388P (ACC→CCC) rhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (23/7);  new base C (0/29);  total (23/36)
Fisher's exact test for biased strand distribution p-value = 1.74e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGC  >  NC_000913/730597‑730892
                                                                                                                                                   |                                                                                                                                                    
gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacc                                                                                                                                                    >  1:969065/1‑150 (MQ=38)
 atgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccat                                                                                                                                                  >  1:214201/1‑151 (MQ=38)
  tgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatc                                                                                                                                                 >  1:848268/1‑151 (MQ=38)
    atAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatcac                                                                                                                                               >  2:419822/1‑151 (MQ=38)
        aTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGAc                                                                                                                                           >  1:418988/1‑151 (MQ=38)
         tACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGAc                                                                                                                                           >  1:586872/1‑150 (MQ=38)
         tACCGGGGCCGGATGGCGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGAc                                                                                                                                           >  1:861962/1‑150 (MQ=38)
                   ggATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccg                                                                                                                                >  2:55739/1‑151 (MQ=40)
                        gtggCCCCCCGTCACAGGGGCCGACCGGAAATTTTTTACCTTTACGACGGGGAGGGGGGGGTAACAGAACAGCTAACCCCGGCAGTCTAAAGTTACACGTACCAGTATGAGAAGGACCGCATCCCCACCCCCGCCAGCCTGAACCGCCGTg                                                                                                                           <  2:644950/151‑1 (MQ=255)
                          ggccccccGTCACCCGGGCCGCCCGGAAATCTGTTCCCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAACCCCGGCAGGCTTAAGTTACACGTATAAGTAGGAGAAGGCCCCCACCCCCATCCCCGACCGCCTGAACCGCCGTGaa                                                                                                                         <  2:998500/145‑1 (MQ=14)
                                         gggCCGACCGGAAATCTGTTCCCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAACCCCGGCAGGCTTAAGTTACACGTATCAGTAGGAAAAGGACCCCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCacac                                                                                                              <  1:397176/147‑1 (MQ=14)
                                         gggCCGACCGAAAATCTGTTACCGTTACGACAGCGACGGGGGGGTGAAAGAACAGCAAACCCCGGCAGGTTTAAGTTACACGTATCAGTATGAGAGGGCCCGCACCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGca                                                                                                                 <  1:1014301/144‑1 (MQ=14)
                                            ccGCCCGGAAATCTGTTCCCGTTACGACAGCGACGGGCGGGTGACAAAACAGCTAACCCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCACCACCATCCCCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGt                                                                                                       <  1:786431/151‑1 (MQ=14)
                                                   gAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCg                                                                                                >  2:823935/1‑151 (MQ=255)
                                                    aaaTCTTTTCCCGTTACGACAGCGGGGGGGGGGGGACAGAACAGCAAACCCCGGCAGGCTTAAGTCACACGTATCAGTAGGAGAAGGACCGCATCACCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCg                                                                                                <  2:135551/150‑1 (MQ=255)
                                                    aaaTCTGTTCCCGTTACGACAGGGAGGGGCGGGTGAAAGAACAGCTAACCCCGGCGGGCTTAAGTTACACGTATCAGTATGAAAAGGACCGCATCCCCATCACCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCg                                                                                                <  2:1008027/150‑1 (MQ=255)
                                                    aaaTCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGCGAAGGCgg                                                                                               >  1:104843/1‑151 (MQ=255)
                                                      aTCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGc                                                                                             >  2:1007436/1‑151 (MQ=255)
                                                       tCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCt                                                                                            >  2:20232/1‑151 (MQ=255)
                                                           ttACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGa                                                                                          >  2:986880/1‑149 (MQ=255)
                                                            tACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAg                                                                                        >  2:990831/1‑150 (MQ=255)
                                                                 ttgcGACAGCGACGGGGGGGGGACAGAAAAGCTAACCCCGGCAGGTTAAAGTTACACGTATCAGTAGGAGAGGGCCCGCATCCCCACCCCCGACGGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCggg                                                                                    <  2:777020/146‑1 (MQ=255)
                                                                    cGACAGCGACGGGGGGGTGACAGAACAGCAAACCCCGGCAGGCTTAATTTACACGTATCAGTAGGAGAAGGCCCGCATCCCCATCACCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGggtgg                                                                                 <  1:549349/149‑1 (MQ=255)
                                                                     gACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGgtggtg                                                                               >  1:99436/1‑150 (MQ=255)
                                                                      aCAGCGACGGGCGGGTGAAAAAACAGCAAAACCCGGCAGGTTTAAGTTACACTTTTCAGTATGAGAGGGCCGGAATCACCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGa                                                                              <  1:74034/150‑1 (MQ=14)
                                                                       cAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGaa                                                                             >  2:194155/1‑150 (MQ=255)
                                                                         gcggcgGGCGGGTAACAGAACAGCTAACCCCGGCAGGCTTAAGTTCCCCGTACCAGTAGGAGAAGGACCGCATCCCCACCACCGACGGCTGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGaa                                                                             <  2:344070/148‑1 (MQ=14)
                                                                          cggcgggcgggTGACAGAACAGCAAAGCCCGGCAGGCTAAAGTTACCGTTACAAGTATGAGAAGCACCGCATCCCAAGCCCCGCCGGCTGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGaaaaa                                                                          <  2:848150/147‑1 (MQ=14)
                                                                              gggcgggTGACAGAACAGCAAACCCCGGCAGGTTTAAGTTACACGTATCAGTATGAGAGGGCCCGCATCCCCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAgg                                                                        <  1:589740/148‑1 (MQ=14)
                                                                              gggcgggTGAAAAACCAGCAAAACCCGGCGGGCTTAAGTTACACGTACCAGTAGAAGAAGAACCGCACCACCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAg                                                                         <  2:645829/147‑1 (MQ=14)
                                                                               ggggggTGACAGAACACCTAAACCGGGCAGCCTTAATTTACACTTCTCAGTAGGAGAAGGACCGCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAaca                                                                    <  1:579188/148‑1 (MQ=14)
                                                                               ggggggTGACAAACCTGCTAACCCCGGAAGGTTTAAGTTACCCGTATCAGTTTGAGAAGGCCCCCTTCCCCAACGCCGACGGCTTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAaca                                                                    <  2:855387/148‑1 (MQ=14)
                                                                               ggggggTGACAAAACAGCTAACCCGGGCAGGCTTAATTTACACGTATCAGTAGGAGAAGGCCCGCATCCCCATCCCCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAaca                                                                    <  2:1030896/148‑1 (MQ=14)
                                                                               ggggggTGACAAAACAGCTAACCCCGGCGGGCTAAATTTACACGTACCAGAATGAGAAGCCCCCCACCCCCATCACCGACAGCCAGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAaca                                                                    <  2:422714/148‑1 (MQ=14)
                                                                               ggcgggTAAAAGAACAGCAAACCCCGGCGGGCTGAAGTTACACGTATCAGTAGGAAAAGACCCCCACCCCCCTCACCGACAGCCGGATCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAaca                                                                    <  2:207026/151‑1 (MQ=14)
                                                                                  gggTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcgc                                                                 >  1:799552/1‑151 (MQ=255)
                                                                                  gggTGACAAAACAGCTAACCCCGGCAGGCTAAGTTTACACGTATCAGTATGAGAAGGCCCGCACCCCCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcg                                                                  <  2:655256/150‑1 (MQ=14)
                                                                                        caaaaCAGCAAACCCCGGCAGCCTTAAGTTACCCGTATCAGTGTAGGAAGTCCCACACCCCCCGCCCCGACGGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGAc                                                             <  2:774007/146‑1 (MQ=14)
                                                                                           aaCAGCAAACCCCGGCGGGCTTAAGTTACACGTACCAGTAGGAGAAGGCCCGCATCCCCACCCCCGACAGCCTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAg                                                        <  2:848243/151‑1 (MQ=14)
                                                                                            agaGCAAAACCCGGCGGGCATAAGTTACACGTAGCAGTATGAGAAGGCCCGCATCCCCAGCTCCAACGGCAAGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGc                                                       <  2:842726/149‑1 (MQ=14)
                                                                                                cTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTc                                                    >  2:804547/1‑150 (MQ=255)
                                                                                                cTAAACCCAGCAGGCTTAAGTTCCACGTTCTAGTAGGAGAAGGGCCACAGCCCCAACCCCGACTGCCGAATCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGc                                                       <  1:990862/147‑1 (MQ=14)
                                                                                                  aaccccGGCGGGCTTAAGTTACACTTATCGGTATGAGAAGGCCCGCATCACCATCACCGACAGCGTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCAc                                                  <  2:375277/147‑1 (MQ=14)
                                                                                                  aaaCCGGGCGGGCTTAAGTTACCCGTATCAGTATGAGAAGGACCGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGc                                                       <  1:902502/145‑1 (MQ=14)
                                                                                                           aGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTc                                          >  1:200085/1‑149 (MQ=255)
                                                                                                             gCTTAAGTTCCCGGACTCAGTAGGAGAAGGACCGCACCACCACCCCCGACAGCTTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGtt                                      <  2:394665/151‑1 (MQ=14)
                                                                                                               ttAATTTACACGTATCAGTATGAGAAGGACCGCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGttt                                     <  1:480594/150‑1 (MQ=255)
                                                                                                                tAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTg                                    >  2:42764/1‑150 (MQ=255)
                                                                                                                            aTCAGTATGAGAAGGCCCGCACCCCCACCCCCGACGGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCa                        <  2:850065/150‑1 (MQ=14)
                                                                                                                              cAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGc                     >  2:17381/1‑151 (MQ=255)
                                                                                                                               attATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCt                    >  1:95473/3‑151 (MQ=255)
                                                                                                                                   tGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAgg                >  1:439360/1‑151 (MQ=255)
                                                                                                                                      gaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAggg               >  1:253793/1‑149 (MQ=255)
                                                                                                                                      gaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGc              >  2:422220/1‑150 (MQ=255)
                                                                                                                                       aaGGACCGCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAgg                <  1:832162/147‑1 (MQ=255)
                                                                                                                                        atacccacatccccatcgcTAACACCCAGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACGCGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAggg               <  1:925641/140‑1 (MQ=14)
                                                                                                                                         ggcccGCACCCCCATCCCCGACAGCCGAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGcac            <  2:586853/146‑1 (MQ=14)
                                                                                                                                         ggcccGCACCCCCATCACCGACAGCGTGAGCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAgg                      <  2:3459/136‑1 (MQ=14)
                                                                                                                                         ggACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGcac            >  2:1045401/1‑149 (MQ=255)
                                                                                                                                            ccGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGa         <  1:90654/149‑1 (MQ=255)
                                                                                                                                            ccGCATCCCCACCACCGACGGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGAc        <  2:313133/150‑1 (MQ=255)
                                                                                                                                            ccGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGCCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACg       >  2:387411/1‑151 (MQ=255)
                                                                                                                                            ccGCATCACCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGAc        <  1:789429/150‑1 (MQ=255)
                                                                                                                                             cGCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACgg      <  1:441970/151‑1 (MQ=255)
                                                                                                                                                 tcaccatcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGACAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGc  >  2:355535/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGC  >  NC_000913/730597‑730892

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