Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,0680TA21.3% 229.2 / 55.8 89pseudogene (348/504 nt)insB1IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base A (10/9);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATC  >  NC_000913/4518935‑4519201
                                                                                                                                     |                                                                                                                                     
cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                                                                                 <  2:474480/139‑1 (MQ=14)
 gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                                                                                 <  1:70828/139‑1 (MQ=255)
  cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCa                                                                                                                                >  2:1392170/1‑139 (MQ=255)
    cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                                                                                 >  2:1270166/1‑136 (MQ=255)
    cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                                                                                 <  1:1270166/136‑1 (MQ=255)
                 cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCaa                                                                                                                               >  2:1150419/1‑125 (MQ=255)
                 cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCaa                                                                                                                               <  1:1150419/125‑1 (MQ=255)
                 cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAgcg                                                                                                                >  2:1447187/1‑139 (MQ=11)
                  gCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGTTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAgcg                                                                                                                <  2:1055418/139‑1 (MQ=255)
                       atGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATATATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                                                                                          <  2:780090/138‑1 (MQ=11)
                         ggCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAg                                                                                                         <  1:1392170/139‑1 (MQ=255)
                         ggCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAg                                                                                                         <  2:15897/139‑1 (MQ=255)
                             acgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggta                                                                                                    <  2:1194532/136‑4 (MQ=11)
                              cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCg                                                                                                                            >  2:1417102/1‑115 (MQ=255)
                              cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCg                                                                                                                            <  1:1417102/115‑1 (MQ=255)
                                  ggAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtataacc                                                                                               >  2:391415/1‑131 (MQ=11)
                                  ggAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATaacc                                                                                               >  2:130611/1‑137 (MQ=11)
                                    aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                                                                                              >  2:1373476/1‑139 (MQ=255)
                                    aGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATaacctg                                                                                             >  2:469154/1‑139 (MQ=11)
                                        tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGt                                                                                                                           >  1:1078235/1‑106 (MQ=255)
                                        tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGt                                                                                                                           <  2:1078235/106‑1 (MQ=255)
                                              gAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGc                                                                                                        >  1:990696/1‑118 (MQ=14)
                                              gAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGc                                                                                                        <  2:990696/118‑1 (MQ=11)
                                                   tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGc                                                                                                                             <  2:663640/93‑1 (MQ=255)
                                                   tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGc                                                                                                                             >  1:663640/1‑93 (MQ=255)
                                                   tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGa                                                                                                      >  2:1185462/1‑116 (MQ=255)
                                                   tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGa                                                                                                      <  1:1185462/116‑1 (MQ=255)
                                                     ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          <  2:1161168/110‑1 (MQ=255)
                                                     ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          >  1:1161168/1‑110 (MQ=255)
                                                      ctgtggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                                                                               >  2:1098316/5‑88 (MQ=17)
                                                      ctgtggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                                                                               <  1:1098316/84‑1 (MQ=17)
                                                          ggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                                                                                          >  1:66323/1‑105 (MQ=11)
                                                          ggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                                                                                          <  2:66323/105‑1 (MQ=11)
                                                           gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCaa                                                                       <  2:636671/139‑1 (MQ=255)
                                                           gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGc                                                                                                        >  1:1463104/1‑106 (MQ=11)
                                                           gCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGc                                                                                                        <  2:1463104/106‑1 (MQ=11)
                                                              tttGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                                                                               >  1:1018789/1‑80 (MQ=17)
                                                              tttGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                                                                               <  2:1018789/80‑1 (MQ=17)
                                                                tGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                                                                                 <  2:69844/76‑1 (MQ=255)
                                                                tGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                                                                                 >  1:69844/1‑76 (MQ=255)
                                                                tGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                                                                                 <  2:668075/76‑1 (MQ=255)
                                                                tGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                                                                                 >  1:668075/1‑76 (MQ=255)
                                                                  aGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTggg                                                                <  2:340647/139‑1 (MQ=255)
                                                                   ggTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGa                                                               <  1:450183/139‑1 (MQ=255)
                                                                       gTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          <  2:1299566/92‑1 (MQ=255)
                                                                       gTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          >  1:1299566/1‑92 (MQ=255)
                                                                          gTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATa                                                                                                  <  1:632742/97‑1 (MQ=255)
                                                                          gTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATa                                                                                                  >  2:632742/1‑97 (MQ=255)
                                                                                    cGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGc                                                                    >  1:978323/1‑117 (MQ=255)
                                                                                    cGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGc                                                                    <  2:978323/117‑1 (MQ=255)
                                                                                    cGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGtt                                              >  2:658996/1‑139 (MQ=255)
                                                                                         ggctggcCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCaa                                         >  2:260192/1‑139 (MQ=255)
                                                                                            gggcCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACTTTAGCAGGAAGCGTTACACTCAGCGCATTGAGCGACATCATGTGAATCTGAGACAATATCTGGCATGGCTGGGACGGAAGTTACTGTCGGTCTCAAAAt                                      <  1:94337/138‑1 (MQ=18)
                                                                                              gcCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGtt                                                                                                                                       <  1:70255/40‑1 (MQ=255)
                                                                                              gcCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGtt                                                                                                                                       >  2:70255/1‑40 (MQ=255)
                                                                                                      aTGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGtt                            <  1:658996/139‑1 (MQ=255)
                                                                                                             aCGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGtt                                                                                                                          >  2:1421697/1‑38 (MQ=255)
                                                                                                             aCGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGtt                                                                                                                          <  1:1421697/38‑1 (MQ=255)
                                                                                                             aCGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGAc                     <  1:260192/139‑1 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          <  2:344315/52‑1 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          >  1:344315/1‑52 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaat                                                                                                <  1:296456/62‑1 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaat                                                                                                >  2:296456/1‑62 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCAt                        >  1:513245/1‑134 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCAt                        <  2:513245/134‑1 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACaa                   >  2:1351265/1‑139 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACaa                   >  2:658940/1‑139 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACaa                   >  1:1065919/1‑139 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACaa                   >  1:1318475/1‑139 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACaa                   >  1:1359455/1‑139 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACa                    >  2:578165/1‑138 (MQ=255)
                                                                                                               gCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACa                    <  1:578165/138‑1 (MQ=255)
                                                                                                                    aaGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCa              >  2:707362/1‑139 (MQ=255)
                                                                                                                    aaGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCa              >  2:814346/1‑139 (MQ=255)
                                                                                                                     aGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCAt             <  2:563899/139‑1 (MQ=255)
                                                                                                                         aaaGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCggg         >  1:1329444/1‑139 (MQ=255)
                                                                                                                         aaaGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCggg         >  2:734507/1‑139 (MQ=255)
                                                                                                                          aaGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGc        <  1:658940/139‑1 (MQ=255)
                                                                                                                          aaGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGc        <  1:1373476/139‑1 (MQ=255)
                                                                                                                          aaGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGc        <  2:1065919/139‑1 (MQ=255)
                                                                                                                              tGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTg                           <  1:939058/116‑1 (MQ=255)
                                                                                                                              tGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTg                           >  2:939058/1‑116 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          >  2:66735/1‑35 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          <  2:1127460/35‑1 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          >  1:1127460/1‑35 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          <  2:890735/35‑1 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          >  1:890735/1‑35 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                          <  1:66735/35‑1 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTgg                                >  1:1279346/1‑109 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTgg                                <  2:1279346/109‑1 (MQ=255)
                                                                                                                                cACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATc  >  2:1171070/1‑139 (MQ=255)
                                                                                                                                     tATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACg                                                             <  1:1004383/75‑1 (MQ=255)
                                                                                                                                     tATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACg                                                             >  2:1004383/1‑75 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATC  >  NC_000913/4518935‑4519201

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