Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,0890AC37.9% 189.9 / 10.8 105Y299D (TAT→GAT) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (52/12);  new base C (0/39);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 4.15e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

ACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACG  >  NC_000913/280941‑281239
                                                                                                                                                     |                                                                                                                                                       
aCCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATAcc                                                                                                                                                       >  2:1118167/1‑151 (MQ=255)
aCCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATAcc                                                                                                                                                       >  2:870452/1‑151 (MQ=255)
  ccGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCGCCCCTTCACCGTATTCCGGGGGCGTTGTGTTGCCGCTGTGCTGCCGCGCTGACGGCTGATACcgc                                                                                                                                                     <  2:268917/151‑1 (MQ=255)
  ccGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACcgc                                                                                                                                                     <  2:518060/151‑1 (MQ=255)
     ttCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACcgcg                                                                                                                                                    >  2:491000/1‑149 (MQ=255)
      tctcCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAg                                                                                                                                                  <  1:1368272/150‑1 (MQ=255)
      tctcCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGc                                                                                                                                                 >  1:1374813/1‑151 (MQ=255)
        tcCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGccc                                                                                                                                               >  2:674410/1‑151 (MQ=255)
         ccccTGAACGACTTGCCTGCGCTCAGACTTACCCCTGTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCCCCCCTTTATCGTATTCCGGGGGCGTTGTGTTTCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGccc                                                                                                                                               <  1:1325575/150‑1 (MQ=255)
          cccTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCgg                                                                                                                                             <  2:1307839/151‑1 (MQ=255)
          cccTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCgg                                                                                                                                             >  2:553703/1‑151 (MQ=255)
          cccTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCgg                                                                                                                                             >  2:1171571/1‑151 (MQ=255)
              gAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCCCCCCTGCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTa                                                                                                                                          <  1:287920/149‑1 (MQ=255)
              gAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTAc                                                                                                                                         <  1:828353/150‑1 (MQ=255)
                        ccTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATcc                                                                                                                               >  1:930619/1‑150 (MQ=255)
                        ccTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATcc                                                                                                                               >  2:605084/1‑150 (MQ=255)
                         cTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATcc                                                                                                                               >  2:805215/1‑149 (MQ=255)
                          tGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAgc                                                                                                                            >  1:835160/1‑151 (MQ=255)
                                cAGCCTTCCCCCTTTCACCCCCAGTTTCCCGCTAATATCCCCTTCCCTGCCCCCCCCCCCTTCCCAGTATCCCGGGGGGTTTTTGTTGCCCCTGTTCGGCCGCGCTGCCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTc                                                                                                                       <  1:549607/150‑1 (MQ=255)
                                 aGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCg                                                                                                                      >  1:811116/1‑150 (MQ=255)
                                   aCTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCgcgc                                                                                                                   >  1:316516/1‑151 (MQ=255)
                                     ttACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCg                                                                                                                  >  1:888893/1‑150 (MQ=255)
                                     tccccccTTTCCCGCTCACCTTTCCGCTGATATCCACTTTCCTCCCCACCCCCCCTTCATCGTATCCGGGGGGGGTTGTGTTGCCGTTGTACTGCCGCGCTGCCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGc                                                                                                                   <  2:808877/146‑1 (MQ=255)
                                        ccccTTCCCCGCTCACCTTCCCGCTGATATCCCCTTCCCTGCCCACCCCCCCTTCATCGTTTCCCGGGGGGTTTGTGTTCCCGCTGTCCTGCCGCGCTGACGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCg                                                                                                                  <  2:1374719/147‑1 (MQ=255)
                                          ccTTTCCCGCCCAGCTTCCCGCTAATATCCTCTTCCCGGACCACCCCCCCTTCATCGTATCCCGGGGGCTTTGTTTTCCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGtt                                                                                                            <  2:532130/151‑1 (MQ=255)
                                          ccTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGt                                                                                                             >  1:1307412/1‑150 (MQ=255)
                                            tttCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTc                                                                                                           >  2:838602/1‑150 (MQ=255)
                                             ttCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                                                                                          >  1:486140/1‑150 (MQ=255)
                                               cACGCGCAGCTTCCCGCTGATTTCCCCTTCCCGGACCACCCCCCCTACATGGTATCCGGGGGGCGTGGGTTTGCCGCTGTCCTGCCGCGCTGACGGCTGATGCCGCGAGCCTGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAg                                                                                                        <  2:930665/150‑1 (MQ=255)
                                                  gCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATCTCCAGCGCCTCGTGCGGGCGTTCAAGGtt                                                                                                     >  2:809834/1‑150 (MQ=255)
                                                   cccAGCTTCCCGCTGATCTCCCCTTCCCTGACCCCCCCCCCTTCATCGTATTCGGGGGGCTTTGTGTCCCCGCTGTACGGCCGCGCTGCCGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTtat                                                                                                   <  1:239108/149‑1 (MQ=255)
                                                   cTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTtat                                                                                                   >  2:108676/1‑151 (MQ=255)
                                                     ccgCTTCCCGCTTATATCCCCTTTCCTGCCCACCCCCCCTCCTGCGTCTTCCGGGGGCGTTGTGTTCCCTCTGTCCTGCCGCGCGGAGGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTtata                                                                                                  <  1:588552/148‑1 (MQ=255)
                                                      aGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTtat                                                                                                   >  1:1232385/1‑148 (MQ=255)
                                                      aGCTTGCCTCTAATATCCCCTTTCCTGCCCATCCCCCTTTCACCGTTTCCGGGGGGCGTGGGGTTTCCGCTGTACTGCCGCGCGGCCGGCTGATCCGGCGAGCCGGGTGCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGa                                                                                                <  2:420363/151‑1 (MQ=255)
                                                        cTTTCCGCTGATATCCACTTTCCTGCCCACCCCCCCTTCATCGTATTCCGGGGGCGTTGTGTTCCCGCTGTACTCCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGa                                                                                                <  2:69567/149‑1 (MQ=255)
                                                          ttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGt                                                                                            >  2:837958/1‑151 (MQ=255)
                                                            ccGCTGATACCCACTTTCCTGACCACCACCCCTTCATCGTATCCCGGGGGGGTTGTGTTGCCCCTGTCCTGCCGCGCTGCCGGCTGATACGGCGAGCCCGTTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGt                                                                                            <  2:1145046/149‑1 (MQ=255)
                                                            ccGCTAATCCCCACTTCCCTGACCATCACCCCTTCAGCGTATTCCGGGGGCGTTGTGTTGCCGCTGTCCTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGa                                                                                                <  2:59253/145‑1 (MQ=255)
                                                              gatGATACCCCCTTTCCTGACCACCCCCCCTTCATCTTATCCCGGGGGCGTGGTGGTGCCGCTGTCCTGCCGGGCGGGCGGCGGAACCCGGGGGCCGGGGGCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGa                                                                                                <  1:887098/141‑1 (MQ=255)
                                                                tGATATCCTTTTCCCCGCCCACCCCCCCTTCATCTTATTCGGGGGGCGTTGTGTTCCCGGTTCCCTGCCGCGCTGGCGGCTGATCCCGCGAGCCCGTTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGc                                                                                        <  1:428642/149‑1 (MQ=255)
                                                                tGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGcc                                                                                       >  2:549029/1‑150 (MQ=255)
                                                                tGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGcc                                                                                       >  1:937127/1‑150 (MQ=255)
                                                                 gATATCCCCTTTCCTGCCCACCCCCCCTTCATCTTATCCCGGGGGCTTTGTGTTCCCGTTGTACTGCCGCGCGGACGGCTGATCCCGCGGCCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCa                                                                                      <  2:1420346/150‑1 (MQ=255)
                                                                  atatCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGt                                                                                    >  1:822982/1‑151 (MQ=255)
                                                                   tatCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCTTCGTGCGGGCGGTCAAGGTTATAGACCGTCCGCCAGTg                                                                                   >  2:1270388/1‑151 (MQ=255)
                                                                    atCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTgg                                                                                  >  2:1132207/1‑151 (MQ=255)
                                                                      cccttttCCGCACCATCCCCCCTTCTTGTTATGCGGGGGGGTTTTTTTTGCCGCGGGCTGGCCGCGCTGGCGGCTGATACCGCGAGCCGGTCACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCa                                                                                      <  1:123659/141‑1 (MQ=255)
                                                                      ccccTTTCCTGACCTCCCCCCCTCCATCGTATTCGGGGGGCTTTGTGTTGCCGTTTTTCCCCCGCGATGACGGTTGATCCCGCGACCCGGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGt                                                                                 <  2:918129/147‑1 (MQ=255)
                                                                      ccccTTCCCTGCCCACCCCCCCTTCATCGTATCCGGGGGGCGTTGTGTTGCCGCTGTCCTGCCGCGCGGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTc                                                                                <  2:857411/148‑1 (MQ=255)
                                                                       ccccttccctacccATCCCCCCTTCATCTTTTCCGGGGGGCGTTGTGTTCCCGCTTTCCTCCCGCGCTGAGGGCTGATCCCGCGAGCCGGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTg                                                                                   <  2:1192804/136‑1 (MQ=255)
                                                                       ccccttCCTGCACCACCCCCCCTTCTGTTTATCCGGGGGGCGTTGTTGTGCCCCGGTCCTCCCCCGCTGCGGGCTGATCCCGCGCGCCCGGTGCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTgg                                                                                  <  1:539515/145‑1 (MQ=255)
                                                                         cTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGa                                                                              >  2:1018089/1‑150 (MQ=255)
                                                                           tcccTGCCCACCCCCCCTCCATCGTATCCCGGGGGCGTTGTGTTGCCGCTTTCCTCCCGCGCTGACGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAg                                                                            <  2:1080395/148‑1 (MQ=255)
                                                                              cTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgc                                                                        >  2:1268848/1‑151 (MQ=255)
                                                                               tGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgcg                                                                       >  2:1200875/1‑151 (MQ=255)
                                                                                   ccccccccccTTCATCTCATGCGGGGGGTGTTGTTTTGCGGGTGCACTGCCGCGCGGACGGGTGAACCCGCGGCCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgcgc                                                                      <  1:175117/143‑1 (MQ=255)
                                                                                   ccacccccccTTCATCGTATTCGGGGGGCTTTGTGTTGCCGCTTTCCTGCCGCGCTGACGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTg                                                                    <  2:1108617/145‑1 (MQ=255)
                                                                                   ccacccccccTCCTTCGTATTCGGGGGGCGTTGTGTTGCCGTTGTACTGCCGCGCGGACGGCTGATCCCGGGGGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCt                                                                     <  1:680744/144‑1 (MQ=255)
                                                                                   ccacccccccTCCATCGTTTCCGGGGGGCGTTGTGTTGCCGCTGTCCTGCCGCGCTGCCGGCTGATCCCGCGAGCCGGTTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCt                                                                     <  2:690535/144‑1 (MQ=255)
                                                                                   ccATCCCCCCTTCATCGTATCCGGGGGGCGTTGTTTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGGGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgcgc                                                                      <  2:627230/149‑1 (MQ=255)
                                                                                           ccTTCATTGTTTTCCGGGGGCGTTGTGTTGCCGCGCTCCTGCGGCGCTGCCGTTGGATCCCGCGGGCCGGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCAc                                                            <  1:519829/151‑1 (MQ=255)
                                                                                             ttCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGCTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTc                                                              >  1:279603/1‑147 (MQ=255)
                                                                                               cATCGTATTCCGGGGGCTTTGTGTGGCGGCTGTACTGCCGCGCGGACGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGc                                                         <  2:125753/150‑1 (MQ=255)
                                                                                                aTCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGGTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCt                                                        >  1:1370802/1‑150 (MQ=255)
                                                                                                 tCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTg                                                       >  1:378413/1‑150 (MQ=255)
                                                                                                   ttttgccGGGGGCGTTGTTTTCCCGCTGGCCGGCCGCGCTGGCGGTTGATACCGCGGCCCGGTCACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCt                                                        <  1:51426/142‑1 (MQ=255)
                                                                                                   gTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCt                                                    >  2:443435/1‑151 (MQ=255)
                                                                                                    ttgcccGGGGGCGTTTTGTTGCCTCTGCCCGGCCGCGCTGACGGCTGATCCCGCGCGCCCGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCg                                                          <  1:139248/140‑1 (MQ=255)
                                                                                                     ttccggggggCGTTTTTTTGCCTCTTTCCTGCCGCGTTGGCGGCTGATCCGGCGACCCGGTTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTg                                                   <  1:164489/145‑1 (MQ=255)
                                                                                                       tCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCg                                                 >  1:1268934/1‑150 (MQ=255)
                                                                                                       ggcGGGGGAGTTGTGTTGCCGTCATAGTGCCGCGCGGACGGCGGGTTCGGGGGGCCCGGAGCCGCCATATCCAGCGGCTGGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGc                                                  <  2:28080/147‑1 (MQ=255)
                                                                                                        cgggggggGTTGTTTTGCCGCTGTCCGGCCGCGCTGCCGGCTGATCCCGGGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGc                                                  <  2:520072/146‑1 (MQ=255)
                                                                                                          gggggCTTTGTGTTGCCGCTGTACGGCCGCGCTGCCGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAAc                                              <  2:722128/150‑1 (MQ=255)
                                                                                                          gggggCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAAcc                                             >  2:349441/1‑151 (MQ=255)
                                                                                                          gggggCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGCACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAAcc                                             >  1:143743/1‑151 (MQ=255)
                                                                                                            ggggGTTGTGTTGCCCTTTTCCTGCCCCGCTGACGCCAAATCCCGCGGCCCCGTTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCAt                                           <  2:1398051/151‑1 (MQ=255)
                                                                                                              gCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCAttt                                         >  2:711451/1‑151 (MQ=255)
                                                                                                               cGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCAtttt                                        >  2:1091856/1‑151 (MQ=255)
                                                                                                                 ttGTGTTGCCGCTGTCCTGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCAtttt                                        <  2:601960/149‑1 (MQ=255)
                                                                                                                 tggtgtTGCCCCTGTCCGGCCGCGCGGGCGGCTGATCCCGGGAGCCCGTTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCAttt                                         <  1:1095654/146‑1 (MQ=255)
                                                                                                                  tgtgTTCCGGCTGTCCGGCCGCGCGGGCGGCGGAACCCGCGAGCCGGGTGCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCAtttt                                        <  1:802176/148‑1 (MQ=255)
                                                                                                                    tgtTGCCGCTTTACCTCCGCCCTGCCGGCTGATCCCGCGGGCCCGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTc                                       <  1:1376166/147‑1 (MQ=255)
                                                                                                                    tgtTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTg                                   >  1:286866/1‑151 (MQ=255)
                                                                                                                         ccGTTGTCCTGCCGCGCTGCCGGCTGACCCCGCGACCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgc                                  <  2:361373/147‑1 (MQ=255)
                                                                                                                         ccGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgcagc                               >  1:550076/1‑150 (MQ=255)
                                                                                                                          ctcTGTCCTGCCGCGGGGCGGGCTGCTCCCGCGCGCAGGTTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgcagca                              <  1:567174/148‑1 (MQ=255)
                                                                                                                          cGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCAc                             >  2:1340529/1‑151 (MQ=255)
                                                                                                                                ccggccgcgtggcgggCTAATCCCGCGCCCGCGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTcc                         <  1:779259/135‑1 (MQ=255)
                                                                                                                                   gCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCt                     >  2:381391/1‑150 (MQ=255)
                                                                                                                                          tGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTg             >  1:601461/1‑151 (MQ=255)
                                                                                                                                          tGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTg             >  1:746785/1‑151 (MQ=255)
                                                                                                                                          ggcggggtGATCCCGCGGCCCCGTTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgca                                 <  1:1182742/124‑1 (MQ=255)
                                                                                                                                           gACGGCTGATCCCGCGGGCCCGGCACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGca                      <  1:1123062/141‑2 (MQ=255)
                                                                                                                                              ggCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGc               <  2:1248999/145‑1 (MQ=255)
                                                                                                                                              ggCTGATCCCGAGACCACGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGc               <  1:1392350/145‑1 (MQ=255)
                                                                                                                                              ggCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCgg          >  2:790546/1‑150 (MQ=255)
                                                                                                                                              ggCTAATCCCGCGGGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCg           <  1:316542/149‑1 (MQ=255)
                                                                                                                                                cTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTg        >  2:1276857/1‑150 (MQ=255)
                                                                                                                                                     aCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACg  >  2:269488/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
ACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACG  >  NC_000913/280941‑281239

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