Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,567 G→A C254Y (TGC→TAC)  yaaU → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,5670GA94.8% 583.7 / 10.0 173C254Y (TGC→TAC) yaaUputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (7/2);  new base A (96/68);  total (103/70)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 176 total aligned reads displayed.

TGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGG  >  NC_000913/46478‑46657
                                                                                         |                                                                                          
tGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCa                                                                                         >  2:4155662/1‑93 (MQ=255)
 gAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAg                                                                                        <  2:2887397/93‑1 (MQ=255)
   agagCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGt                                                                                      <  2:2144298/93‑1 (MQ=255)
    gagCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTg                                                                                     >  1:3836555/1‑93 (MQ=255)
      gcagcTGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGAt                                                                                   <  2:4713844/93‑1 (MQ=255)
         gcCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATccc                                                                                <  1:4388338/93‑1 (MQ=255)
            gcagcaAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAAt                                                                             <  2:684644/93‑1 (MQ=255)
            gcagcaAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAAt                                                                             <  2:285026/93‑1 (MQ=255)
             cagcaAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATg                                                                            >  1:798968/1‑93 (MQ=255)
              agcaAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGt                                                                           >  2:1970501/1‑93 (MQ=255)
               gcaAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGtt                                                                          <  1:1042262/93‑1 (MQ=255)
                caAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTc                                                                         >  2:4439905/1‑93 (MQ=255)
                caAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTc                                                                         >  2:3722963/1‑93 (MQ=255)
                caAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTc                                                                         <  1:1904856/93‑1 (MQ=255)
                caAACCCGTTTTCGCGATCTGTTTAATCGCCCCCATTTTCCTTTTGTTCTGTTTGTTTCCGCCATCTGGACCTACCAGGTGATCCCAATCTTc                                                                         >  2:90626/1‑93 (MQ=255)
                 aaaCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCg                                                                        >  2:4469921/1‑93 (MQ=255)
                 aaaCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCg                                                                        >  2:4187446/1‑93 (MQ=255)
                      cGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCAtt                                                                   <  1:3164778/93‑1 (MQ=255)
                      cGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCAtt                                                                   <  2:3557729/93‑1 (MQ=255)
                      cGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCAtt                                                                   >  1:2595169/1‑93 (MQ=255)
                       gTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCAttt                                                                  >  1:5163340/1‑93 (MQ=255)
                       gTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCAttt                                                                  >  1:1460371/1‑93 (MQ=255)
                       gTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCAttt                                                                  >  1:4953213/1‑93 (MQ=255)
                        ttttCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTa                                                                 >  1:5680906/1‑93 (MQ=255)
                        ttttCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTa                                                                 <  1:4337620/93‑1 (MQ=255)
                        ttttCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTa                                                                 >  2:2361857/1‑93 (MQ=255)
                          ttCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTaca                                                               <  2:579415/93‑1 (MQ=255)
                            cgcgATCTGTTTAATCGCCGCTATTTTTCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACAcc                                                             <  2:204647/93‑1 (MQ=255)
                            cgcgATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACAcc                                                             <  1:941136/93‑1 (MQ=255)
                             gcgATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCt                                                            >  2:5604174/1‑93 (MQ=255)
                             gcgATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCt                                                            >  1:637960/1‑93 (MQ=255)
                              cgATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCtt                                                           <  1:4305046/93‑1 (MQ=255)
                              cgATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCtt                                                           <  2:6037281/93‑1 (MQ=255)
                               gATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCttt                                                          <  2:2333033/93‑1 (MQ=255)
                               gATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCttt                                                          >  2:6028995/1‑93 (MQ=255)
                                aTCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTg                                                         >  2:5353269/1‑93 (MQ=255)
                                aTCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTg                                                         >  1:1161400/1‑93 (MQ=255)
                                aTCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTg                                                         <  1:2930380/93‑1 (MQ=255)
                                aTCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTAAACCttgg                                                         >  1:667545/1‑91 (MQ=255)
                                 tCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTgg                                                        >  1:5282017/1‑93 (MQ=255)
                                   tGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGcc                                                      >  2:618307/1‑93 (MQ=255)
                                    gTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCTCCATTTATACCTTTGGccc                                                     >  1:3688167/1‑93 (MQ=255)
                                    gTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGccc                                                     >  1:4181732/1‑93 (MQ=255)
                                      ttAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGc                                                   >  2:2616947/1‑93 (MQ=255)
                                      ttAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGc                                                   >  1:476921/1‑93 (MQ=255)
                                         aTCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCaaa                                                <  2:2509445/93‑1 (MQ=255)
                                         aTCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCaaa                                                >  2:1940268/1‑93 (MQ=255)
                                          tCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAAt                                               <  2:2816238/93‑1 (MQ=255)
                                          tCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAAt                                               <  1:917649/93‑1 (MQ=255)
                                           cgccgcCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATc                                              <  1:2268844/93‑1 (MQ=255)
                                            gccgccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCg                                             >  1:1742764/1‑93 (MQ=255)
                                             ccgccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGt                                            <  2:321990/93‑1 (MQ=255)
                                               gccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTg                                          >  1:5513644/1‑93 (MQ=255)
                                               gccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTg                                          <  2:5785278/93‑1 (MQ=255)
                                                ccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTgg                                         >  2:3094316/1‑93 (MQ=255)
                                                ccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTgg                                         >  1:858414/1‑93 (MQ=255)
                                                ccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTgg                                         >  1:2742791/1‑93 (MQ=255)
                                                ccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTgg                                         >  2:2508453/1‑93 (MQ=255)
                                                 cATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGt                                        >  1:3933844/1‑93 (MQ=255)
                                                  aTTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGtt                                       >  1:1527560/1‑93 (MQ=255)
                                                  aTTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGtt                                       >  2:2566205/1‑93 (MQ=255)
                                                    tttCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTttg                                     <  1:2771404/93‑1 (MQ=255)
                                                    tttCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTttg                                     <  1:2070549/93‑1 (MQ=255)
                                                      tCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTttgtt                                   <  1:4475186/93‑1 (MQ=255)
                                                       ccTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTttgttg                                  >  2:2511534/1‑93 (MQ=255)
                                                        cTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCGCAATGTTCGCCATTTAGACCTTTGGCCTGCAAATCGTTGGtgtggttg                                 >  1:3985859/1‑89 (MQ=255)
                                                        cTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTgg                                 >  2:2280338/1‑93 (MQ=255)
                                                         ttttGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTggg                                >  2:5645823/1‑93 (MQ=255)
                                                         ctgtGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTggg                                <  2:3185907/90‑1 (MQ=255)
                                                           ttGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGAt                              >  1:213841/1‑93 (MQ=255)
                                                           ttGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGAt                              <  2:413815/93‑1 (MQ=255)
                                                           ttGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGAt                              >  1:464578/1‑93 (MQ=255)
                                                           ttGTTCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGAt                              >  2:4403463/1‑93 (MQ=255)
                                                             tttCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTg                            <  1:6128014/92‑1 (MQ=255)
                                                             gTTCTGTTTGTTGCCTCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCGTTGGCCCGCTAATCGCTGTCTTGTTCTGATTg                            >  1:869132/1‑93 (MQ=255)
                                                              ttCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTgg                           <  1:4147052/93‑1 (MQ=255)
                                                               tCTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAGATCGTTGGTTTGTTGGGATTggg                          >  1:5501599/1‑93 (MQ=255)
                                                                cTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTgggg                         <  2:3387202/93‑1 (MQ=255)
                                                                cTGTTTGTTGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTgggg                         <  2:1790878/93‑1 (MQ=255)
                                                                  gtttgttGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGt                       >  1:3538892/1‑93 (MQ=255)
                                                                   tttgttGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGtt                      <  2:1594927/93‑1 (MQ=255)
                                                                   tttgttGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGtt                      >  1:2231680/1‑93 (MQ=255)
                                                                   tttgttGCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGtt                      >  2:1988191/1‑93 (MQ=255)
                                                                    ttgttgCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTg                     <  1:1060401/93‑1 (MQ=255)
                                                                        tgCCGCCATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCaa                 >  2:5711944/1‑93 (MQ=255)
                                                                          ccgccATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCaaaa               <  2:1125238/93‑1 (MQ=255)
                                                                           cgccATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCaaaaa              <  1:5240054/93‑1 (MQ=255)
                                                                           cgccATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCaaaaa              <  1:483504/93‑1 (MQ=255)
                                                                             ccATCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAAcg            >  2:2300801/1‑93 (MQ=255)
                                                                               aTCTGGACCTACGAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAAcgcg          <  1:257781/93‑1 (MQ=255)
                                                                               aTCTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAAcgcg          >  1:600065/1‑93 (MQ=255)
                                                                                 cTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTTGGTTGTTGGGATTGGGGGGTGGCAAAAAcgctga        >  2:5412811/1‑90 (MQ=255)
                                                                                 cTGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGc        >  2:3820209/1‑93 (MQ=255)
                                                                                  tGGACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCa       >  1:1698498/1‑93 (MQ=255)
                                                                                    gACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACt     <  2:524878/93‑1 (MQ=255)
                                                                                    gACCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACt     >  2:1731110/1‑93 (MQ=255)
                                                                                     gccTACCAGGTGATCCCAAGGTTCGCCATTTACCCCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTa    <  2:3728276/92‑1 (MQ=255)
                                                                                     aCCTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTa    <  2:48223/93‑1 (MQ=255)
                                                                                      ccTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAg   <  2:4063924/93‑1 (MQ=255)
                                                                                      ccTACCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAg   >  1:5083043/1‑93 (MQ=255)
                                                                                         |                                                                                          
TGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGG  >  NC_000913/46478‑46657

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