Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,309,164 A→G R54R (CGT→CGC)  rpsO ← 30S ribosomal subunit protein S15

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,309,1640AG98.7% 492.3 / ‑3.3 155R54R (CGT→CGC) rpsO30S ribosomal subunit protein S15
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (83/70);  total (83/72)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 169 total aligned reads displayed.

CGCTCGATGAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCACGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCTGAACTTCGG  >  NC_000913/3309075‑3309252
                                                                                         |                                                                                        
cgcTCGATGAGCTCGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCgcga                                                                                       <  1:1502229/93‑1 (MQ=255)
 gcTCGATGAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGACACCATGCGCAGCAGACCgcgac                                                                                      >  2:2683527/1‑93 (MQ=255)
 gcTCGATGAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCgcgac                                                                                      <  1:2882090/93‑1 (MQ=255)
 gcTCGATGAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCgcgac                                                                                      <  2:4642854/93‑1 (MQ=255)
 gcTCGATGAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCgcgac                                                                                      <  2:2553035/93‑1 (MQ=255)
       tGAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCtgt                                                                                <  1:281659/93‑1 (MQ=255)
       tGAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCtgt                                                                                <  1:3871210/93‑1 (MQ=255)
        gAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTGGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCtggg                                                                               >  1:208954/1‑91 (MQ=255)
        gAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCtgtg                                                                               <  1:3093393/93‑1 (MQ=255)
        gAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCtgtg                                                                               <  1:563482/93‑1 (MQ=255)
        gAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCACGACGGCtgtg                                                                               <  1:3304537/93‑1 (MQ=255)
         aGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCATACCGCGACGGCTgtgg                                                                              <  1:2684946/93‑1 (MQ=255)
          gCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTgtggt                                                                             >  1:4723483/1‑93 (MQ=255)
           cTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTgtggtg                                                                            <  2:4163444/93‑1 (MQ=255)
            tGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGa                                                                           >  1:1038576/1‑93 (MQ=255)
            tGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGa                                                                           >  2:260524/1‑93 (MQ=255)
            tGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGa                                                                           >  2:2057006/1‑93 (MQ=255)
             gggTGTGACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGAt                                                                          >  2:1372272/1‑93 (MQ=255)
             gggTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGAt                                                                          >  1:3336469/1‑93 (MQ=255)
              ggTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATc                                                                         >  1:2166704/1‑93 (MQ=255)
              ggTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATc                                                                         <  2:1994941/93‑1 (MQ=255)
              ggTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATc                                                                         >  1:3072094/1‑93 (MQ=255)
                 gtAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCttt                                                                      <  1:1840607/93‑1 (MQ=255)
                 gtAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCttt                                                                      >  2:3700909/1‑93 (MQ=255)
                  tAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGAACGAGACGGCTGTGGCGATGtttt                                                                     >  1:3359605/1‑93 (MQ=255)
                   aaCGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCttttt                                                                    >  2:1339664/1‑93 (MQ=255)
                    aCGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCtttttt                                                                   >  1:738959/1‑93 (MQ=255)
                      gTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTtgt                                                                 >  2:148628/1‑93 (MQ=255)
                      gTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTtgt                                                                 <  1:2466207/93‑1 (MQ=255)
                          tACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGctc                                                             <  2:633111/93‑1 (MQ=255)
                          tACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGctc                                                             >  2:433802/1‑93 (MQ=255)
                          tACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGctc                                                             <  2:3111124/93‑1 (MQ=255)
                          tACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGctc                                                             <  2:1239419/93‑1 (MQ=255)
                          tACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGctc                                                             >  1:2861679/1‑93 (MQ=255)
                           aCGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGctct                                                            >  2:37188/1‑93 (MQ=255)
                           aCGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGctct                                                            >  2:1961712/1‑93 (MQ=255)
                           aCGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGctct                                                            <  2:626387/93‑1 (MQ=255)
                              tCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCa                                                         >  2:3431937/1‑93 (MQ=255)
                                tttACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCaaa                                                       >  1:348262/1‑93 (MQ=255)
                                tttACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCaaa                                                       >  1:169017/1‑93 (MQ=255)
                                tttACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCaaa                                                       <  1:3273853/93‑1 (MQ=255)
                                tttACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCaaa                                                       >  1:122785/1‑93 (MQ=255)
                                  tACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGt                                                     <  1:1682840/93‑1 (MQ=255)
                                   aCGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTg                                                    <  2:2910764/93‑1 (MQ=255)
                                   aCGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCACGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTg                                                    <  1:1965204/93‑1 (MQ=255)
                                        tCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCt                                               >  1:4110771/1‑93 (MQ=255)
                                        tCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCt                                               <  1:4325948/93‑1 (MQ=255)
                                        tCAGGTAATCGAGAAGTTTACGACGCTGAGAAACCATGAGCAGCGGCCCCAGCCGGCTGTGGTGATCTTTATTGTGCTATGCAAAGTGGCCCt                                               <  2:2149345/93‑1 (MQ=255)
                                          aGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGc                                             >  2:4368187/1‑93 (MQ=255)
                                          aGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGc                                             >  1:1512219/1‑93 (MQ=255)
                                           ggTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCa                                            <  1:2045765/93‑1 (MQ=255)
                                            gTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAg                                           <  1:4342325/93‑1 (MQ=255)
                                             tAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAgg                                          >  1:3738783/1‑93 (MQ=255)
                                              aGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAggt                                         <  2:3096915/93‑1 (MQ=255)
                                               gTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCGAAGTGGCCCTGCAggtg                                        <  1:473361/93‑1 (MQ=255)
                                               gTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAggtg                                        >  1:3813427/1‑93 (MQ=255)
                                                tCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAggtgg                                       >  1:418657/1‑93 (MQ=255)
                                                 cGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAggtggt                                      <  2:2499143/93‑1 (MQ=255)
                                                  gAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGtt                                     >  1:1403630/1‑93 (MQ=255)
                                                   aGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTg                                    >  2:4491417/1‑93 (MQ=255)
                                                    gCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGa                                   >  2:537749/1‑93 (MQ=255)
                                                      aGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATc                                 >  2:3513113/1‑93 (MQ=255)
                                                      aGTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATc                                 <  1:2611723/93‑1 (MQ=255)
                                                       gTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGctct                                >  1:2836669/1‑93 (MQ=255)
                                                       gTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCt                                <  2:1405865/93‑1 (MQ=255)
                                                       gTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCt                                >  1:861338/1‑93 (MQ=255)
                                                       gTTTACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCt                                <  1:2887256/93‑1 (MQ=255)
                                                       gTTTACGACGCTAAGAAACCATGCGAAGCAGACCGCGTCGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCt                                <  2:2244743/93‑1 (MQ=255)
                                                         ttACGACGCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCtgt                              <  2:2194029/93‑1 (MQ=255)
                                                           acgacgCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGc                            <  2:47858/93‑1 (MQ=255)
                                                              acgCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGt                         <  2:3977392/93‑1 (MQ=255)
                                                               cgCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTc                        <  1:2508805/93‑1 (MQ=255)
                                                                gCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTca                       <  1:581881/93‑1 (MQ=255)
                                                                gCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTca                       >  2:2779657/1‑93 (MQ=255)
                                                                gCTGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTca                       <  1:723824/93‑1 (MQ=255)
                                                                  tGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTcagc                     >  1:2766251/1‑93 (MQ=255)
                                                                  tGAGAAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTcagc                     >  2:2021630/1‑93 (MQ=255)
                                                                   gagaAACCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTcagca                    >  1:1650530/1‑93 (MQ=255)
                                                                      aaaCCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTg                 <  2:449900/93‑1 (MQ=255)
                                                                      aaaCCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTg                 >  1:2726419/1‑93 (MQ=255)
                                                                      aaaCCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTg                 >  2:4575681/1‑93 (MQ=255)
                                                                       aaCCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGc                >  1:4475212/1‑93 (MQ=255)
                                                                        aCCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCt               >  2:2607145/1‑93 (MQ=255)
                                                                        aCCATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCt               >  1:4649450/1‑93 (MQ=255)
                                                                         ccATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTa              <  2:3372816/93‑1 (MQ=255)
                                                                         ccATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTa              <  1:3550296/93‑1 (MQ=255)
                                                                         ccATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTa              <  1:3550297/93‑1 (MQ=255)
                                                                         ccATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTa              >  1:1025011/1‑93 (MQ=255)
                                                                          cATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTAc             <  2:4287219/93‑1 (MQ=255)
                                                                          cATGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTAc             <  2:843519/93‑1 (MQ=255)
                                                                           aTGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTAcc            >  1:1161150/1‑93 (MQ=255)
                                                                           aTGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTAcc            <  1:902624/93‑1 (MQ=255)
                                                                            tGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCt           >  1:499199/1‑93 (MQ=255)
                                                                            tGCGCAGCAGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCt           >  1:1246357/1‑93 (MQ=255)
                                                                               gcagcaGACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCTGaa        >  1:4570094/1‑93 (MQ=255)
                                                                                cagcagACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTAAGCAGTGCTACCTGAAc       >  2:3629318/1‑93 (MQ=255)
                                                                                 agcagACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCTGAACt      <  1:3174832/93‑1 (MQ=255)
                                                                                  gcagACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCTGAACtt     >  2:138119/1‑93 (MQ=255)
                                                                                  gcagACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCTGAACtt     <  2:956256/93‑1 (MQ=255)
                                                                                     gACCGCGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCTGAACTTCgg  <  1:4288299/93‑1 (MQ=255)
                                                                                         |                                                                                        
CGCTCGATGAGCTGGGTGTAACGTGCTACGTCTTTACGTTTCAGGTAGTCGAGCAGTTTACGACGCTGAGAAACCATGCGCAGCAGACCACGACGGCTGTGGTGATCTTTTTTGTGCTCTGCAAAGTGGCCCTGCAGGTGGTTGATCTGTGCAGTCAGCAGTGCTACCTGAACTTCGG  >  NC_000913/3309075‑3309252

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