Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,341,621 G→T 6.5% L266L (CTG→CTT yraO → transporter YraO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,341,6210GT 6.5% 195.7 / 2.5 61L266L (CTG→CTTyraOtransporter YraO
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base T (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

TACATCAAGGAAATCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGATGGTCGATGCAT  >  NC_002947/2341472‑2341760
                                                                                                                                                     |                                                                                                                                           
tACATCAAGGAAATCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtg                                                                                                                                             >  2:405866/1‑150 (MQ=255)
 aCATCAAGGAAATCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtga                                                                                                                                            <  1:405866/150‑1 (MQ=255)
   aTCAAGGAAATCCTCGGCGACACTCCTCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtgatg                                                                                                                                          <  2:379158/150‑1 (MQ=255)
       aGGAAATCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtgat                                                                                                                                           <  1:360995/145‑1 (MQ=255)
       aGGAAATCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtgat                                                                                                                                           >  2:360995/1‑145 (MQ=255)
         gAAATCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAAc                                                                                                                                    >  2:1376487/1‑150 (MQ=255)
          aaaTCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACt                                                                                                                                   >  1:770709/1‑150 (MQ=255)
            aTCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTAc                                                                                                                                 <  2:628984/150‑1 (MQ=255)
               cTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCaa                                                                                                                                     >  1:1060118/1‑143 (MQ=255)
               cTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCaa                                                                                                                                     <  2:1060118/143‑1 (MQ=255)
                 cGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCa                                                                                                                            >  2:1137943/1‑150 (MQ=255)
                     gACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtgatg                                                                                                                                          >  1:143282/1‑132 (MQ=255)
                     gACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtgatg                                                                                                                                          <  2:143282/132‑1 (MQ=255)
                           ccGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtgat                                                                                                                                           >  1:1147046/1‑125 (MQ=255)
                           ccGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACtgat                                                                                                                                           <  2:1147046/125‑1 (MQ=255)
                                 aaGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCCGGGTACGGTCCTGATGATGGTCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGaaa                                                                                                            >  1:28870/1‑150 (MQ=255)
                                 aaGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGaaa                                                                                                            <  1:33077/150‑1 (MQ=255)
                                   gcgcCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAaga                                                                                                          >  1:839400/1‑150 (MQ=255)
                                   gcgcCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAaga                                                                                                          >  2:907264/1‑150 (MQ=255)
                                     gcCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAgc                                                                                                        <  2:322202/150‑1 (MQ=255)
                                       ccGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCATCTACCCGCAGCTGGACCTGCAGAAAGAgcgc                                                                                                      >  1:1189467/1‑150 (MQ=255)
                                       ccGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAgcgc                                                                                                      >  2:1432784/1‑150 (MQ=255)
                                             cTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTc                                                                                                <  2:770709/150‑1 (MQ=255)
                                                gCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGcc                                                                                             <  1:1137943/150‑1 (MQ=255)
                                                  cTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCa                                                                                           >  2:186074/1‑150 (MQ=255)
                                                    aTGTGAACCTGGTGCTGGTGATCGCACCGCTGAACTGGAGTTTGGTGGTGTTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCAtt                                                                                         <  1:1598922/150‑1 (MQ=255)
                                                         aaCCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGc                                                                                    >  1:1101416/1‑150 (MQ=255)
                                                          aCCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCggg                                                                                     >  1:320767/1‑148 (MQ=255)
                                                          aCCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCggg                                                                                     <  2:320767/148‑1 (MQ=255)
                                                                     gtgATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGAcc                                                                        >  1:1287740/1‑150 (MQ=255)
                                                                      tgATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCg                                                                       >  2:1116568/1‑150 (MQ=255)
                                                                        aTCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCgtg                                                                     <  1:1116568/150‑1 (MQ=255)
                                                                         tCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCgtgg                                                                    >  1:1218847/1‑150 (MQ=255)
                                                                            ccGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGtgc                                                                 >  1:1552332/1‑150 (MQ=255)
                                                                               tgatgaCCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGtgctgc                                                              <  1:1551255/150‑1 (MQ=255)
                                                                                gatgaCCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTgctgct                                                             >  1:65820/1‑150 (MQ=255)
                                                                                 atgaCCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGctgctg                                                            <  2:839400/150‑1 (MQ=255)
                                                                                   gaCCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGgt                                                          <  2:1552332/150‑1 (MQ=255)
                                                                                     ccGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGAc                                                                                                                     <  1:152570/89‑1 (MQ=255)
                                                                                     ccGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGAc                                                                                                                     >  2:152570/1‑89 (MQ=255)
                                                                                     ccGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATAATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGgtgt                                                        >  1:378502/1‑150 (MQ=255)
                                                                                         gTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCa                                                                             >  2:733978/1‑125 (MQ=255)
                                                                                         gTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCa                                                                             <  1:733978/125‑1 (MQ=255)
                                                                                              tGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCgg                                               <  2:1519088/150‑1 (MQ=255)
                                                                                              tGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCgg                                               <  1:186074/150‑1 (MQ=255)
                                                                                               gCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCggg                                              >  2:1280170/1‑150 (MQ=255)
                                                                                                cTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGc                                             <  1:1376487/150‑1 (MQ=255)
                                                                                                   gggCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCgg                                               >  1:683143/1‑145 (MQ=255)
                                                                                                   gggCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCgg                                               <  2:683143/145‑1 (MQ=255)
                                                                                                     gCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATAATAGCGGAAGAGGGCAGCGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCtt                                        >  2:458007/1‑150 (MQ=255)
                                                                                                        gATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGc                                     >  1:1367842/1‑150 (MQ=255)
                                                                                                               gcgGCCATCCTGTTCATGATCGGTTTCGTCGCCGCAGTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTAtt                              <  2:385913/150‑1 (MQ=255)
                                                                                                                    cATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTa                                <  2:1453625/143‑1 (MQ=255)
                                                                                                                    cATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTa                                >  1:1453625/1‑143 (MQ=255)
                                                                                                                    cATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGAGAAGCTCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTtctc                         >  1:430859/1‑150 (MQ=255)
                                                                                                                      tCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCgg                       <  1:907264/150‑1 (MQ=255)
                                                                                                                      tCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCgg                       <  2:1218847/150‑1 (MQ=255)
                                                                                                                              aTGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGAc                                                                                                                     <  1:688340/48‑1 (MQ=255)
                                                                                                                              aTGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGAc                                                                                                                     >  2:688340/1‑48 (MQ=255)
                                                                                                                                gATCGGTTTCGTCGCCGCACTTATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAgc                                                                                                        <  1:630599/59‑1 (MQ=255)
                                                                                                                                gATCGGTTTCGTCGCCGCACTTATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAgc                                                                                                        >  2:630599/1‑59 (MQ=255)
                                                                                                                                gATCGGTTTCGTCGCCGCACTTATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAgc                                                                                                        >  2:630617/1‑59 (MQ=255)
                                                                                                                                gATCGGTTTCGTCGCCGCACTTATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAgc                                                                                                        <  1:630617/59‑1 (MQ=255)
                                                                                                                                gATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGAt             >  2:292893/1‑150 (MQ=255)
                                                                                                                                   cGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGATGGt          <  2:65820/150‑1 (MQ=255)
                                                                                                                                    ggTTTCGTCGCCGCACTGATGATCAACTAACCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGATGGTc         >  2:1102584/1‑150 (MQ=255)
                                                                                                                                    gctttCGTCGTCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTAGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGATGGTc         <  1:908960/148‑1 (MQ=255)
                                                                                                                                      tttCGTCGCCGCACTGATGATCATCTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCACGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGATGGTCGa       <  2:1189467/150‑1 (MQ=255)
                                                                                                                                      tttCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGATGGTCGa       <  2:1287740/150‑1 (MQ=255)
                                                                                                                                           tcgCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGATGGTCGATGCAt  <  2:1367842/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
TACATCAAGGAAATCCTCGGCGACAACCCGCACAAGCGCCCGCGCCTTGCCTATGTGAACCTGGTGCTGGTGATCGCCGTGATGACCGCGTTGGTGCTGGGGCTGATGCACGCGGCCATCCTGTTCATGATCGGTTTCGTCGCCGCACTGATGATCAACTACCCGCAGCTGGACCTGCAGAAAGAGCGCATCCTCGCCCATTCGGGCAACGCCATGACCGTGGTGCTGCTGGTGTTCGCCGCGGGCATCTTTGCCGGTATTTTCTCGGGTACCAAGATGGTCGATGCAT  >  NC_002947/2341472‑2341760

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