Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,426,592 T→A 10.5% 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,5920TA10.5% 432.6 / 31.0 134intergenic (+35/+36)rplI/yjfZ50S ribosomal subunit protein L9/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (58/61);  new base A (5/9);  total (64/70)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 141 total aligned reads displayed.

AGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCT  >  NC_000913/4426444‑4426740
                                                                                                                                                    |                                                                                                                                                    
agCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAt                                                                                                                                                      >  3:295475/1‑149 (MQ=255)
 gCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAtt                                                                                                                                                     >  4:523206/1‑149 (MQ=255)
      gTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTc                                                                                                                                                 >  1:769549/1‑148 (MQ=255)
      gTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTc                                                                                                                                                 <  2:769549/148‑1 (MQ=255)
      gTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCg                                                                                                                                                >  2:343600/1‑149 (MQ=255)
        tCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGc                                                                                                                                              <  3:166043/149‑1 (MQ=255)
        tCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGc                                                                                                                                              <  3:5246/149‑1 (MQ=255)
         cGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCg                                                                                                                                             <  3:917961/149‑1 (MQ=255)
           tCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtt                                                                                                                                           >  3:328103/1‑149 (MQ=255)
           tCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtt                                                                                                                                           <  3:766596/149‑1 (MQ=255)
           tCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtt                                                                                                                                           <  2:905691/149‑1 (MQ=255)
             tGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtttt                                                                                                                                         <  2:760378/149‑1 (MQ=255)
              gCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTACAGGTTCACTGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGttttt                                                                                                                                        <  3:542928/149‑2 (MQ=255)
               ccGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGc                                                                                                                                       <  3:279986/149‑1 (MQ=255)
                 gAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCtt                                                                                                                                     <  2:740437/149‑1 (MQ=255)
                           cTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTc                                                                                                                           <  2:581291/149‑1 (MQ=255)
                               gTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCa                                                                                                                       <  1:885721/149‑1 (MQ=255)
                                tACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAg                                                                                                                      <  4:609002/149‑1 (MQ=255)
                                tACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAg                                                                                                                      <  3:523206/149‑1 (MQ=255)
                                   cACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTa                                                                                                                   >  2:363342/1‑149 (MQ=255)
                                         cGAACACGAAGTGAGCTTCCAGGTTCACAGCCAAGTATTCGCGAAAGTGATCGTACACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGTGTTTTGCTTTCTATTTTTCGTCAGGTAGTAGtt                                                                                                             <  1:315764/149‑1 (MQ=255)
                                          gAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTg                                                                                                                                                    <  3:32183/109‑1 (MQ=255)
                                          gAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTg                                                                                                                                                    >  4:32183/1‑109 (MQ=255)
                                          gAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGttt                                                                                                            >  1:566856/1‑149 (MQ=255)
                                          gAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGttt                                                                                                            <  4:557005/149‑1 (MQ=255)
                                          gAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGttt                                                                                                            <  1:849972/149‑1 (MQ=255)
                                              acGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGACGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCa                                                                                                        >  1:852683/1‑149 (MQ=255)
                                               cGAAGTGAGCTTCCAGGTTCACAGCGTAGTATTCGCGAAAGTGATCATAAACGTAGTAGCGGAAAAATTCCTTATTCAACGAGACGTAACATCGCCGAGCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCACGTATTAGTTTCGCaa                                                                                                       <  2:494584/149‑1 (MQ=255)
                                               cGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCaa                                                                                                       <  4:153531/149‑1 (MQ=255)
                                                gatgTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCACTGGTCGGCGTTTTGCTTTCTAttt                                                                                                                              >  4:355233/4‑125 (MQ=255)
                                                gAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAttt                                                                                                                              <  3:355233/125‑1 (MQ=255)
                                                   gTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGACTAATTCGTTATTCAACGAGACGTAAAAACGCCGAGCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAg                                                                                                   >  2:294210/1‑149 (MQ=255)
                                                    tGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGa                                                                                                  >  1:761294/1‑149 (MQ=255)
                                                             aGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAc                                                                                         >  2:750199/1‑149 (MQ=255)
                                                                ttCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAt                                                                                      >  2:676939/1‑149 (MQ=255)
                                                                     aGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTctct                                                                                 >  3:752763/1‑149 (MQ=255)
                                                                     aGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTctct                                                                                 <  4:746139/149‑1 (MQ=255)
                                                                        gAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtgt                                                                              <  1:54344/149‑1 (MQ=255)
                                                                            tATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTcc                                                                          >  1:672123/1‑149 (MQ=255)
                                                                            tATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTcc                                                                          <  3:832057/149‑1 (MQ=255)
                                                                             aTTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCa                                                                         <  1:1071701/149‑1 (MQ=255)
                                                                              ttCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCAt                                                                        >  2:316412/1‑149 (MQ=255)
                                                                              ttCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCAt                                                                        >  1:66997/1‑149 (MQ=255)
                                                                                cgcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCTACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATaa                                                                      <  4:298596/149‑1 (MQ=255)
                                                                                cgcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATaa                                                                      <  4:959439/149‑1 (MQ=255)
                                                                                     aaGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAg                                                                 <  4:199709/149‑1 (MQ=255)
                                                                                     aaGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAg                                                                 <  4:408693/149‑1 (MQ=255)
                                                                                       gTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa                                                               <  3:140450/149‑1 (MQ=255)
                                                                                        tGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGAt                                                              <  2:649376/149‑1 (MQ=255)
                                                                                             gTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa                                                               <  2:446341/143‑1 (MQ=255)
                                                                                             gTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGttt                                                         <  2:663307/149‑1 (MQ=255)
                                                                                             gTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGttt                                                         >  4:358953/1‑149 (MQ=255)
                                                                                              tAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTg                                                        >  1:655947/1‑149 (MQ=255)
                                                                                              tAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTg                                                        >  3:680655/1‑149 (MQ=255)
                                                                                                aaCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCt                                                      >  4:23147/1‑149 (MQ=255)
                                                                                                 aCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCCTTCTCtt                                                                                <  2:964109/122‑1 (MQ=255)
                                                                                                 aCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCtt                                                                                >  1:964109/1‑122 (MQ=255)
                                                                                                   gtagtaGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCAt                                                                        <  3:21732/128‑1 (MQ=255)
                                                                                                   gtagtaGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCAt                                                                        >  4:21732/1‑128 (MQ=255)
                                                                                                   gtagtaGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGc                                                   <  2:910776/149‑1 (MQ=255)
                                                                                                     agtagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGt                                                 >  3:207427/1‑149 (MQ=255)
                                                                                                         gCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGttt                                                         <  2:543454/137‑1 (MQ=255)
                                                                                                             aataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGataca                                         >  2:180322/1‑149 (MQ=255)
                                                                                                                aatTCATTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGatacatac                                      <  1:343600/149‑1 (MQ=255)
                                                                                                                  ttCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCTGGATGGTTTGCTGGCGTTAGATACATACGa                                    >  1:715265/1‑149 (MQ=255)
                                                                                                                  ttCGTTATTCAACGAGACGTAAAAACGCCGACCATCGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGa                                    <  4:704585/149‑1 (MQ=255)
                                                                                                                      ttATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtttt                                                                                                                                         >  4:189133/1‑44 (MQ=255)
                                                                                                                      ttATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtttt                                                                                                                                         <  3:189133/44‑1 (MQ=255)
                                                                                                                       tATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCa                               <  3:352720/149‑1 (MQ=255)
                                                                                                                        aTTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCaa                              >  3:378989/1‑149 (MQ=255)
                                                                                                                          tCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACt                            >  1:411913/1‑149 (MQ=255)
                                                                                                                            aaCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACt                            >  2:1067112/1‑147 (MQ=255)
                                                                                                                            aaCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACt                            <  1:1067112/147‑1 (MQ=255)
                                                                                                                             aCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTTAGGTATTAGTTTCGCAAGTAGATCCATTCACCGTTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGAt                         <  1:1069568/149‑1 (MQ=255)
                                                                                                                             aCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGAt                         <  2:184745/149‑1 (MQ=255)
                                                                                                                              cGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATg                        >  3:207537/1‑149 (MQ=255)
                                                                                                                                 gaCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGtct                     >  4:580969/1‑149 (MQ=255)
                                                                                                                                  aCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGtctc                    <  1:395440/149‑1 (MQ=255)
                                                                                                                                  aCGTAAAAACCCCGACCATTGGTCGGTGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGACCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACTTACGAGCGCAACTGATGtctc                    <  1:791025/149‑1 (MQ=255)
                                                                                                                                   cGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCg                   <  2:272212/149‑1 (MQ=255)
                                                                                                                                    gTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGc                  <  3:713136/149‑1 (MQ=255)
                                                                                                                                    gTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGc                  <  3:580969/149‑1 (MQ=255)
                                                                                                                                    agcaaaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGc                  >  2:512063/4‑149 (MQ=255)
                                                                                                                                     tAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCt                 >  2:784488/1‑149 (MQ=255)
                                                                                                                                      aaaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTc                >  1:502323/1‑149 (MQ=255)
                                                                                                                                        aaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCt              >  2:686473/1‑149 (MQ=255)
                                                                                                                                        aaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTCGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGACCGCAACTGATGTCTCGCTCCt              >  1:131311/1‑149 (MQ=255)
                                                                                                                                        aaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCt              >  4:139873/1‑149 (MQ=255)
                                                                                                                                         aaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTg             >  3:581662/1‑149 (MQ=255)
                                                                                                                                         aaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                                                                           >  4:842440/1‑55 (MQ=255)
                                                                                                                                         aaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                                                                           <  3:842440/55‑1 (MQ=255)
                                                                                                                                           cGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCg           >  4:988363/1‑149 (MQ=255)
                                                                                                                                            gCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCgg          <  1:512063/149‑1 (MQ=255)
                                                                                                                                               gACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCAc       >  2:165033/1‑149 (MQ=255)
                                                                                                                                                  cATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcg    >  3:390972/1‑149 (MQ=255)
                                                                                                                                                   aTTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcgc   >  1:282579/1‑149 (MQ=255)
                                                                                                                                                   aTTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcgc   >  3:693182/1‑149 (MQ=255)
                                                                                                                                                    ttGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa                                                               <  1:825844/88‑1 (MQ=255)
                                                                                                                                                    ttGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcgct  <  4:214762/149‑1 (MQ=255)
                                                                                                                                                    ttGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGCCAGGa                                                               <  1:314249/88‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCT  >  NC_000913/4426444‑4426740

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: