Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 868,422 (C)8→7 coding (903/1872 nt) gsiA → glutathione ABC transporter ATP binding subunit GsiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913868,4150C.98.4% 251.9 / ‑2.1 61coding (896/1872 nt)gsiAglutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (24/36);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGT  >  NC_000913/868279‑868547
                                                                                                                                        |                                                                                                                                    
gAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccc                                                                                                                                   >  1:1300404/1‑140 (MQ=255)
 aaCAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccc                                                                                                                                  >  2:231225/1‑140 (MQ=255)
      aTTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTCCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccccccatc                                                                                                                             <  1:820161/140‑4 (MQ=255)
        ttttCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGa                                                                                                                          <  2:154408/140‑1 (MQ=255)
         tttCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAg                                                                                                                         <  1:221277/140‑1 (MQ=255)
            cATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCa                                                                                                                       <  1:1176073/140‑1 (MQ=255)
                cACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGaaaa                                                                                                                  >  2:54650/1‑140 (MQ=255)
                cACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGaaaa                                                                                                                  <  2:451981/140‑1 (MQ=255)
                       aCATCCTTACCCCGGTGCGCTGTTAGCTGCTGTTCCGCAACTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGgtggtg                                                                                                          <  1:470695/140‑1 (MQ=255)
                             ttACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATgg                                                                                                     <  1:1410351/140‑1 (MQ=255)
                               acacCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCg                                                                                                   >  1:691180/1‑140 (MQ=255)
                                 acCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 <  2:791003/140‑1 (MQ=255)
                                   ccGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  1:215783/140‑1 (MQ=255)
                                   ccGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  1:836221/140‑1 (MQ=255)
                                   ccGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAcc                                                                                               <  2:1137173/140‑1 (MQ=255)
                                    cGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACggcgaa                                                                                                           >  1:943320/1‑123 (MQ=255)
                                    cGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACggcgaa                                                                                                           <  2:943320/127‑5 (MQ=255)
                                    cGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCt                                                                                              >  2:459702/1‑140 (MQ=255)
                                    cGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGGGGTGGATGGCGAACCt                                                                                              >  2:92286/1‑140 (MQ=255)
                                     gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATgg                                                                                                     <  1:431193/132‑1 (MQ=255)
                                     gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATgg                                                                                                     >  2:431193/1‑132 (MQ=255)
                                     gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTg                                                                                             >  1:411089/1‑140 (MQ=255)
                                      tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 <  1:1101846/135‑1 (MQ=255)
                                      tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                                                                 >  2:1101846/1‑135 (MQ=255)
                                      tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                                                            >  1:583796/1‑140 (MQ=255)
                                      tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                                                            >  1:1249624/1‑140 (MQ=255)
                                      tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                                                            >  1:888574/1‑140 (MQ=255)
                                       gcgcTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                                                                           >  2:300252/1‑140 (MQ=255)
                                       gcgcTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                                                                           <  2:1285277/140‑1 (MQ=255)
                                        cgcTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGttt                                                                                          >  1:893316/1‑140 (MQ=255)
                                            gTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACg                                                                                      <  1:744356/140‑1 (MQ=255)
                                            gTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACg                                                                                      <  1:509514/140‑1 (MQ=255)
                                            gTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACg                                                                                      <  2:1249624/140‑1 (MQ=255)
                                             ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                                                                     >  1:69472/1‑140 (MQ=255)
                                             ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                                                                     <  2:143172/140‑1 (MQ=255)
                                                 ctgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGc                                                                                 <  1:300252/140‑1 (MQ=255)
                                                    ctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                                                              <  1:54650/140‑1 (MQ=255)
                                                    ctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                                                              <  2:69472/140‑1 (MQ=255)
                                                       ttCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATc                                                                           >  1:1053329/1‑140 (MQ=255)
                                                       ttCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTCTACGAGTGCGTAATc                                                                           >  1:803492/1‑140 (MQ=255)
                                                         ccGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCtt                                                                         <  2:691180/140‑1 (MQ=255)
                                                         ccGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCtt                                                                         >  1:1324341/1‑140 (MQ=255)
                                                                   gTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTCTACGAGTGCGTAATCTTGTCACCCGtt                                                               <  2:803492/140‑1 (MQ=255)
                                                                                  tAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                <  1:459702/140‑1 (MQ=255)
                                                                                  tAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                <  1:1123968/140‑1 (MQ=255)
                                                                                  tAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg                                                <  1:330601/140‑1 (MQ=255)
                                                                                     aTTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGtt                                             >  2:280883/1‑140 (MQ=255)
                                                                                     aTTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGtt                                             >  2:495763/1‑140 (MQ=255)
                                                                                                    tCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTaa                              <  2:1127117/140‑1 (MQ=255)
                                                                                                     cccGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAAc                             >  2:211571/1‑140 (MQ=255)
                                                                                                     cccGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAAc                             <  1:211571/140‑1 (MQ=255)
                                                                                                           gATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGa                       <  2:1324341/140‑1 (MQ=255)
                                                                                                              atCGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGt                    >  1:201846/1‑140 (MQ=255)
                                                                                                                cGCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGc                  <  1:1104428/140‑1 (MQ=255)
                                                                                                                 gCTTGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCa                 <  2:411089/140‑1 (MQ=255)
                                                                                                                   ttGAACATCCAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATg               <  2:893316/140‑1 (MQ=255)
                                                                                                                           ccAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTgag       >  2:804696/1‑140 (MQ=255)
                                                                                                                            cAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCTGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCg            >  2:1319515/1‑134 (MQ=255)
                                                                                                                            cAGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCg            <  1:1319515/134‑1 (MQ=255)
                                                                                                                             aGCGAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGaa     >  1:512524/1‑140 (MQ=255)
                                                                                                                                gAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  1:495763/140‑1 (MQ=255)
                                                                                                                                gAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  2:888574/140‑1 (MQ=255)
                                                                                                                                gAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  2:201846/140‑1 (MQ=255)
                                                                                                                                gAAACAGGCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGt  <  1:280883/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
GAACAGATTTTTCATGCACCGCAACATCCTTACACCCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTGTTGAATCGCGTAACGCGGGAAGTGCATGCCGTTGAGAAAGT  >  NC_000913/868279‑868547

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