Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 962,373 A→T 52.9% I127F (ATT→TTT)  rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913962,3730AT52.9% 3.1 / 103.7 67I127F (ATT→TTT) rpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base A (16/15);  new base T (21/14);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-02

GTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATC  >  NC_000913/962233‑962516
                                                                                                                                            |                                                                                                                                               
gTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCg                                                                                                                                         >  4:265035/1‑149 (MQ=255)
 tGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGt                                                                                                                                        >  3:33981/1‑149 (MQ=255)
  gAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGtt                                                                                                                                       >  1:77537/1‑149 (MQ=255)
    aaCTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTcc                                                                                                                                     >  2:664018/1‑149 (MQ=255)
     aCTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCt                                                                                                                                    <  3:265035/149‑1 (MQ=255)
     aCTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCt                                                                                                                                    <  1:526603/149‑1 (MQ=255)
       tctGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGc                                                                                                                                  >  3:386652/1‑149 (MQ=255)
       tctGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGc                                                                                                                                  >  3:145608/1‑149 (MQ=255)
            tgTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGt                                                                                                                             >  1:599522/1‑149 (MQ=255)
                ccGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTtctc                                                                                                                         >  1:164094/1‑149 (MQ=255)
                 cGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTctct                                                                                                                        <  4:386652/149‑1 (MQ=255)
                 cGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTctct                                                                                                                        <  1:522993/149‑1 (MQ=255)
                        aaGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGAc                                                                                                                 >  2:540925/1‑149 (MQ=255)
                          gCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGt                                                                                                               >  2:45421/1‑149 (MQ=255)
                            tAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTc                                                                                                             >  2:603287/1‑149 (MQ=255)
                                 gTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCg                                                                                                        >  4:293238/1‑149 (MQ=255)
                                    aCGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTg                                                                                                     >  1:320109/1‑149 (MQ=255)
                                      gAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCg                                                                                                   >  4:526381/1‑149 (MQ=255)
                                         gCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGa                                                                                                >  2:657895/1‑149 (MQ=255)
                                            tGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGacac                                                                                             >  1:47396/1‑149 (MQ=255)
                                              gATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACActc                                                                                           <  4:672587/149‑1 (MQ=255)
                                              gATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTc                                                                                                                                              <  3:310917/98‑1 (MQ=255)
                                              gATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTc                                                                                                                                              >  4:310917/1‑98 (MQ=255)
                                              gATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACActc                                                                                           <  2:599522/149‑1 (MQ=255)
                                               aTCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACActct                                                                                          <  3:246323/149‑1 (MQ=255)
                                                      tGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTg                                                                                   <  2:646920/149‑1 (MQ=255)
                                                       ggAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTgg                                                                                  >  1:344633/1‑149 (MQ=255)
                                                        gAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGa                                                                                 >  3:103781/1‑149 (MQ=255)
                                                          aaaaGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAg                                                                               <  3:293238/149‑1 (MQ=255)
                                                               cTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCaaa                                                                          <  4:671100/149‑1 (MQ=255)
                                                                ttACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAag                                                                         <  4:145608/149‑1 (MQ=255)
                                                                    gAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCt                                                                     >  1:532758/1‑149 (MQ=255)
                                                                       gATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTCTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGa                                                                  >  4:241113/1‑149 (MQ=255)
                                                                         tGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAAt                                                                <  2:77065/149‑1 (MQ=255)
                                                                             gAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTaa                                                            >  3:468433/1‑149 (MQ=255)
                                                                             gAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTaa                                                            >  3:257416/1‑149 (MQ=255)
                                                                              aaaCTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTaaa                                                           <  1:45421/149‑1 (MQ=255)
                                                                                  tGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTaa                                                       >  4:444570/1‑149 (MQ=255)
                                                                                  tGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTaa                                                       >  1:281412/1‑149 (MQ=255)
                                                                                        cGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGc                                                 >  1:322500/1‑149 (MQ=255)
                                                                                         ggTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCt                                                >  2:251300/1‑149 (MQ=255)
                                                                                         ggTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCt                                                <  2:344633/149‑1 (MQ=255)
                                                                                         ggTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTCTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCt                                                <  3:241113/149‑1 (MQ=255)
                                                                                            gtTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGa                                             <  3:526381/149‑1 (MQ=255)
                                                                                                  aaCGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGaa                                       >  1:565783/1‑149 (MQ=255)
                                                                                                  aaCGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGaa                                       <  1:66002/149‑1 (MQ=255)
                                                                                                   aCGGCAAAGTTAAGGGCGGCTTGACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAATCTGGATCAGAAg                                      >  2:680573/1‑149 (MQ=255)
                                                                                                    cGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAgc                                     <  2:532758/149‑1 (MQ=255)
                                                                                                      gCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTCTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGGTCAGAAgcgc                                   >  4:645129/1‑149 (MQ=255)
                                                                                                      gCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAgcgc                                   <  1:228048/149‑1 (MQ=255)
                                                                                                            ttAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAg                                                          <  4:441571/120‑1 (MQ=255)
                                                                                                            ttAAGGGCGGCTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAg                                                          >  3:441571/1‑120 (MQ=255)
                                                                                                                    ggcTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACgttgttgt                     <  4:103781/149‑1 (MQ=255)
                                                                                                                      cTTCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGttt                   <  2:565783/149‑1 (MQ=255)
                                                                                                                       ttCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTtc                  >  4:666851/1‑149 (MQ=255)
                                                                                                                        tCACTGTTGAGCTGAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTtct                 >  1:575115/1‑149 (MQ=255)
                                                                                                                         cACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTtctc                <  2:320109/149‑1 (MQ=255)
                                                                                                                          aCTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCtt                                                                    <  3:332863/96‑1 (MQ=255)
                                                                                                                          aCTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCtt                                                                    >  4:332863/1‑96 (MQ=255)
                                                                                                                                   cTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTaaa                                                           >  2:68398/1‑96 (MQ=255)
                                                                                                                                   cTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTaaa                                                           <  1:68398/96‑1 (MQ=255)
                                                                                                                                    tGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGtt     <  3:666851/149‑1 (MQ=255)
                                                                                                                                     gAACGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTa    >  1:517113/1‑149 (MQ=255)
                                                                                                                                      aaCGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTAt   >  2:377656/1‑149 (MQ=255)
                                                                                                                                      aaCGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGaa                                                                                <  4:637081/72‑1 (MQ=255)
                                                                                                                                      aaCGGTATGCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGaa                                                                                >  3:637081/1‑72 (MQ=255)
                                                                                                                                       aCGGTTTTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATc  <  2:575115/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATC  >  NC_000913/962233‑962516

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