Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,112,787 C→T intergenic (+32/‑180) yiiR → / → yiiS DUF805 family putative inner membrane protein/UPF0381 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,112,7870CT95.5% 266.6 / 2.8 89intergenic (+32/‑180)yiiR/yiiSDUF805 family putative inner membrane protein/UPF0381 family protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (44/41);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

TGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTCCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTTGTTTGTTCTTCCCTTAACGGCGCA  >  NC_000913/4112668‑4112917
                                                                                                                       |                                                                                                                                  
tGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACg                                                                                                                 >  1:803647/1‑139 (MQ=255)
    tttGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGc                                                                                                             >  2:439480/1‑139 (MQ=255)
    tttGTCGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACg                                                                                                                 >  2:1180029/1‑135 (MQ=255)
    tttGTCGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACg                                                                                                                 <  1:1180029/135‑1 (MQ=255)
     ttgttgGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCa                                                                                                            >  1:1190944/1‑139 (MQ=255)
               acaAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCa                                                                                                                   >  1:595811/1‑122 (MQ=255)
               acaAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCa                                                                                                                   <  2:595811/122‑1 (MQ=255)
                  aGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACg                                                                                                                 >  1:1358429/1‑121 (MQ=255)
                  aGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACg                                                                                                                 <  2:1358429/121‑1 (MQ=255)
                     cGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGa                                                                                            >  1:1448345/1‑139 (MQ=255)
                               tatGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCt                                                                                  >  2:625795/1‑139 (MQ=255)
                                atGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTg                                                                                 <  1:625795/139‑1 (MQ=255)
                                atGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTg                                                                                 <  1:439480/139‑1 (MQ=255)
                                  ggTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGaa                                                                               >  1:582545/1‑139 (MQ=255)
                                  ggTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTGATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGaa                                                                               >  2:1124397/1‑139 (MQ=255)
                                   gTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAAc                                                                              >  2:475275/1‑139 (MQ=255)
                                      aaGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAActct                                                                           <  2:803647/139‑1 (MQ=255)
                                      aaGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAActct                                                                           <  1:1724176/139‑1 (MQ=255)
                                       aGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATTTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  1:1069529/139‑1 (MQ=255)
                                        gATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTc                                                                         <  2:582545/139‑1 (MQ=255)
                                              cAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGaaaa                                                                                         >  2:1415652/1‑117 (MQ=255)
                                              cAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGaaaa                                                                                         <  1:1415652/117‑1 (MQ=255)
                                                ggATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACg                                                                                                                 <  2:1683588/91‑1 (MQ=255)
                                                ggATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACg                                                                                                                 >  1:1683588/1‑91 (MQ=255)
                                                    gtgAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          >  1:1518701/1‑126 (MQ=255)
                                                    gtgAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  2:1518701/126‑1 (MQ=255)
                                                       aaGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  2:486028/123‑1 (MQ=255)
                                                       aaGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          >  1:486028/1‑123 (MQ=255)
                                                           ataAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCaaaa                                                      >  1:938777/1‑139 (MQ=255)
                                                               aGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          >  2:446485/1‑115 (MQ=255)
                                                               aGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  1:446485/115‑1 (MQ=255)
                                                                   gACAATAAATCTAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAACAAAATGGt                                              <  1:1124397/139‑1 (MQ=255)
                                                                   gACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGt                                              <  2:1448345/139‑1 (MQ=255)
                                                                    aCAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCaaa                                                                                                          >  1:1195488/1‑78 (MQ=255)
                                                                    aCAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCaaa                                                                                                          <  2:1195488/78‑1 (MQ=255)
                                                                    aCAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGt                                              <  2:306403/138‑1 (MQ=255)
                                                                    aCAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGt                                              >  1:306403/1‑138 (MQ=255)
                                                                    aCAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTa                                             <  1:475275/139‑1 (MQ=255)
                                                                    aCAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTa                                             >  1:943180/1‑139 (MQ=255)
                                                                         aaaTCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCa                                        <  2:1190944/139‑1 (MQ=255)
                                                                           aTCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATa                                      >  1:722657/1‑139 (MQ=255)
                                                                                    ttAACCTTTAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTa                             >  1:26294/1‑139 (MQ=255)
                                                                                          tttAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGt                                              >  2:938036/1‑116 (MQ=255)
                                                                                          tttAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGt                                              <  1:938036/116‑1 (MQ=255)
                                                                                           ttAATATATCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACttgttt                      >  1:1513774/1‑139 (MQ=255)
                                                                                                atatCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTTGTTTGttct                 <  2:943180/139‑1 (MQ=255)
                                                                                                  atCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATg                                                                                                              <  2:1445373/44‑1 (MQ=255)
                                                                                                  atCTCTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATg                                                                                                              >  1:1445373/1‑44 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCa                                                                                                            >  2:1428190/1‑45 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCa                                                                                                            <  1:1428190/45‑1 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGa                                                                                <  1:1021497/73‑1 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGa                                                                                >  2:1021497/1‑73 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTTGTTTGttct                 <  2:1123618/136‑1 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTTGTTTGttct                 >  1:1123618/1‑136 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTTGTTTGTTCTTcc              >  1:814220/1‑139 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTCCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCtt                                                                    <  2:77099/85‑1 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTCCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCtt                                                                    >  1:77099/1‑85 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTCCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCAt                                       <  1:548180/114‑1 (MQ=255)
                                                                                                   tctcTTACTGCTATTTTTTTCCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCAt                                       >  2:548180/1‑114 (MQ=255)
                                                                                                    ctctTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAActct                                                                           <  1:203870/77‑1 (MQ=255)
                                                                                                    ctctTACTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAActct                                                                           >  2:203870/1‑77 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTg                                                                                 >  2:148961/1‑66 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTg                                                                                 <  1:148961/66‑1 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAActct                                                                           <  1:868309/72‑1 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAActct                                                                           >  2:868309/1‑72 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          >  2:1546425/1‑73 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  2:254193/73‑1 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  1:1546425/73‑1 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          >  1:254193/1‑73 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCa                                                               >  1:135871/1‑84 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCa                                                               <  2:135871/84‑1 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATc                                                          <  2:1664782/89‑1 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATc                                                          >  1:1664782/1‑89 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCaaaaaaatggt                                               <  1:573185/100‑5 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCaaaaaaatggt                                               >  2:573185/1‑96 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTtgtttgtt                   >  1:1425049/1‑128 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTtgtttgtt                   <  2:1425049/128‑1 (MQ=255)
                                                                                                         aCTGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTTGTTTGTTCTTCCCTTAAc        >  1:748353/1‑139 (MQ=255)
                                                                                                           tGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  2:1505973/71‑1 (MQ=255)
                                                                                                           tGCTATTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          >  1:1505973/1‑71 (MQ=255)
                                                                                                               aTTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          >  2:500236/1‑67 (MQ=255)
                                                                                                               aTTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  1:500236/67‑1 (MQ=255)
                                                                                                               aTTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          <  1:1618723/67‑1 (MQ=255)
                                                                                                               aTTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCtt                                                                          >  2:1618723/1‑67 (MQ=255)
                                                                                                               aTTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTtgtttgtt                   >  1:874739/1‑122 (MQ=255)
                                                                                                               aTTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTtgtttgtt                   <  2:874739/122‑1 (MQ=255)
                                                                                                               aTTTTTTTTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTTGTTTGTTCTTCCCTTAACGGCGCa  <  2:748353/139‑1 (MQ=255)
                                                                                                                 ttttttTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCaaaaaaaa                                                  <  1:1663496/89‑1 (MQ=255)
                                                                                                                 ttttttTCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCaaaaaaaa                                                  >  2:1663496/1‑89 (MQ=255)
                                                                                                                       |                                                                                                                                  
TGCGTTTGTTGGCACACAAGGCGAAAATAAATATGGTAAAGATACTCAGGATGTGAAGTATAAAGCCGACAATAAATCAAGTAATTAACCTTTAATATATCTCTTACTGCTATTTTTTTCCGGGATTGTTATTTCCACGATGCAAATGACAATACCCGGAAAATCCTTCTGAACTCTTCACCTTAAGCAATATCAAAAAAAATGGTAATCCATAAGATCATTACTTGTTTGTTCTTCCCTTAACGGCGCA  >  NC_000913/4112668‑4112917

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