Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,513,930 A→T 100% L413F (TTA→TTT)  USA300HOU_RS07375 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,513,9300AT99.6% 777.5 / ‑4.9 217L413F (TTA→TTT) USA300HOU_RS07375hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (111/104);  minor base C (1/0);  total (113/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 217 total aligned reads displayed.

CTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCG  >  USA300TCH1516_ALE/1513843‑1514016
                                                                                       |                                                                                      
ctctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTAtt                                                                                     >  2:1075993/1‑91 (MQ=255)
ctctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTAtt                                                                                     >  1:714973/1‑91 (MQ=255)
   ttGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTa                                                                                  <  1:3845397/91‑1 (MQ=255)
   ttGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTa                                                                                  >  1:2812593/1‑91 (MQ=255)
      gATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTa                                                                               <  1:2912816/91‑1 (MQ=255)
      gATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTa                                                                               >  1:2177328/1‑91 (MQ=255)
       aTAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAg                                                                              >  1:3349058/1‑91 (MQ=255)
        tAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAGAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAgg                                                                             <  1:3138163/91‑1 (MQ=255)
         aaaGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGa                                                                            <  1:856412/91‑1 (MQ=255)
          aaGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGat                                                                           >  1:3197743/1‑91 (MQ=255)
          aaGCGTACAAAAAATACTGCAATAAAAGTAACAAAAAAAAAAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGat                                                                           <  2:3925260/91‑1 (MQ=255)
             cGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGattat                                                                        <  2:2832158/91‑1 (MQ=255)
               tACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGAttatta                                                                      >  2:2443696/1‑91 (MQ=255)
                aCAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTaa                                                                     >  1:1450763/1‑91 (MQ=255)
                 cAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAAt                                                                    <  1:1897052/91‑1 (MQ=255)
                 cAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAAt                                                                    <  1:2686279/91‑1 (MQ=255)
                 aaaaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAAt                                                                    <  2:1988313/90‑1 (MQ=255)
                  aaaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTGGGATTATTAATg                                                                   <  1:3053584/91‑1 (MQ=255)
                   aaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGc                                                                  <  2:1884842/91‑1 (MQ=255)
                   aaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTGATTGTATTAGGATTATTAATGc                                                                  >  2:3189061/1‑91 (MQ=255)
                    aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGcc                                                                 >  2:2763712/1‑91 (MQ=255)
                     aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCa                                                                <  2:682540/91‑1 (MQ=255)
                     aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCa                                                                >  1:1159653/1‑91 (MQ=255)
                      aaTACTGCAATAAAAGTACCACATCGCGTAAAACAACCAGCATCTATATCACTCCTAATCTGTTTTATTGTATTAGGCTTATTAATGccct                                                               >  1:1254335/1‑89 (MQ=255)
                      aaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCat                                                               <  2:655819/91‑1 (MQ=255)
                        tACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatta                                                             <  1:626450/91‑1 (MQ=255)
                         aCTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                            >  2:1728414/1‑91 (MQ=255)
                         aCTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                            >  2:400557/1‑91 (MQ=255)
                         aCTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                            >  1:3233609/1‑91 (MQ=255)
                         aCTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                            >  2:1851779/1‑91 (MQ=255)
                          cTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                           >  1:151868/1‑91 (MQ=255)
                          cTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                           >  2:2057610/1‑91 (MQ=255)
                          cTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                           >  1:3119799/1‑91 (MQ=255)
                           tGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTAttt                                                          <  1:1313428/91‑1 (MQ=255)
                             cAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTgg                                                        <  1:2958294/91‑1 (MQ=255)
                              aaTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGa                                                       >  2:45898/1‑91 (MQ=255)
                                 aaaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                    <  1:1586602/91‑1 (MQ=255)
                                 aaaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                    <  2:2355923/91‑1 (MQ=255)
                                  aaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATACTATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTAt                                                   <  2:3786570/91‑1 (MQ=255)
                                  aaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTAt                                                   <  1:3963389/91‑1 (MQ=255)
                                   aaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATc                                                  <  2:1020280/91‑1 (MQ=255)
                                    aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                 >  1:3368828/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                               >  2:1286655/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                               >  2:2666072/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                               <  2:3508373/91‑1 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                               >  1:1801998/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                               >  2:2324204/1‑91 (MQ=255)
                                        ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                             >  1:2607060/1‑91 (MQ=255)
                                        ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                             >  1:2545325/1‑91 (MQ=255)
                                        ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                             >  1:243383/1‑91 (MQ=255)
                                        ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                             >  1:1691126/1‑91 (MQ=255)
                                        ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAAGATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                             >  2:3259647/1‑91 (MQ=255)
                                          acaTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATc                                           <  1:100349/91‑1 (MQ=255)
                                          acaTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATc                                           <  1:2422811/91‑1 (MQ=255)
                                            aTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTAACCt                                         >  1:3550348/1‑91 (MQ=255)
                                             tCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCtt                                        >  2:565195/1‑91 (MQ=255)
                                               gcgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTTGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGt                                      >  2:801758/1‑91 (MQ=255)
                                               gcgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGt                                      <  1:2463172/91‑1 (MQ=255)
                                                cgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTa                                     <  1:791525/91‑1 (MQ=255)
                                                  tAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAtt                                   <  1:492202/91‑1 (MQ=255)
                                                   aaaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAttt                                  >  2:3427558/1‑91 (MQ=255)
                                                   aaaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAttt                                  <  1:2151340/91‑1 (MQ=255)
                                                    aaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTta                                 >  1:2680456/1‑91 (MQ=255)
                                                      aaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTtata                               >  1:3693943/1‑91 (MQ=255)
                                                        aaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAt                             <  2:2226732/91‑1 (MQ=255)
                                                        aaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAt                             <  2:3583733/91‑1 (MQ=255)
                                                         aCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                            <  1:1674742/91‑1 (MQ=255)
                                                          ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                           >  1:2937659/1‑91 (MQ=255)
                                                          ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                           >  1:1120799/1‑91 (MQ=255)
                                                            aGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGaa                         <  2:1369497/91‑1 (MQ=255)
                                                             gCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATCATTGAAt                        >  2:1104862/1‑91 (MQ=255)
                                                             gCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAt                        <  1:1114687/91‑1 (MQ=255)
                                                             gCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAt                        <  1:3041169/91‑1 (MQ=255)
                                                             gCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAt                        >  1:3436181/1‑91 (MQ=255)
                                                                tCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTaa                     >  1:3017098/1‑91 (MQ=255)
                                                                 cTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAat                    >  2:1840328/1‑91 (MQ=255)
                                                                 cTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAat                    <  1:456399/91‑1 (MQ=255)
                                                                  tataTCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAata                   >  2:848980/1‑91 (MQ=255)
                                                                  tataTCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAata                   >  2:2697207/1‑91 (MQ=255)
                                                                   atatCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAatat                  >  1:2328837/1‑91 (MQ=255)
                                                                   atatCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAatat                  <  1:3125154/91‑1 (MQ=255)
                                                                    tatCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatt                 >  2:1283243/1‑91 (MQ=255)
                                                                     atCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatta                <  1:1355637/91‑1 (MQ=255)
                                                                     atCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATCTGGATTATCACTTATCCTGATATTTATAATTGAATTAAtatta                >  2:3670852/1‑91 (MQ=255)
                                                                     atCACTCATAATATGTCTTATTGTACTAGGATTATTAATGCCATTATTTGGATTATCCCTTATCCTTGTATTTATAATTGAATTAAtatta                >  2:2850459/1‑91 (MQ=255)
                                                                      tCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtattat               >  1:3857561/1‑91 (MQ=255)
                                                                      tCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtattat               >  2:2123437/1‑91 (MQ=255)
                                                                         cTCATAGTGTGTTTTATTGTCTTAGGATTAGTGATTCCATTATTTGGATTAACACTTAGCCTTGTATTTATAATTGAATTAATATtatata            <  2:1591601/91‑1 (MQ=255)
                                                                          tCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTatatat           >  1:901664/1‑91 (MQ=255)
                                                                            ataataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTa         >  2:1459646/1‑91 (MQ=255)
                                                                            ataataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTa         <  2:2775035/91‑1 (MQ=255)
                                                                            ataataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTa         <  2:3206851/91‑1 (MQ=255)
                                                                            ataataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTa         >  2:3023967/1‑91 (MQ=255)
                                                                             taataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaa        <  2:1805255/91‑1 (MQ=255)
                                                                             taataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaa        >  2:1389294/1‑91 (MQ=255)
                                                                               atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAg      <  1:3071463/91‑1 (MQ=255)
                                                                                tatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGa     <  1:2093607/91‑1 (MQ=255)
                                                                                 atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAt    >  2:3131130/1‑91 (MQ=255)
                                                                                 atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAt    <  2:2900289/91‑1 (MQ=255)
                                                                                   gTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcg  >  2:1167618/1‑91 (MQ=255)
                                                                                       |                                                                                      
CTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCG  >  USA300TCH1516_ALE/1513843‑1514016

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