Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 917,330 C→A 7.8% intergenic (+89/+241) sufB_2 → / ← USA300HOU_RS04545 FeS cluster assembly protein SufB/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE917,3300CA 7.8% 772.8 / 36.1 226intergenic (+89/+241)sufB_2/USA300HOU_RS04545FeS cluster assembly protein SufB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (115/93);  new base A (5/13);  total (120/106)
Fisher's exact test for biased strand distribution p-value = 2.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

AGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAAAA  >  USA300TCH1516_ALE/917240‑917419
                                                                                          |                                                                                         
  cTTTAAGCCGCGTTGTTAAGCCATTATTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAGCAAACAGCCACAGTAAtg                                                                                         <  2:942181/91‑1 (MQ=255)
  cTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACtg                                                                                         >  2:1072647/1‑91 (MQ=255)
    ttAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACAGTAAtgtg                                                                                       <  2:3666391/91‑1 (MQ=255)
     tAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACAGTAATGTgg                                                                                      <  1:1338316/91‑1 (MQ=255)
     tAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACAGTAATGTgg                                                                                      <  1:204270/91‑1 (MQ=255)
      aaaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGc                                                                                     >  1:1526021/1‑91 (MQ=255)
       aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCt                                                                                    <  1:298846/91‑1 (MQ=255)
       aaCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCt                                                                                    >  1:1991727/1‑91 (MQ=255)
        aCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCtt                                                                                   >  2:1050336/1‑91 (MQ=255)
        aCCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCtt                                                                                   <  2:1390634/91‑1 (MQ=255)
         ccGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttt                                                                                  <  1:1167/91‑1 (MQ=255)
         ccGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttt                                                                                  >  1:2605748/1‑91 (MQ=255)
           gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttttt                                                                                <  1:569872/91‑1 (MQ=255)
           gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttttt                                                                                <  1:3815239/91‑1 (MQ=255)
           gcgTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttttt                                                                                >  1:1355663/1‑91 (MQ=255)
             gttgttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCttttttt                                                                              <  1:3516211/91‑1 (MQ=255)
                gttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGtt                                                                           >  2:1126008/1‑91 (MQ=255)
                gttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCCTTTTTTgct                                                                           >  2:2725028/1‑89 (MQ=255)
                gatAAGCCATTCTTGACTGCCGAAAATGGCTATTGATACCATTTTGCTACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGtt                                                                           >  2:1385073/3‑91 (MQ=255)
                 ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCATATTACTGTGGCTTTTTTTGttt                                                                          >  2:3120305/1‑91 (MQ=255)
                 ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGttt                                                                          >  2:2121106/1‑91 (MQ=255)
                 ttAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGttt                                                                          >  1:741574/1‑91 (MQ=255)
                   aaGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTta                                                                        >  2:359518/1‑91 (MQ=255)
                   aaGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTta                                                                        >  1:1134059/1‑91 (MQ=255)
                     gCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTtata                                                                      <  1:3378610/91‑1 (MQ=255)
                        aTTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACt                                                                   <  2:2854917/91‑1 (MQ=255)
                         ttCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTa                                                                  <  2:3281742/91‑1 (MQ=255)
                          tCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaa                                                                 >  2:3679735/1‑91 (MQ=255)
                           cTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaaa                                                                <  1:3056915/91‑1 (MQ=255)
                           cTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTaaa                                                                <  2:2919558/91‑1 (MQ=255)
                            ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                               <  1:3754082/91‑1 (MQ=255)
                            ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                               >  1:3031916/1‑91 (MQ=255)
                            ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                               >  2:1024747/1‑91 (MQ=255)
                            ttGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAAt                                                               >  2:1217060/1‑91 (MQ=255)
                             tGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATc                                                              <  2:1869458/91‑1 (MQ=255)
                              gACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCg                                                             >  1:3312586/1‑91 (MQ=255)
                                   ccGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTg                                                        <  1:1498387/91‑1 (MQ=255)
                                    cGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTga                                                       >  1:2324081/1‑91 (MQ=255)
                                     gAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgat                                                      <  1:918592/91‑1 (MQ=255)
                                        aaTGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTgataga                                                   <  2:621788/91‑1 (MQ=255)
                                             cTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGc                                              <  1:1473196/91‑1 (MQ=255)
                                              tATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCt                                             <  2:1669817/91‑1 (MQ=255)
                                              tATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCt                                             >  2:1870105/1‑91 (MQ=255)
                                               aTTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                            <  1:2402131/91‑1 (MQ=255)
                                               aTTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                            <  1:1011823/91‑1 (MQ=255)
                                               aTTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                            >  1:3832776/1‑91 (MQ=255)
                                               aTTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                            >  2:3134501/1‑91 (MQ=255)
                                               aTTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCtt                                            >  2:2608011/1‑91 (MQ=255)
                                                ttGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCttt                                           >  2:18426/1‑91 (MQ=255)
                                                  gATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGt                                         >  2:2365291/1‑91 (MQ=255)
                                                     aCCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                      <  1:74845/91‑1 (MQ=255)
                                                     aCCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                      <  1:824385/91‑1 (MQ=255)
                                                     aCCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                      <  1:2542620/91‑1 (MQ=255)
                                                     aCCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                      <  1:2079731/91‑1 (MQ=255)
                                                     aCCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACt                                      <  1:1964012/91‑1 (MQ=255)
                                                       cATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACtta                                    >  1:517442/1‑91 (MQ=255)
                                                        aTTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttat                                   >  2:3227208/1‑91 (MQ=255)
                                                        aTTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttat                                   >  2:2529694/1‑91 (MQ=255)
                                                         ttttGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTCTAACTAAATCGGATTGATAGATAAGCTTTGTACttatt                                  >  1:3387762/1‑91 (MQ=255)
                                                          tttGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttattt                                 >  1:2820427/1‑91 (MQ=255)
                                                            tGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttatttat                               >  2:1809170/1‑91 (MQ=255)
                                                            tGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttatttat                               >  2:2261146/1‑91 (MQ=255)
                                                            tGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACttatttat                               <  2:2444437/91‑1 (MQ=255)
                                                             gATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTtata                              >  2:546454/1‑91 (MQ=255)
                                                                  taaaTATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                         <  2:2813395/89‑1 (MQ=255)
                                                                  tGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                         <  2:885401/91‑1 (MQ=255)
                                                                  tGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                         <  2:347859/91‑1 (MQ=255)
                                                                  tGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                         <  2:1492406/91‑1 (MQ=255)
                                                                  tGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGt                         >  1:2258237/1‑91 (MQ=255)
                                                                      tataACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGa                     <  1:2140592/91‑1 (MQ=255)
                                                                       ataACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAt                    <  1:1825452/91‑1 (MQ=255)
                                                                       ataACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAt                    <  1:2590901/91‑1 (MQ=255)
                                                                        taACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                   <  1:2261146/91‑1 (MQ=255)
                                                                        taACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtt                   <  1:858036/91‑1 (MQ=255)
                                                                          aCAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAtttt                 <  2:14813/91‑1 (MQ=255)
                                                                           cAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttttt                <  2:2051420/91‑1 (MQ=255)
                                                                           cAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGAttttt                <  2:2312594/91‑1 (MQ=255)
                                                                            aaaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAg                          <  2:2029292/80‑1 (MQ=255)
                                                                            aaaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGActtttt               >  1:3616843/1‑91 (MQ=255)
                                                                             aaaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTg              >  1:1300056/1‑91 (MQ=255)
                                                                              aaaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGa             >  2:274464/1‑91 (MQ=255)
                                                                               aaGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAt            >  1:1689565/1‑91 (MQ=255)
                                                                                aGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGAtt           >  2:1533068/1‑91 (MQ=255)
                                                                                 gCCACAGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTg          >  2:734130/1‑91 (MQ=255)
                                                                                   cacaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATCGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTTATTGgt        >  2:480262/1‑91 (MQ=255)
                                                                                   cacaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgt        >  1:388092/1‑91 (MQ=255)
                                                                                   cacaTAACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTATGAACTTATTTATATCCGTTCGATTTTTTGATCGgt        >  2:3507095/1‑91 (MQ=255)
                                                                                   cacaGTAATGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgt        >  1:2487412/1‑91 (MQ=255)
                                                                                    acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg       >  1:41025/1‑91 (MQ=255)
                                                                                    acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg       >  1:3785905/1‑91 (MQ=255)
                                                                                    acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg       <  1:2061111/91‑1 (MQ=255)
                                                                                    acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg       >  1:210120/1‑91 (MQ=255)
                                                                                    acaTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtg       >  1:2599875/1‑91 (MQ=255)
                                                                                     caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt      >  1:3479524/1‑91 (MQ=255)
                                                                                     caTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGgtgt      >  2:2970542/1‑91 (MQ=255)
                                                                                      aTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGGTAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTa     <  1:3080950/91‑1 (MQ=255)
                                                                                      aTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTa     <  1:2808957/91‑1 (MQ=255)
                                                                                      aTTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTa     <  1:1072647/91‑1 (MQ=255)
                                                                                       ttACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaa    <  2:524521/91‑1 (MQ=255)
                                                                                       ttACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTaa    >  1:2783354/1‑91 (MQ=255)
                                                                                          |                                                                                         
AGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTACTTATTTATATCAGTTCGATTTTTTGATTGGTGTAAAA  >  USA300TCH1516_ALE/917240‑917419

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: