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.7% 241.9 / ‑3.7 75coding (9/594 nt)bdcRtranscriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (39/35);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGT  >  NC_000913/4473995‑4474264
                                                                                                                                          |                                                                                                                                     
gCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTa                                                                                                                                      >  1:844507/1‑139 (MQ=255)
 cACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCTCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTGCTGTCAACAGGTTTTGTAATGGTCACTaa                                                                                                                                     >  1:1037473/1‑139 (MQ=255)
                  ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGtt                                                                                                                    >  2:374840/1‑139 (MQ=255)
                            tGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgtcgtc                                                                                                          >  2:1141435/1‑139 (MQ=255)
                                             gTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACCTGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAGAAAACCATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCt                                                                                         <  1:1026667/139‑1 (MQ=255)
                                               aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAgg                                                                                                               <  1:645053/115‑1 (MQ=255)
                                               aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATCACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAgg                                                                                                               >  2:645053/1‑115 (MQ=255)
                                                      cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTcg                                                                                                                        <  1:209359/99‑1 (MQ=255)
                                                      cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTcg                                                                                                                        >  2:209359/1‑99 (MQ=255)
                                                        cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAt                                                                              <  2:844507/139‑1 (MQ=255)
                                                        cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAt                                                                              >  2:1213884/1‑139 (MQ=255)
                                                                 tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAtctc                                                                           >  1:797600/1‑133 (MQ=255)
                                                                 tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAtctc                                                                           <  2:797600/133‑1 (MQ=255)
                                                                   cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCg                                                                                              <  2:1149158/112‑1 (MQ=255)
                                                                   cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCg                                                                                              >  1:1149158/1‑112 (MQ=255)
                                                                   cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaa                                                                   >  2:352894/1‑139 (MQ=255)
                                                                        tttATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTg                                                              >  1:163651/1‑139 (MQ=255)
                                                                           attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  1:535477/1‑112 (MQ=255)
                                                                           attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  2:535477/112‑1 (MQ=255)
                                                                              atGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                        <  1:1213884/139‑1 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                         <  1:230185/71‑1 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                         >  2:230185/1‑71 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  2:739856/1‑106 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:739856/106‑1 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAc                                                     >  1:1532362/1‑139 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTAGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAc                                                     >  2:546968/1‑139 (MQ=255)
                                                                                  aaCATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              >  2:1153300/1‑81 (MQ=255)
                                                                                  aaCATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              <  1:1153300/81‑1 (MQ=255)
                                                                                      tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                >  2:1430033/1‑139 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTcg                                                                                                                        <  1:612154/64‑1 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTcg                                                                                                                        >  2:612154/1‑64 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                >  1:631212/1‑88 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                <  2:631212/88‑1 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                             >  2:69313/1‑139 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                             >  1:899968/1‑139 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              >  1:427323/1‑72 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              <  2:427323/72‑1 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 >  2:54375/1‑133 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 <  1:54375/133‑1 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGACAAAAGTGCTTTTTCACCAAAAAGGtt                                           >  2:437026/1‑139 (MQ=255)
                                                                                               aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                       <  2:1532362/139‑1 (MQ=255)
                                                                                               aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                       <  1:352894/139‑1 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              >  2:1365499/1‑67 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                              <  1:1365499/67‑1 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  1:250675/1‑91 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  2:250675/91‑1 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGACAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                      <  1:437026/139‑1 (MQ=255)
                                                                                                  cTGTTTTAACTTTCTGTCAACAGGTTTTGTTATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaaa                                               <  2:918917/128‑1 (MQ=255)
                                                                                                  cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAATTGCTTTTTCACCaaaaa                                               >  1:918917/1‑128 (MQ=255)
                                                                                                     ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                         <  2:736615/51‑1 (MQ=255)
                                                                                                     ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                 >  1:1029412/1‑139 (MQ=255)
                                                                                                     ttttAACTTTCTGTCAACAGGATTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                         >  1:736615/1‑51 (MQ=255)
                                                                                                       ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                <  1:321527/122‑1 (MQ=255)
                                                                                                       ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                >  2:321527/1‑122 (MQ=255)
                                                                                                       ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCa                               >  2:143264/1‑139 (MQ=255)
                                                                                                         aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAgt                             <  1:775677/139‑1 (MQ=255)
                                                                                                             ttCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                     >  2:57444/1‑63 (MQ=255)
                                                                                                             ttCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                     <  1:57444/63‑1 (MQ=255)
                                                                                                               cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAgg                                                                                                               >  2:23618/1‑51 (MQ=255)
                                                                                                               cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAgg                                                                                                               <  1:23618/51‑1 (MQ=255)
                                                                                                               cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  2:869065/1‑76 (MQ=255)
                                                                                                               cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:869065/76‑1 (MQ=255)
                                                                                                                    aaCAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt                  <  1:69313/139‑1 (MQ=255)
                                                                                                                      cAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa                <  2:972160/139‑1 (MQ=255)
                                                                                                                       aGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTAc               <  2:163651/139‑1 (MQ=255)
                                                                                                                         gTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTg             >  2:15727/1‑139 (MQ=255)
                                                                                                                           tttGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAt           <  2:1029412/139‑1 (MQ=255)
                                                                                                                                aaTGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt      <  1:15727/139‑1 (MQ=255)
                                                                                                                                   ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCt                                   <  2:1132398/107‑1 (MQ=255)
                                                                                                                                   ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCt                                   >  1:1132398/1‑107 (MQ=255)
                                                                                                                                   ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt      <  2:441892/136‑1 (MQ=255)
                                                                                                                                   ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt      >  1:441892/1‑136 (MQ=255)
                                                                                                                                    gTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGTCGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGt                                                               <  1:212866/78‑1 (MQ=255)
                                                                                                                                    gTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                           <  2:1160920/50‑1 (MQ=255)
                                                                                                                                    gTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                           >  1:1160920/1‑50 (MQ=255)
                                                                                                                                    gTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGt                                                               >  2:212866/1‑78 (MQ=255)
                                                                                                                                    gTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGt  >  2:237028/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                     
GCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGT  >  NC_000913/4473995‑4474264

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