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

GCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGT  >  NC_000913/4473995‑4474264
                                                                                                                                          |                                                                                                                                     
gCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑a                                                                                                                                      >  1:844507/1‑139 (MQ=255)
 cACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCTCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTGCTGTCAACAGGTTTTGTAATGGTCACT‑aa                                                                                                                                     >  1:1037473/1‑139 (MQ=255)
                  ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGtt                                                                                                                    >  2:374840/1‑139 (MQ=255)
                            tGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgtcgtc                                                                                                          >  2:1141435/1‑139 (MQ=255)
                                             gTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACCTGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAGAAAACCATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCt                                                                                         <  1:1026667/139‑1 (MQ=255)
                                               aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAgg                                                                                                               <  1:645053/115‑1 (MQ=255)
                                               aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATCACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAgg                                                                                                               >  2:645053/1‑115 (MQ=255)
                                                      cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTcg                                                                                                                        <  1:209359/99‑1 (MQ=255)
                                                      cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTcg                                                                                                                        >  2:209359/1‑99 (MQ=255)
                                                        cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAt                                                                              <  2:844507/139‑1 (MQ=255)
                                                        cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAt                                                                              >  2:1213884/1‑139 (MQ=255)
                                                                 tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAtctc                                                                           >  1:797600/1‑133 (MQ=255)
                                                                 tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAtctc                                                                           <  2:797600/133‑1 (MQ=255)
                                                                   cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCg                                                                                              <  2:1149158/112‑1 (MQ=255)
                                                                   cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCg                                                                                              >  1:1149158/1‑112 (MQ=255)
                                                                   cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaa                                                                   >  2:352894/1‑139 (MQ=255)
                                                                        tttATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTg                                                              >  1:163651/1‑139 (MQ=255)
                                                                           attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  1:535477/1‑112 (MQ=255)
                                                                           attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  2:535477/112‑1 (MQ=255)
                                                                              atGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                        <  1:1213884/139‑1 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAAtctc                                                                                                                         <  1:230185/71‑1 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAAtctc                                                                                                                         >  2:230185/1‑71 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  2:739856/1‑106 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:739856/106‑1 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAc                                                     >  1:1532362/1‑139 (MQ=255)
                                                                                 tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTAGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAc                                                     >  2:546968/1‑139 (MQ=255)
                                                                                  aaCATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              >  2:1153300/1‑81 (MQ=255)
                                                                                  aaCATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              <  1:1153300/81‑1 (MQ=255)
                                                                                      tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                >  2:1430033/1‑139 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTcg                                                                                                                        <  1:612154/64‑1 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTcg                                                                                                                        >  2:612154/1‑64 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                >  1:631212/1‑88 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                <  2:631212/88‑1 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                             >  2:69313/1‑139 (MQ=255)
                                                                                         aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                             >  1:899968/1‑139 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              >  1:427323/1‑72 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              <  2:427323/72‑1 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 >  2:54375/1‑133 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 <  1:54375/133‑1 (MQ=255)
                                                                                           tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGACAAAAGTGCTTTTTCACCAAAAAGGtt                                           >  2:437026/1‑139 (MQ=255)
                                                                                               aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                       <  2:1532362/139‑1 (MQ=255)
                                                                                               aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                       <  1:352894/139‑1 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              >  2:1365499/1‑67 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              <  1:1365499/67‑1 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  1:250675/1‑91 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  2:250675/91‑1 (MQ=255)
                                                                                                aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGACAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                      <  1:437026/139‑1 (MQ=255)
                                                                                                  cTGTTTTAACTTTCTGTCAACAGGTTTTGTTATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaaa                                               <  2:918917/128‑1 (MQ=255)
                                                                                                  cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAATTGCTTTTTCACCaaaaa                                               >  1:918917/1‑128 (MQ=255)
                                                                                                     ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAAtctc                                                                                                                         <  2:736615/51‑1 (MQ=255)
                                                                                                     ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                 >  1:1029412/1‑139 (MQ=255)
                                                                                                     ttttAACTTTCTGTCAACAGGATTTGTAATGGTCACTAAA‑AAAACAAtctc                                                                                                                         >  1:736615/1‑51 (MQ=255)
                                                                                                       ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                <  1:321527/122‑1 (MQ=255)
                                                                                                       ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                >  2:321527/1‑122 (MQ=255)
                                                                                                       ttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCa                               >  2:143264/1‑139 (MQ=255)
                                                                                                         aaCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAgt                             <  1:775677/139‑1 (MQ=255)
                                                                                                             ttCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                     >  2:57444/1‑63 (MQ=255)
                                                                                                             ttCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                     <  1:57444/63‑1 (MQ=255)
                                                                                                               cTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAgg                                                                                                               >  2:23618/1‑51 (MQ=255)
                                                                                                               cTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAgg                                                                                                               <  1:23618/51‑1 (MQ=255)
                                                                                                               cTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  2:869065/1‑76 (MQ=255)
                                                                                                               cTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:869065/76‑1 (MQ=255)
                                                                                                                    aaCAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt                  <  1:69313/139‑1 (MQ=255)
                                                                                                                      cAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa                <  2:972160/139‑1 (MQ=255)
                                                                                                                       aGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTAc               <  2:163651/139‑1 (MQ=255)
                                                                                                                         gTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTg             >  2:15727/1‑139 (MQ=255)
                                                                                                                           tttGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAt           <  2:1029412/139‑1 (MQ=255)
                                                                                                                                aaTGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt      <  1:15727/139‑1 (MQ=255)
                                                                                                                                   ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCt                                   <  2:1132398/107‑1 (MQ=255)
                                                                                                                                   ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCt                                   >  1:1132398/1‑107 (MQ=255)
                                                                                                                                   ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt      <  2:441892/136‑1 (MQ=255)
                                                                                                                                   ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt      >  1:441892/1‑136 (MQ=255)
                                                                                                                                    gTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGTCGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGt                                                               <  1:212866/78‑1 (MQ=255)
                                                                                                                                    gTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                           <  2:1160920/50‑1 (MQ=255)
                                                                                                                                    gTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                           >  1:1160920/1‑50 (MQ=255)
                                                                                                                                    gTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGt                                                               >  2:212866/1‑78 (MQ=255)
                                                                                                                                    gTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGt  >  2:237028/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                     
GCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGT  >  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: ‑