Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,138 (A)6→7 coding (15/594 nt) bdcR → transcriptional repressor for divergent bdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,1321.A98.4% 206.0 / ‑4.1 64coding (9/594 nt)bdcRtranscriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (36/27);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTT  >  NC_000913/4474001‑4474261
                                                                                                                                    |                                                                                                                                    
aTACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaa                                                                                                                                    >  1:1105857/1‑134 (MQ=255)
aTACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaa                                                                                                                                    <  2:1105857/134‑1 (MQ=255)
          tGCCACCGAGGATGAGAACTGTCTTACCTGTACAAGCGCCCATAAGGACTCCTTGATTGATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCATTAAAAAAACAATCAcgcg                                                                                                                     >  2:1731747/1‑139 (MQ=255)
                         gaACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaaaaaaaca                                                                                                                            <  1:1120751/116‑4 (MQ=255)
                         gaACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaaaaaaaca                                                                                                                            >  2:1120751/1‑113 (MQ=255)
                               tCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGtt                                                                                                >  1:64807/1‑139 (MQ=255)
                                 ttACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTc                                                                                              >  1:1232466/1‑139 (MQ=255)
                                    ccTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaaaaaaaca                                                                                                                            <  1:1245328/105‑4 (MQ=255)
                                    ccTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaaaaaaaca                                                                                                                            >  2:1245328/1‑102 (MQ=255)
                                         aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                             >  1:1815803/1‑116 (MQ=255)
                                         aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                             <  2:1815803/116‑1 (MQ=255)
                                          aaaGCGCCCATAAGGACTACTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     <  2:1232466/139‑1 (MQ=255)
                                                           tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                        >  1:1015507/1‑87 (MQ=255)
                                                           tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                        <  2:1015507/87‑1 (MQ=255)
                                                             cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                        <  2:1228350/85‑1 (MQ=255)
                                                             cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAg                                                                                                     <  2:785855/104‑1 (MQ=255)
                                                             cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAg                                                                                                     >  1:785855/1‑104 (MQ=255)
                                                             cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     <  2:1320253/120‑1 (MQ=255)
                                                             cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     >  1:1320253/1‑120 (MQ=255)
                                                             cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑Tctctctgcggc                                                                >  2:1760584/1‑132 (MQ=255)
                                                             cTTGATTTATCATATAACATGCATTACAAAACTGTTTAAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                        >  1:1228350/1‑85 (MQ=255)
                                                                     attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑tctc                                                                        >  1:184646/1‑123 (MQ=255)
                                                                     attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑tctc                                                                        <  2:184646/123‑1 (MQ=255)
                                                                       tatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     <  1:1146786/110‑1 (MQ=255)
                                                                       tatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     >  2:1146786/1‑110 (MQ=255)
                                                                           tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCAc                                                  >  2:1613768/1‑139 (MQ=255)
                                                                           tAACATACATTCCAAAACTGTCTTCACTTTCTGTAACCAGGCTTTAAAATGCTCACTAACAAAACAATCCCGCGCTCCAGCTCGTCCCAAACGTTTCGCTCCTGAGCAGGCAA‑‑TCACCCCACCAAAAGTGCTATTTCAc                                                  >  1:679928/1‑139 (MQ=255)
                                                                               aTGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                              <  1:1613768/139‑1 (MQ=255)
                                                                                tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                             >  2:1484485/1‑139 (MQ=255)
                                                                                tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                             >  1:1688157/1‑139 (MQ=255)
                                                                                tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTATAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                             >  2:1117909/1‑139 (MQ=255)
                                                                                 gCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                              <  1:1760584/139‑1 (MQ=255)
                                                                                 gCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCtttt                                                      >  1:436168/1‑129 (MQ=255)
                                                                                 gCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCtttt                                                      <  2:436168/129‑1 (MQ=255)
                                                                                   aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCAcc                                                 >  2:309993/1‑132 (MQ=255)
                                                                                   aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCAcc                                                 <  1:309993/132‑1 (MQ=255)
                                                                                     tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                        >  1:1704761/1‑139 (MQ=255)
                                                                                         aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                    <  1:1117909/139‑1 (MQ=255)
                                                                                          aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                   >  1:1375689/1‑139 (MQ=255)
                                                                                            cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                        <  1:518425/54‑1 (MQ=255)
                                                                                            cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                        >  2:518425/1‑54 (MQ=255)
                                                                                            cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGc                                 >  1:1734173/1‑139 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                               <  2:844472/76‑1 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                               >  1:844472/1‑76 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                   >  1:1628616/1‑134 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                   <  2:1628616/134‑1 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                              >  1:1446722/1‑139 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                              >  1:1444413/1‑139 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaa                                                                                                                                    <  2:969917/39‑1 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaa                                                                                                                                    >  1:969917/1‑39 (MQ=255)
                                                                                               ttttAACTTTCTGTCAACAGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGGTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTc                             >  2:794107/1‑139 (MQ=255)
                                                                                                       ttCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTg                      >  2:936491/1‑139 (MQ=255)
                                                                                                         cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                               >  1:1353224/1‑66 (MQ=255)
                                                                                                         cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                               <  2:1353224/66‑1 (MQ=255)
                                                                                                              aaCAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt               >  2:1757405/1‑139 (MQ=255)
                                                                                                              aaCAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt               <  1:1484485/139‑1 (MQ=255)
                                                                                                               aCAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGtt              <  2:1704761/139‑1 (MQ=255)
                                                                                                                     tttGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAt        <  2:1734173/139‑1 (MQ=255)
                                                                                                                     tttGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAt        <  1:936491/139‑1 (MQ=255)
                                                                                                                      ttGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa             <  1:520996/133‑1 (MQ=255)
                                                                                                                      ttGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa             >  2:520996/1‑133 (MQ=255)
                                                                                                                      ttGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAtt       >  1:1817000/1‑139 (MQ=255)
                                                                                                                      ttGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAtt       >  2:576929/1‑139 (MQ=255)
                                                                                                                           aTGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCtt  <  2:1817000/139‑1 (MQ=255)
                                                                                                                             ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                    >  2:755856/1‑41 (MQ=255)
                                                                                                                             ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                    <  1:755856/41‑1 (MQ=255)
                                                                                                                             ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     <  1:1638030/56‑1 (MQ=255)
                                                                                                                             ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     >  2:1638030/1‑56 (MQ=255)
                                                                                                                             ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCttttt                                                     <  1:241925/86‑1 (MQ=255)
                                                                                                                             ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCttttt                                                     >  2:241925/1‑86 (MQ=255)
                                                                                                                             ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCtt  >  2:811878/1‑137 (MQ=255)
                                                                                                                             ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCtt  <  1:811878/137‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
ATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAA‑‑TCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTT  >  NC_000913/4474001‑4474261

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