Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,653,200 A→G 5.1% M591V (ATG→GTG)  fbp → Fructose‑1,6‑bisphosphatase class 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,653,2000AG 5.1% 695.2 / 15.0 216M591V (ATG→GTG) fbpFructose‑1,6‑bisphosphatase class 3
Reads supporting (aligned to +/- strand):  ref base A (120/85);  new base G (2/9);  total (122/94)
Fisher's exact test for biased strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 219 total aligned reads displayed.

AAGATGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGCGTCGTAGATA  >  USA300TCH1516_ALE/2653110‑2653290
                                                                                          |                                                                                          
aaGATGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTa                                                                                            >  1:921519/1‑91 (MQ=255)
aaGATGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTa                                                                                            >  2:3511060/1‑91 (MQ=255)
  gATGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATg                                                                                          >  2:285933/1‑91 (MQ=255)
  gATGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATg                                                                                          >  2:286856/1‑91 (MQ=255)
   aTGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGc                                                                                         <  2:528306/91‑1 (MQ=255)
   aTGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGc                                                                                         <  2:2175701/91‑1 (MQ=255)
    tGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCa                                                                                        >  1:844520/1‑91 (MQ=255)
    tGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCa                                                                                        >  1:2110582/1‑91 (MQ=255)
    tGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCa                                                                                        >  2:2170250/1‑91 (MQ=255)
     gtttGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCaa                                                                                       <  2:275961/89‑1 (MQ=255)
      cTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAAc                                                                                      >  2:552213/1‑91 (MQ=255)
       ttGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACt                                                                                     >  1:1147305/1‑91 (MQ=255)
        tGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAAcga                                                                                    <  2:3858624/91‑3 (MQ=255)
        tGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTc                                                                                    <  1:2785904/91‑1 (MQ=255)
        tGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTc                                                                                    <  1:419517/91‑1 (MQ=255)
        tGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATGGTATAATTCATTCGGTATGCAACTc                                                                                    <  2:2701754/91‑1 (MQ=255)
         gTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTGTGCAAcgag                                                                                   <  1:3561135/91‑4 (MQ=255)
         gTCATTGATGGCGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACATTGTTGTATAATTCATTCGGAATGCAACTCg                                                                                   <  2:1161636/91‑1 (MQ=255)
          tCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGt                                                                                  >  2:1389704/1‑91 (MQ=255)
           cATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGtt                                                                                 >  2:3547796/1‑91 (MQ=255)
           cATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGtt                                                                                 >  1:3570559/1‑91 (MQ=255)
             ttGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTGTGCAACGAGTTGc                                                                               <  2:1420159/91‑1 (MQ=255)
              tGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGca                                                                              >  2:1668502/1‑91 (MQ=255)
                   gtggCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACca                                                                         <  2:244264/91‑1 (MQ=255)
                    tggCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACcaa                                                                        <  1:53097/91‑1 (MQ=255)
                     ggCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACcaac                                                                       <  1:239268/91‑1 (MQ=255)
                     ggCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACcaac                                                                       >  2:772779/1‑91 (MQ=255)
                         tttCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTCTTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAAtt                                                                   >  1:3074783/1‑91 (MQ=255)
                         tttCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAAtt                                                                   >  1:145584/1‑91 (MQ=255)
                             aaaaGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAAc                                                               >  1:2387248/1‑91 (MQ=255)
                             aaaaGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAAc                                                               >  2:1282945/1‑91 (MQ=255)
                              aaaGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACg                                                              <  2:1051111/91‑1 (MQ=255)
                               aaGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGc                                                             <  2:1517026/91‑1 (MQ=255)
                               aaGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGc                                                             <  1:3821819/91‑1 (MQ=255)
                                 gCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTa                                                           >  2:3077642/1‑91 (MQ=255)
                                 gCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTa                                                           >  2:713658/1‑91 (MQ=255)
                                 gCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTa                                                           >  2:906242/1‑91 (MQ=255)
                                   ttATCAGTCAACAACTGGCCTTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTaaa                                                         >  2:1026670/1‑91 (MQ=255)
                                   ttATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTaaa                                                         >  1:1084217/1‑91 (MQ=255)
                                   ttATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTaaa                                                         >  1:1083974/1‑91 (MQ=255)
                                   ttATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTaaa                                                         >  1:3255457/1‑91 (MQ=255)
                                    tATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAg                                                        >  2:1561214/1‑91 (MQ=255)
                                    tATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAg                                                        >  2:48883/1‑91 (MQ=255)
                                     aTCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGa                                                       <  1:1557646/91‑1 (MQ=255)
                                      tCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGaa                                                      >  2:159283/1‑91 (MQ=255)
                                      tCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGaa                                                      >  1:2212362/1‑91 (MQ=255)
                                      tCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGaa                                                      >  2:859866/1‑91 (MQ=255)
                                       cAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGaaa                                                     >  1:2544536/1‑91 (MQ=255)
                                       cAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGAACACCAACAATTTAACGCTAAAGaaa                                                     >  1:1871256/1‑91 (MQ=255)
                                        aGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGaaaa                                                    <  2:1951141/91‑1 (MQ=255)
                                         gTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGaaaaa                                                   <  1:463020/91‑1 (MQ=255)
                                          tCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGaaaaaa                                                  <  1:2660821/91‑1 (MQ=255)
                                            aacaacTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATa                                                >  2:3166495/1‑91 (MQ=255)
                                             acaacTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAAAAAAAAAAAAAAAATAc                                               <  2:2808777/91‑1 (MQ=255)
                                             acaacTGGCATTGCAGGCTATACACGATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATAc                                               <  1:3600927/91‑1 (MQ=255)
                                              caacTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACt                                              <  2:2065164/91‑1 (MQ=255)
                                                acTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACttt                                            >  1:2551580/1‑91 (MQ=255)
                                                 cTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACttcc                                           >  1:3796320/1‑89 (MQ=255)
                                                 cTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTc                                           >  1:2123987/1‑91 (MQ=255)
                                                 cTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTc                                           <  1:1097938/91‑1 (MQ=255)
                                                  tGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTcc                                          <  2:3470893/91‑1 (MQ=255)
                                                   ggCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCg                                         >  1:2663573/1‑91 (MQ=255)
                                                    gCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGa                                        >  2:902141/1‑91 (MQ=255)
                                                       ttGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAgg                                     >  1:3721407/1‑91 (MQ=255)
                                                        tGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGt                                    >  2:360648/1‑91 (MQ=255)
                                                            ggCTATATACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAAAAAAAATACTTTCCGAAGGTCTCg                                <  2:3440247/91‑1 (MQ=255)
                                                            aattataCAATAGTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCg                                >  1:1837825/4‑91 (MQ=255)
                                                             gCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGa                               >  1:3234334/1‑91 (MQ=255)
                                                                ataCACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGa                            >  1:3423784/1‑91 (MQ=255)
                                                                    acTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACtt                        >  2:1177657/1‑91 (MQ=255)
                                                                    acTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACtt                        <  1:2616458/91‑1 (MQ=255)
                                                                      tATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTc                      >  1:1903201/1‑91 (MQ=255)
                                                                      tATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTc                      >  1:256135/1‑91 (MQ=255)
                                                                           tataATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATg                             >  2:966494/1‑79 (MQ=255)
                                                                           tataATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATg                             <  1:966494/79‑1 (MQ=255)
                                                                            ataATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATaaa                <  2:2187689/91‑1 (MQ=255)
                                                                            ataATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATaaa                <  1:1967638/91‑1 (MQ=255)
                                                                            ataATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATaaa                <  2:2002930/91‑1 (MQ=255)
                                                                            ataATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATaaa                <  2:1674648/91‑1 (MQ=255)
                                                                             taATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAg               <  2:3244343/91‑1 (MQ=255)
                                                                              aaTTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAgc              <  1:693749/91‑1 (MQ=255)
                                                                              aaTTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAgc              <  1:285933/91‑1 (MQ=255)
                                                                              aaTTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAgc              <  1:286856/91‑1 (MQ=255)
                                                                               attcattcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAGAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAgcg             >  2:3417449/1‑91 (MQ=255)
                                                                               attcattcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAgcg             >  2:2581212/1‑91 (MQ=255)
                                                                               attcattcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAgcg             >  1:1648018/1‑91 (MQ=255)
                                                                                tccattcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAgcgc            <  1:1424202/89‑1 (MQ=255)
                                                                                 tcattcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGcgcg           >  1:822375/1‑91 (MQ=255)
                                                                                 tcattcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGcgcg           >  1:1304977/1‑91 (MQ=255)
                                                                                   attcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGcgtc         >  2:2051326/1‑91 (MQ=255)
                                                                                    ttcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGcgtcg        >  1:98505/1‑91 (MQ=255)
                                                                                    ttcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGcgtcg        <  2:3775579/91‑1 (MQ=255)
                                                                                     tcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGcgtcgt       >  1:2790536/1‑91 (MQ=255)
                                                                                     tcGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGcgtcgt       <  1:1667830/91‑1 (MQ=255)
                                                                                       ggTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGCGTCGTAg     >  1:1212136/1‑91 (MQ=255)
                                                                                        gTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGCGTCGTAGa    >  2:1888400/1‑91 (MQ=255)
                                                                                        gTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGCGTCGTAGa    <  2:1147305/91‑1 (MQ=255)
                                                                                         tATGCAACTCGTTGGACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGCGTCGTAGAt   >  2:3295605/1‑91 (MQ=255)
                                                                                         tATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGCGTCGTAGAt   <  1:5842/91‑1 (MQ=255)
                                                                                          aTGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGCGTCGTAGATa  <  1:1707047/91‑1 (MQ=255)
                                                                                          |                                                                                          
AAGATGCTTGTCATTGATGGTGGCTTTTCAAAAGCTTATCAGTCAACAACTGGCATTGCAGGCTATACACTATTGTATAATTCATTCGGTATGCAACTCGTTGCACACCAACAATTTAACGCTAAAGAAAAAATACTTTCCGAAGGTATCGATGAACTTTCTATAAAGCGCGTCGTAGATA  >  USA300TCH1516_ALE/2653110‑2653290

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