Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 961,562 A→G K115K (AAA→AAG rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913961,5620AG94.3% 497.6 / 14.5 177K115K (AAA→AAGrpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base A (8/2);  new base G (106/61);  total (114/63)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 180 total aligned reads displayed.

CCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCT  >  NC_000913/961472‑961651
                                                                                          |                                                                                         
ccGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGt                                                                                         >  2:3983256/1‑93 (MQ=255)
ccGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGt                                                                                         >  2:2083109/1‑93 (MQ=255)
ccGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGt                                                                                         <  1:1539323/93‑1 (MQ=255)
 cGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGtt                                                                                        >  1:921546/1‑93 (MQ=255)
  gTGAGAAAGCTAAGCGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACAGGTGTTATCAACGGCAAGGTTa                                                                                       <  2:2994004/93‑1 (MQ=255)
      gaAAGCTACACGTCACGAAGCCTGGATCACGCTAGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAggg                                                                                   <  1:2204255/93‑1 (MQ=255)
      gaAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAggg                                                                                   >  1:3430272/1‑93 (MQ=255)
      gaAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAggg                                                                                   >  1:904918/1‑93 (MQ=255)
      gaAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAggg                                                                                   <  1:731901/93‑1 (MQ=255)
       aaaGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGggc                                                                                  <  1:151327/93‑1 (MQ=255)
        aaGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGggcg                                                                                 >  1:4011929/1‑93 (MQ=255)
        aaGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGggcg                                                                                 <  1:1886861/93‑1 (MQ=255)
         aGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGggcgg                                                                                >  1:406728/1‑93 (MQ=255)
         aGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGggcgg                                                                                >  2:2067230/1‑93 (MQ=255)
         aGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGCTATCAACGGCAAGGTTAAGggcgg                                                                                >  1:4681606/1‑93 (MQ=255)
          gCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGggcggc                                                                               <  1:3234755/93‑1 (MQ=255)
          gCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGggcggc                                                                               >  1:1288920/1‑93 (MQ=255)
           cTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCt                                                                              >  1:4309888/1‑93 (MQ=255)
                cGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACt                                                                         <  1:1001074/93‑1 (MQ=255)
                 gTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTCCCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTAACtt                                                                        >  2:3704630/1‑92 (MQ=255)
                 gTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTg                                                                        >  1:1046444/1‑93 (MQ=255)
                    aCGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTg                                                                     >  2:50941/1‑93 (MQ=255)
                    aCGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTg                                                                     >  2:2112975/1‑93 (MQ=255)
                    aCGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTg                                                                     >  2:3074349/1‑93 (MQ=255)
                     cGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGa                                                                    <  1:2540875/93‑1 (MQ=255)
                      gAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAg                                                                   >  2:4218024/1‑93 (MQ=255)
                       aaGCCTGTATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGc                                                                  <  2:2752976/93‑1 (MQ=255)
                       aaGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGc                                                                  <  1:2093759/93‑1 (MQ=255)
                       aaGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGc                                                                  >  1:2236286/1‑93 (MQ=255)
                         gCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTg                                                                >  1:263040/1‑93 (MQ=255)
                          ccTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGa                                                               >  1:1514585/1‑93 (MQ=255)
                          ccTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGa                                                               >  2:4532182/1‑93 (MQ=255)
                           cTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGaa                                                              <  2:1899483/93‑1 (MQ=255)
                           cTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGaa                                                              >  2:1118976/1‑93 (MQ=255)
                           cTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGaa                                                              <  1:1525291/93‑1 (MQ=255)
                           cTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGAGGCTTCACTGTTGAGCCGaa                                                              >  2:4441193/1‑93 (MQ=255)
                            tGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAAc                                                             >  2:4082263/1‑93 (MQ=255)
                             ggATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACg                                                            >  2:383214/1‑93 (MQ=255)
                             ggATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACg                                                            <  1:3318850/93‑1 (MQ=255)
                             ggATCACGCTGGAAAAAGCCTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACg                                                            >  2:4579366/1‑93 (MQ=255)
                                tCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTa                                                         >  2:1840417/1‑93 (MQ=255)
                                tCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTa                                                         >  1:2464125/1‑93 (MQ=255)
                                 cACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTAt                                                        >  1:653981/1‑93 (MQ=255)
                                 cACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTAt                                                        <  2:3767402/93‑1 (MQ=255)
                                  aCGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTAtt                                                       <  1:3553134/93‑1 (MQ=255)
                                  aCGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTAtt                                                       <  2:323785/93‑1 (MQ=255)
                                  aCGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGAAtt                                                       >  1:752037/1‑93 (MQ=255)
                                  aCGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTGAAGGGCGGCTTCACTGTTGAGTTGAACGGTAtt                                                       >  2:172716/1‑93 (MQ=255)
                                  aCGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTAtt                                                       >  1:1457589/1‑93 (MQ=255)
                                    gCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCg                                                     >  2:3131723/1‑93 (MQ=255)
                                     cTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGt                                                    <  2:2582409/93‑1 (MQ=255)
                                     cTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGt                                                    <  2:3912683/93‑1 (MQ=255)
                                      tGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTg                                                   >  1:1820661/1‑93 (MQ=255)
                                        gAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCg                                                 >  1:3101339/1‑93 (MQ=255)
                                        gAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCg                                                 >  1:3720714/1‑93 (MQ=255)
                                        gAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCg                                                 <  1:4145959/93‑1 (MQ=255)
                                         aaaaaGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGt                                                >  1:1841873/1‑93 (MQ=255)
                                          aaaaGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGtt                                               >  1:698849/1‑93 (MQ=255)
                                          aaaaGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGtt                                               <  1:3006862/93‑1 (MQ=255)
                                          aaaaGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGGTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGtt                                               >  2:2831628/1‑93 (MQ=255)
                                           aaaGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTc                                              >  2:3965153/1‑93 (MQ=255)
                                            aaGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTcc                                             <  1:2727475/93‑1 (MQ=255)
                                             aGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCt                                            >  2:2901325/1‑93 (MQ=255)
                                              gCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTg                                           >  1:720128/1‑93 (MQ=255)
                                              gCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTg                                           >  1:4466706/1‑93 (MQ=255)
                                                  aCGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAg                                       >  2:3354618/1‑93 (MQ=255)
                                                   cGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAgg                                      >  2:3634417/1‑93 (MQ=255)
                                                    gAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGt                                     >  2:1132148/1‑93 (MQ=255)
                                                      aGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTtc                                   <  2:1940771/93‑1 (MQ=255)
                                                      aGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTtc                                   >  1:2698778/1‑93 (MQ=255)
                                                         tGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTctct                                >  2:2036393/1‑93 (MQ=255)
                                                         tGCTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTctct                                <  1:3963962/93‑1 (MQ=255)
                                                           cTGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTgg                              <  1:3865709/93‑1 (MQ=255)
                                                            tGAAACTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGt                             >  1:2887573/1‑93 (MQ=255)
                                                               aaCTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGa                          >  2:3539019/1‑93 (MQ=255)
                                                                aCTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGAc                         <  1:1492049/93‑1 (MQ=255)
                                                                aCTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGAc                         <  1:1562884/93‑1 (MQ=255)
                                                                aCTGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGAc                         >  1:2303460/1‑93 (MQ=255)
                                                                  tGTTACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGt                       >  1:4576604/1‑93 (MQ=255)
                                                                    tttccGGGGGTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCGTGCCAGGTTCTCTGGTAGACGTTc                     <  2:4518083/90‑1 (MQ=255)
                                                                    ttACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTc                     >  2:1242007/1‑93 (MQ=255)
                                                                     tACCGGTGTTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCg                    >  2:2007585/1‑93 (MQ=255)
                                                                          gtgtTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCgg               >  2:2471641/1‑93 (MQ=255)
                                                                          gtgtTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCgg               >  2:2471640/1‑93 (MQ=255)
                                                                           tgtTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCggg              >  1:347526/1‑92 (MQ=255)
                                                                           tgtTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGt              <  1:2737969/93‑1 (MQ=255)
                                                                            gtTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCgggg             >  2:3541794/1‑91 (MQ=255)
                                                                            gtTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTg             >  2:1750243/1‑93 (MQ=255)
                                                                            gtTATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTg             <  1:4192127/93‑1 (MQ=255)
                                                                              tATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCg           >  2:1542505/1‑93 (MQ=255)
                                                                              tATCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCg           >  2:1129247/1‑93 (MQ=255)
                                                                                tCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTg         >  2:3589044/1‑93 (MQ=255)
                                                                                tCAACGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTg         >  2:4107037/1‑93 (MQ=255)
                                                                                  aaCGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGac       <  1:4300937/93‑1 (MQ=255)
                                                                                    cGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGacac     >  1:2827767/1‑93 (MQ=255)
                                                                                    cGGCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGacac     <  2:594755/93‑1 (MQ=255)
                                                                                     ggCAAGGTTAAGGGCGGCTTCACTGTTGGGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACAct    >  2:838960/1‑93 (MQ=255)
                                                                                      gCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTTCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACActc   <  1:536128/93‑1 (MQ=255)
                                                                                      gCAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACActc   <  1:3533523/93‑1 (MQ=255)
                                                                                       cAAGGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACActct  >  1:4298887/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAAGCTTACGAAGATGCTGAAACTGTTACCGGTGTTATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCT  >  NC_000913/961472‑961651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: