Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,023,245 A→G 6.2% E170E (GAA→GAG tatC → TatABCE protein translocation system subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,023,2450AG 6.2% 198.3 / 2.2 66E170E (GAA→GAGtatCTatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base G (2/2);  total (37/29)
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

GCATGGCATTCGCCTACTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTG  >  NC_000913/4023097‑4023389
                                                                                                                                                    |                                                                                                                                                
gcatggcatTCGCCTACTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGaa                                                                                                                                                  <  4:265089/149‑1 (MQ=255)
gcatggcatTCGCCTACTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGaa                                                                                                                                                  <  1:214139/149‑1 (MQ=255)
gcatggcatTCGCCTACTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGaa                                                                                                                                                  <  7:267858/149‑1 (MQ=255)
    ggcatTCGCCTACTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGTCGATGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGc                                                                                                                                              >  2:69543/1‑149 (MQ=255)
     gcatTCGCCTACTTTGGGGGCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGGAGGTATCCACCGACATCTCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGc                                                                                                                                              <  2:264347/148‑1 (MQ=255)
           gCCTACTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAg                                                                                                                                        <  2:253370/148‑1 (MQ=255)
              tACTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAAt                                                                                                                                    >  2:77189/1‑149 (MQ=255)
              tACTTTGTGGGCTTTCGGCTGGCATTTGGCTTGCGTGCCAATACCGCGCTGGTAGGGGTGCAGGTATCCACCGTGTTCGACCTCTATTTAAGCTGCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTGCTTTGAAGTGACGGTAGCAAt                                                                                                                                    <  5:26319/149‑1 (MQ=255)
                cTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAAt                                                                                                                                    <  4:140256/147‑1 (MQ=255)
                 tttGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATtgt                                                                                                                                 >  2:182037/1‑149 (MQ=255)
                    gtgGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCTTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGATGTGCCGGTAGCAATTGtgct                                                                                                                              >  6:273366/1‑149 (MQ=255)
                    gtgGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGtgc                                                                                                                               <  6:105614/148‑1 (MQ=255)
                     tgGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGtgctg                                                                                                                             >  7:189051/1‑149 (MQ=255)
                        tCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGctgctg                                                                                                                          >  7:301071/1‑149 (MQ=255)
                          tttCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCtactgtg                                                                                                                        >  6:161913/1‑149 (MQ=255)
                            tCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCtgctgtgct                                                                                                                      >  8:142099/1‑149 (MQ=255)
                            tCCGCTGGAATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCAACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTACTTTGAAGTGCCGGTAGCAATTGTGCtgctgtgct                                                                                                                      >  3:96352/1‑149 (MQ=255)
                              cGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCtgctgtgctg                                                                                                                     <  4:268560/148‑1 (MQ=255)
                                cTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGAt                                                                                                                  <  1:103729/149‑1 (MQ=255)
                                   gCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATggg                                                                                                               <  8:301071/149‑1 (MQ=255)
                                    cATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATgggg                                                                                                              <  6:198801/149‑1 (MQ=255)
                                        tGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGAtt                                                                                                           <  7:142099/148‑1 (MQ=255)
                                            ttCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATACACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAGGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTc                                                                                                      >  6:241205/1‑149 (MQ=255)
                                                ttGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGCCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGc                                                                                                    >  2:148323/1‑147 (MQ=255)
                                                  gCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTACTGTGCTGGATGGGGATTACCTCGCCag                                                                                                 >  3:120662/1‑148 (MQ=255)
                                                   ccAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCag                                                                                                 >  8:20887/1‑147 (MQ=255)
                                                    cAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGAGTTTGGTGTCTCCTTTGAAGTTCCGGTAGCAATTGTGCTGCTGTGCTGGAGGGGGATTACCTCGCCag                                                                                                 >  6:61219/1‑146 (MQ=255)
                                                    cAATACCGCGCCGGAAGGGGTGCAGGTATCAACCGACATCGACAGCGATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTATCCTTTGAAATGCCTGTAGCAATTGTGCTGCTGTGATGGATGGGGATTACCTCGCCag                                                                                                 >  4:201971/1‑146 (MQ=255)
                                                         ccGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACtt                                                                                          <  8:189051/148‑1 (MQ=255)
                                                          cgcgCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGCCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTAc                                                                                        <  1:148323/149‑1 (MQ=255)
                                                                ggAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCGTCGGTATGGCGCTGTTTATGGCGTTTGGTGGCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCTGAAGCCTTACGCata                                                                                   >  2:7173/1‑146 (MQ=255)
                                                                     gggTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGc                                                                             <  4:181809/149‑1 (MQ=255)
                                                                           aGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTAt                                                                       >  3:86593/1‑149 (MQ=255)
                                                                             gTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTAtgt                                                                     >  5:48958/1‑149 (MQ=255)
                                                                                tCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAGGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCt                                                                  >  7:304837/1‑149 (MQ=255)
                                                                                   aCCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCtggt                                                               >  7:283280/1‑149 (MQ=255)
                                                                                    ccGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCGGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTCGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCttgtt                                                              >  5:305568/1‑149 (MQ=255)
                                                                                          tCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCa                                                        >  4:298138/1‑149 (MQ=255)
                                                                                          tCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGAGACTCCTTTGAAGTGCAGGTAGCAATTGTGCTGCTGTGCTGGAAGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCa                                                        >  2:242504/1‑149 (MQ=255)
                                                                                             ccAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTc                                                     >  7:304363/1‑149 (MQ=255)
                                                                                             ccAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTc                                                     <  7:20887/149‑1 (MQ=255)
                                                                                             ccAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTACTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTc                                                     <  4:120662/149‑1 (MQ=255)
                                                                                               aGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGt                                                   >  1:161636/1‑149 (MQ=255)
                                                                                                    tttAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCg                                              <  8:283280/149‑1 (MQ=255)
                                                                                                     ttAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTAt                                                                       >  8:77843/1‑123 (MQ=255)
                                                                                                     ttAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTAt                                                                       <  7:77843/123‑1 (MQ=255)
                                                                                                          cTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGt                                        >  8:125654/1‑149 (MQ=255)
                                                                                                           ttCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGtt                                       >  2:31998/1‑149 (MQ=255)
                                                                                                           ttCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGGATGTGCTGGTTGGTGCATTCGTTGTCGGGATGtt                                       >  3:211070/1‑149 (MQ=255)
                                                                                                             cGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGtt                                       <  1:31998/147‑1 (MQ=255)
                                                                                                             cGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTg                                      <  3:298138/148‑1 (MQ=255)
                                                                                                              gTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGc                                     <  1:77189/148‑1 (MQ=255)
                                                                                                                     cgcTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGAcgccg                             >  4:239001/1‑149 (MQ=255)
                                                                                                                          tttATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGATGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGccgccg                          >  5:304320/1‑147 (MQ=255)
                                                                                                                                gCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCTAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGtctt                  >  6:25199/1‑149 (MQ=255)
                                                                                                                                gCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCTGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGtct                   <  3:79519/148‑1 (MQ=255)
                                                                                                                                gCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTACTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGtct                   <  5:161913/148‑1 (MQ=255)
                                                                                                                                  gTTGGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGAGGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGtcttct                <  3:201971/149‑1 (MQ=255)
                                                                                                                                        tgtCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCaaaaaa                                                                                >  4:68706/1‑79 (MQ=255)
                                                                                                                                        tgtCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCaaaaaa                                                                                <  3:68706/79‑1 (MQ=255)
                                                                                                                                          tctcCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACg        >  3:170013/1‑149 (MQ=255)
                                                                                                                                              catagAAGTGCCGGTAGCAATTGTCCTCCTGTGCTGGAAGGGGATTACCTCGCCAGAAGACTTAAGAAAAAAACGACCGTATGTGCTCGTTCGAGCATACTTTGTCGGCATGTTACTCACGCCGCAGGATGACTACTCGCAAACGCTGt    >  4:282847/5‑149 (MQ=255)
                                                                                                                                              cTTTGCAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCgg                         >  2:258971/1‑128 (MQ=255)
                                                                                                                                              cTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGa                                           <  1:171169/110‑1 (MQ=255)
                                                                                                                                              cTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCgg                         <  1:258971/128‑1 (MQ=255)
                                                                                                                                                 tGAGGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGtt   <  5:241205/147‑1 (MQ=255)
                                                                                                                                                 tGAGGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTg  <  8:304837/148‑1 (MQ=255)
                                                                                                                                                  gAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCtggttggt                                                           <  2:31308/90‑1 (MQ=255)
                                                                                                                                                  gAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCtggttggt                                                           >  1:31308/1‑90 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCATGGCATTCGCCTACTTTGTGGTCTTTCCGCTGGCATTTGGCTTCCTTGCCAATACCGCGCCGGAAGGGGTGCAGGTATCCACCGACATCGCCAGCTATTTAAGCTTCGTTATGGCGCTGTTTATGGCGTTTGGTGTCTCCTTTGAAGTGCCGGTAGCAATTGTGCTGCTGTGCTGGATGGGGATTACCTCGCCAGAAGACTTACGCAAAAAACGCCCGTATGTGCTGGTTGGTGCATTCGTTGTCGGGATGTTGCTGACGCCGCCGGATGTCTTCTCGCAAACGCTGTTG  >  NC_000913/4023097‑4023389

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