Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,426,592 T→A 18.3% intergenic (+35/+36) rplI → / ← yjfZ 50S ribosomal subunit protein L9/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,426,5920TA18.3% 424.6 / 82.9 175intergenic (+35/+36)rplI/yjfZ50S ribosomal subunit protein L9/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (71/72);  new base A (14/18);  total (85/90)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

AGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCT  >  NC_000913/4426498‑4426688
                                                                                              |                                                                                                
aGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTc                                                                                             <  1:705776/100‑1 (MQ=255)
    tCCAGGTTCACAGCGAAGTATTCGCGAAAGTTATTGTAAACGGAATAGCTGAGTAATTCGTTATTCAACGAGACGTAAAAACGCCGACTAATGGTCGGCg                                                                                         <  1:4063734/100‑1 (MQ=255)
    tCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCg                                                                                         <  1:3019079/100‑1 (MQ=255)
    tCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCg                                                                                         <  1:3785666/100‑1 (MQ=255)
    tCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCg                                                                                         <  2:2304620/100‑1 (MQ=255)
     ccAGGTTCACAGCGAAGTATTCGCGAAAGTGATCTTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCCATGCTCGGCGt                                                                                        <  2:2205756/100‑1 (MQ=255)
     ccAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTCAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGt                                                                                        <  1:2947404/100‑1 (MQ=255)
     ccAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGt                                                                                        <  1:2467011/100‑1 (MQ=255)
       aGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGttt                                                                                      <  1:1329164/100‑1 (MQ=255)
        ggTTCACATCGAAGTATTCGCGAAAGTGATCGTAGACGTAGTACCTGATAAATTCGATATGCAACGAGACGTAAAAACGCCGACCAATGGTCGGAGtttt                                                                                     <  2:2189943/100‑1 (MQ=255)
         gTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAATGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTg                                                                                    <  1:2957063/100‑1 (MQ=255)
          ttCACAGCGAAGTATTCGCGAAAGTGATCGTATACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGc                                                                                   <  2:2945390/100‑1 (MQ=255)
          ttCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtt                                                                                       >  2:1064039/1‑96 (MQ=255)
          ttCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGc                                                                                   >  1:3227094/1‑100 (MQ=255)
          ttCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGc                                                                                   <  1:112272/100‑1 (MQ=255)
          ttCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGttttta                                                                                   <  1:865100/100‑3 (MQ=255)
          ttCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATGCAACCAGACGCAAAAACGCCGACCAATGGTCGGCGttttta                                                                                   <  1:3149392/100‑3 (MQ=38)
          ttCACAGCGAAGTATTCGCAAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGCCGTAAAAACGCCGACCATTGGTCGGCGTTTTGc                                                                                   >  2:2084355/1‑100 (MQ=255)
           tCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCt                                                                                  >  2:3034759/1‑100 (MQ=255)
           tCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttttac                                                                                  <  2:1898826/100‑4 (MQ=255)
            cacaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCtt                                                                                 >  2:1028035/1‑100 (MQ=255)
            cacaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttttacg                                                                                 <  1:2063359/100‑5 (MQ=255)
             acaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCttt                                                                                >  1:3784513/1‑100 (MQ=255)
             acaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCttt                                                                                <  2:3229590/100‑1 (MQ=255)
             acaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCttt                                                                                >  1:1646734/1‑100 (MQ=255)
              caGCGAAGTATTCGAGAAAGTGATCGTAAAAGTAGTAGATGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGAACATTGGTCGGCGTATTGCtt                                                                                 >  2:495776/1‑98 (MQ=255)
                gCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTa                                                                             >  2:539452/1‑100 (MQ=255)
                 cGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAt                                                                            <  2:3289660/100‑1 (MQ=255)
                  gAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAtt                                                                           >  1:288775/1‑100 (MQ=255)
                  gAAGTATTCGCGAAAGTGATCGTAAACGTAGAAGCTGAATAATTCGTTATTCAACGAACCGTAACAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAtt                                                                           >  1:2824156/1‑100 (MQ=255)
                   aaGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAttt                                                                          <  1:1359505/100‑1 (MQ=255)
                   aaGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAttt                                                                          <  1:4161006/100‑1 (MQ=255)
                     gTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAttttt                                                                        <  1:2151372/100‑1 (MQ=255)
                      tATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTc                                                                       >  1:2896924/1‑100 (MQ=255)
                       agtCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCg                                                                      <  2:384144/98‑1 (MQ=255)
                       aTTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCg                                                                      >  2:4029394/1‑100 (MQ=255)
                       aTTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCg                                                                      >  2:918662/1‑100 (MQ=255)
                          cgcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCa                                                                   >  2:677136/1‑100 (MQ=255)
                           gcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTCTTTTTCGTCAg                                                                  <  2:1238495/100‑1 (MQ=255)
                            cgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAgg                                                                 >  1:2452606/1‑100 (MQ=255)
                             gAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAgg                                                                 <  2:3176151/99‑1 (MQ=255)
                                aGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAtt                                                             <  2:3996672/100‑1 (MQ=255)
                                  tGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGTTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAg                                                           <  2:1808478/100‑1 (MQ=255)
                                    aTCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGtt                                                         >  2:104470/1‑100 (MQ=255)
                                    aTCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGtt                                                         <  1:1028035/100‑1 (MQ=255)
                                    aTCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGtt                                                         >  2:1572431/1‑100 (MQ=255)
                                      cGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                       >  1:3204147/1‑100 (MQ=255)
                                      cGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                       <  2:4119756/100‑1 (MQ=255)
                                       gTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCg                                                      <  1:577756/100‑1 (MQ=255)
                                         aaaCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCa                                                    <  1:554370/100‑1 (MQ=255)
                                           aCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAg                                                  <  2:1661127/100‑1 (MQ=255)
                                              tagtagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAg                                               >  1:1739475/1‑100 (MQ=255)
                                               agtagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGa                                              <  1:2289316/100‑1 (MQ=255)
                                               agtagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGa                                              <  1:3300985/100‑1 (MQ=255)
                                                  agCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATcc                                           <  2:2511614/100‑1 (MQ=255)
                                                      gAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTTTTTTTTTTTTTTTTTCGCAAGTAGATCCAtt                                        >  2:1502225/1‑99 (MQ=255)
                                                        ataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAc                                     <  1:603339/100‑1 (MQ=255)
                                                         taatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAcc                                    <  2:3204133/100‑1 (MQ=255)
                                                         taatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAcc                                    >  1:2426874/1‑100 (MQ=255)
                                                         taatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAcc                                    <  2:3984786/100‑1 (MQ=255)
                                                          aatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCa                                   <  2:1480661/100‑1 (MQ=255)
                                                           atTCGTTATTCAACGAGACGTAAAAACGCCTACCATTGGTCGGAGTTTGGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGGAGATGCATTCAccc                                   >  1:1719917/1‑98 (MQ=255)
                                                            ttCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAtt                                 <  1:2051788/100‑1 (MQ=255)
                                                             tCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATtc                                >  1:1419335/1‑100 (MQ=255)
                                                               gTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATtctc                              >  2:2714650/1‑100 (MQ=255)
                                                                 tATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCtt                            >  1:3529490/1‑100 (MQ=255)
                                                                  aTTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtg                           <  1:4029421/100‑1 (MQ=255)
                                                                  aTTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtg                           <  1:895825/100‑1 (MQ=255)
                                                                  aTTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtg                           <  1:841371/100‑1 (MQ=255)
                                                                  aTTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtg                           <  2:604172/100‑1 (MQ=255)
                                                                   ttCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtgt                          >  1:1979961/1‑100 (MQ=255)
                                                                    tCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtgtg                         <  1:539453/100‑1 (MQ=255)
                                                                    tCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtgtg                         <  1:3985281/100‑1 (MQ=255)
                                                                      aaCGAGACGTAAAAACCCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTc                       <  2:2452601/100‑1 (MQ=255)
                                                                       aCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTcc                      <  2:2997574/100‑1 (MQ=255)
                                                                         gagaCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCAt                    <  1:4094662/100‑1 (MQ=255)
                                                                            aagTACAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTCTTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAg                 <  1:2718028/98‑1 (MQ=255)
                                                                             cGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGt                >  1:1198586/1‑100 (MQ=255)
                                                                             cGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGt                >  1:3598277/1‑100 (MQ=255)
                                                                              agcaaaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTc               >  1:2595685/4‑100 (MQ=255)
                                                                               tAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAATTCa              >  1:2774761/1‑100 (MQ=255)
                                                                               tAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCa              <  2:748178/100‑1 (MQ=255)
                                                                               tAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCa              >  1:2216835/1‑100 (MQ=255)
                                                                               tAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCa              >  1:1887089/1‑100 (MQ=255)
                                                                                caaaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAg             >  1:4184944/2‑100 (MQ=255)
                                                                                aaaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATGAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTCTCCATAAGTCAg             <  2:2484209/100‑1 (MQ=255)
                                                                                 aaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAgg            >  1:3811296/1‑100 (MQ=255)
                                                                                 aaaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAgg            >  1:3257542/1‑100 (MQ=255)
                                                                                   aaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGAt          <  1:2304492/100‑1 (MQ=255)
                                                                                   aaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGAt          >  2:1225190/1‑100 (MQ=255)
                                                                                   aaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCa              >  2:949751/1‑96 (MQ=255)
                                                                                    aCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATg         >  1:983392/1‑100 (MQ=255)
                                                                                     cGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATgg        >  2:1267063/1‑100 (MQ=255)
                                                                                      gCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGt       >  1:3900198/1‑100 (MQ=255)
                                                                                        caaCCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAATAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGttt     >  1:2838535/3‑100 (MQ=255)
                                                                                        cGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGttt     >  2:4152460/1‑100 (MQ=255)
                                                                                        cGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGttt     <  1:3431479/100‑1 (MQ=255)
                                                                                         gACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTg    >  2:3401933/1‑100 (MQ=255)
                                                                                         gACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTg    >  1:1680397/1‑100 (MQ=255)
                                                                                          aCCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGc   <  1:3680939/100‑1 (MQ=255)
                                                                                              |                                                                                                
AGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCT  >  NC_000913/4426498‑4426688

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