Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,685,130 A→G 7.5% intergenic (‑142/+78) USA300HOU_RS08395 ← / ← rpsU hypothetical protein/30S ribosomal protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,685,1300AG 7.5% 580.8 / 25.1 187intergenic (‑142/+78)USA300HOU_RS08395/rpsUhypothetical protein/30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base A (90/82);  new base G (5/9);  total (96/91)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 195 total aligned reads displayed.

AGTTATAAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTT  >  USA300TCH1516_ALE/1685040‑1685219
                                                                                          |                                                                                         
aGTTATAAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACata                                                                                           <  1:2568349/91‑1 (MQ=255)
 gTTATAAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACatat                                                                                          <  2:213208/91‑1 (MQ=255)
 gTTATAAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACatat                                                                                          >  1:2304925/1‑91 (MQ=255)
  ttATAAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATAtt                                                                                         <  1:1606491/91‑1 (MQ=255)
    ataAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAg                                                                                       >  1:981143/1‑91 (MQ=255)
    ataAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAg                                                                                       >  2:963350/1‑91 (MQ=255)
    ataAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAg                                                                                       >  2:1926848/1‑91 (MQ=255)
    ataAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAg                                                                                       >  1:2707175/1‑91 (MQ=255)
     taAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAgt                                                                                      <  2:1821886/91‑1 (MQ=255)
      aaGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAgtg                                                                                     <  1:32771/91‑1 (MQ=255)
      aaGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAgtg                                                                                     >  2:3170297/1‑91 (MQ=255)
      aaGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAgtg                                                                                     >  1:2238268/1‑91 (MQ=255)
        gAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTc                                                                                   <  1:144016/91‑1 (MQ=255)
         aaCCTATAGTATATCAATTTGATTATTTATAAACTATATATCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTct                                                                                  <  2:3873191/91‑2 (MQ=255)
          aCCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAGATTAGTGTCAg                                                                                 >  1:3515032/1‑91 (MQ=255)
          aCCGATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTCGTGTCAg                                                                                 >  2:3439946/1‑91 (MQ=255)
           ccTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGa                                                                                >  2:3374537/1‑91 (MQ=255)
           ccTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGa                                                                                >  1:662352/1‑91 (MQ=255)
           ccTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGa                                                                                <  1:2057299/91‑1 (MQ=255)
           ccTAGAGAATATCAATGTGATTATTTATAAACTATATAACTTTTTTTAAATAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGa                                                                                <  2:3535152/91‑1 (MQ=255)
            cTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGac                                                                               <  2:3280104/91‑1 (MQ=39)
              ataGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAATATTGACAACAAAAAAACATCACATATTAGTGTCAGacac                                                                             >  2:1152518/1‑91 (MQ=255)
              ataGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGacac                                                                             <  2:79893/91‑1 (MQ=255)
               taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                            >  1:60109/1‑91 (MQ=255)
               taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                            >  1:2811547/1‑91 (MQ=255)
               taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                            >  1:3535605/1‑91 (MQ=255)
               taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                            <  2:2709095/91‑1 (MQ=255)
                aGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTa                                                                           <  1:2360378/91‑1 (MQ=255)
                 gAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTaa                                                                          >  2:1973016/1‑91 (MQ=255)
                 gAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTaa                                                                          <  1:391593/91‑1 (MQ=255)
                  aaTATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAAc                                                                         >  2:1329100/1‑91 (MQ=255)
                  aaTATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAAc                                                                         >  2:3707914/1‑91 (MQ=255)
                  aaTATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAAc                                                                         >  1:1876607/1‑91 (MQ=255)
                   atatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACt                                                                        <  1:2196352/91‑1 (MQ=255)
                    tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAACAAACATCACAAGTTAGTGTCTGACACTaatat                                                                       <  2:3234340/91‑4 (MQ=25)
                    tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                       >  1:2200398/1‑91 (MQ=255)
                      tCAATTTGACTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                     <  1:3711886/91‑1 (MQ=255)
                         aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                  <  1:1972530/91‑1 (MQ=255)
                         aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                  >  2:2412331/1‑91 (MQ=255)
                          tttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtg                                                                 >  2:1078387/1‑91 (MQ=255)
                           ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                <  2:2846160/91‑1 (MQ=255)
                           ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                >  1:2646473/1‑91 (MQ=255)
                           ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgt                                                                <  2:1621450/91‑1 (MQ=21)
                            tGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtt                                                               <  1:2484165/91‑1 (MQ=255)
                              attattTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttg                                                             >  2:2201339/1‑91 (MQ=255)
                                tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAACAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgtttttt                                                           <  2:2372221/91‑4 (MQ=17)
                                tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                           >  1:424290/1‑91 (MQ=255)
                                tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                           <  2:100752/91‑1 (MQ=255)
                                 atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                          <  2:2945187/91‑1 (MQ=255)
                                   ttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTta                                                        <  1:3520068/91‑1 (MQ=255)
                                       aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                    >  1:821142/1‑91 (MQ=255)
                                       aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                    >  2:3630612/1‑91 (MQ=255)
                                       aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                    >  1:1028962/1‑91 (MQ=255)
                                         aCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaat                                                  >  2:3689631/1‑91 (MQ=255)
                                          cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                 >  1:959862/1‑91 (MQ=255)
                                            atataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaa                                               <  2:1058499/91‑1 (MQ=255)
                                             tataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                              <  1:2367910/91‑1 (MQ=255)
                                              ataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAata                                             >  1:3189337/1‑91 (MQ=255)
                                                aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                           >  2:1506254/1‑91 (MQ=255)
                                                 aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                          <  2:2097135/91‑1 (MQ=255)
                                                   tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                        >  2:2671941/1‑91 (MQ=255)
                                                   tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                        <  1:1885028/91‑1 (MQ=255)
                                                   tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                        <  2:690710/91‑1 (MQ=255)
                                                   tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                        >  1:2492539/1‑91 (MQ=255)
                                                    ttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtg                                       >  1:1555053/1‑91 (MQ=255)
                                                     tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                      <  2:504109/91‑1 (MQ=255)
                                                     tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                      <  2:731145/91‑1 (MQ=255)
                                                     tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                      <  2:1337013/91‑1 (MQ=255)
                                                     tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                      <  2:3025846/91‑1 (MQ=255)
                                                        ttaattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGa                                   <  2:672020/91‑1 (MQ=255)
                                                           attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                <  1:1518933/91‑4 (MQ=255)
                                                           attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                <  1:1457154/91‑4 (MQ=255)
                                                           attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                >  1:549567/1‑91 (MQ=255)
                                                             taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                              <  2:1537002/91‑1 (MQ=255)
                                                              aaaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGt                             <  2:179900/91‑1 (MQ=255)
                                                              aaaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGt                             <  2:1733172/91‑1 (MQ=255)
                                                                 aaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCaa                          <  1:1053155/91‑1 (MQ=255)
                                                                 aaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCaa                          <  1:3766841/91‑1 (MQ=255)
                                                                   ttGACAACAAAAAAACATCACATACTAGTGCCAGACACTAACTTGTGATGCTCGTTTTATATAATTAATATTCGTGTTGTGGGACTCAAcc                        >  2:1876374/1‑90 (MQ=255)
                                                                    tGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACag                       <  1:3646827/91‑1 (MQ=255)
                                                                     gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                      >  1:3402684/1‑91 (MQ=255)
                                                                     gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                      >  1:3401826/1‑91 (MQ=255)
                                                                      acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                     >  1:2013209/1‑91 (MQ=255)
                                                                      acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                     <  1:2412331/91‑1 (MQ=255)
                                                                       caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                    >  1:351905/1‑91 (MQ=255)
                                                                         acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                  >  1:2512042/1‑91 (MQ=255)
                                                                         acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAAGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                  >  2:2853538/1‑91 (MQ=255)
                                                                          caaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                 <  1:1163305/91‑1 (MQ=255)
                                                                           aaaaaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTatt                >  2:2980108/1‑91 (MQ=255)
                                                                               aaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatt            >  2:2753457/1‑91 (MQ=255)
                                                                                 aCATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattat          >  2:2344484/1‑91 (MQ=25)
                                                                                  cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt         <  1:169656/91‑1 (MQ=255)
                                                                                  cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt         <  1:243586/91‑1 (MQ=255)
                                                                                   aTCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTAttt        >  1:1345593/1‑91 (MQ=25)
                                                                                    tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg       >  2:1042001/1‑91 (MQ=255)
                                                                                    tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg       >  2:2855094/1‑91 (MQ=255)
                                                                                     cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa      >  2:3798167/1‑91 (MQ=255)
                                                                                      acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa     >  1:3327893/1‑91 (MQ=255)
                                                                                      acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa     <  1:3707914/91‑1 (MQ=255)
                                                                                        atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt   >  2:2486038/1‑91 (MQ=255)
                                                                                          |                                                                                         
AGTTATAAGAACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTT  >  USA300TCH1516_ALE/1685040‑1685219

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