Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 694,071 G→T 13.7% A50E (GCA→GAA)  USA300HOU_RS03355 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE694,0710GT13.7% 165.5 / 6.5 73A50E (GCA→GAA) USA300HOU_RS03355hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base T (9/2);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTATCGGTTCATCACCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGATAATTTAATTGATTTGAACCATCAAAAAAAGGTGTAAAC  >  USA300TCH1516_ALE/693934‑694211
                                                                                                                                         |                                                                                                                                            
ccTAACCGTTCATCACCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCa                                                                                                                                           <  2:1075944/141‑1 (MQ=255)
   tgcGGTTCATCACCTCAATAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCAttt                                                                                                                                        <  1:3136/139‑1 (MQ=255)
      ggTTCATCACCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTcac                                                                                                                                     >  1:54801/1‑141 (MQ=255)
      ggTTCATCACCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACAACTTATGGGTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTcac                                                                                                                                     >  1:568269/1‑141 (MQ=255)
         tcatcaCCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACaaa                                                                                                                                  >  1:812178/1‑141 (MQ=255)
         tcatcaCCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACaaa                                                                                                                                  >  2:802697/1‑141 (MQ=255)
         tcatcaCCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACaaa                                                                                                                                  >  1:83464/1‑141 (MQ=255)
           atcaCCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATCACCATCCTTTTTTGCCATTTCACAAACt                                                                                                                                >  1:638927/1‑141 (MQ=255)
                cTCAAAAAATTCTAAATCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTaa                                                                                                                           >  1:889603/1‑141 (MQ=255)
                      aaaTTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACg                                                                                                                     <  2:695413/141‑1 (MQ=255)
                        aTTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTCTAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATTTAGGaa                                                                                                                   >  2:31145/1‑141 (MQ=255)
                        aTTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTTTTTTTTTTTAACTTCTAATATACGaa                                                                                                                   >  2:241107/1‑141 (MQ=255)
                        aTTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGaa                                                                                                                   >  1:992968/1‑141 (MQ=255)
                          tCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGACTTTTTACAAACTTTTAATAAActtaaa                                                                                                                 >  2:1201727/1‑136 (MQ=255)
                              aGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATg                                                                                                             <  1:1053211/141‑1 (MQ=255)
                              aGTCTTCAACTGAAATATTATACATCGTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATg                                                                                                             <  1:641275/141‑1 (MQ=255)
                               gtcataaaTTGAAATATGATACATCTTACAAATAAATGCATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGATAATATGc                                                                                                            <  2:143419/135‑1 (MQ=255)
                                  ttCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGCTAATCATTCTTTTTTGCCATTTCACAAACTTCTCATATACTAAAATATGCaaa                                                                                                         >  1:435462/1‑141 (MQ=255)
                                     aaTTCAAATATCATACTGCTAACAAATAAATATATTAGACACACCATATTGTTTTAAACATTTATAATCTTTATAAGTAAATCGATAATCCTTCTCTTTTGACATTTCACAAACTTCTAACATCCGTAAATATGCCAAAtt                                                                                                      <  2:576288/141‑1 (MQ=255)
                                      aTTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATta                                                                                                     <  1:395045/141‑1 (MQ=255)
                                        tGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTTTATTTTTCAAACTTCTAATATACGAAAATATTCAAAATtata                                                                                                   >  2:202555/1‑141 (MQ=255)
                                             tatGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACAttt                                                                                              >  1:896002/1‑141 (MQ=255)
                                             aatGATACATCTTACAAAATAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACAttt                                                                                              <  2:2236/140‑1 (MQ=255)
                                              atGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACAtttt                                                                                             >  2:1172503/1‑141 (MQ=255)
                                               tGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTa                                                                                            >  1:406977/1‑141 (MQ=255)
                                               tGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTa                                                                                            >  1:180601/1‑141 (MQ=255)
                                                 aTACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACa                                                                                          >  1:267555/1‑141 (MQ=255)
                                                 aTACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCAAAAATCTTTCCTTTTAGGTAGTTCACAATCTTCTAATTTTCGAAAATATGCAAAATTATACATTTTACa                                                                                          <  2:1151535/141‑1 (MQ=255)
                                                   cactCTAACAAA‑AAATGCATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAAt                                                                                       <  2:143513/138‑1 (MQ=255)
                                                      tCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTTCATTTTTCTAACTTCTAATATACGAAAATATGCAAAATTATACATTTTATAAAttt                                                                                     >  2:99976/1‑141 (MQ=255)
                                                       cTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAAtttt                                                                                    >  1:1065256/1‑141 (MQ=255)
                                                        ttACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAAttttt                                                                                   <  1:1107050/141‑1 (MQ=255)
                                                         tACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTc                                                                                  >  1:936145/1‑141 (MQ=255)
                                                         tACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTc                                                                                  >  1:961718/1‑141 (MQ=255)
                                                          aCAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTTCTTTTCTCAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCt                                                                                 >  2:54028/1‑141 (MQ=255)
                                                          aCAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCt                                                                                 >  2:159213/1‑141 (MQ=255)
                                                          aCAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCCTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCt                                                                                 >  2:240892/1‑141 (MQ=255)
                                                          aCAAATAAATGTATTAAACACACCATATTGTTTTAAAAATTTACAATCCTTATAAGTAAATCGATAATCATTCTTTTTGGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCt                                                                                 <  2:185067/141‑1 (MQ=255)
                                                               taaatGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTCCATTTCTCAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAg                                                                            >  2:1064909/1‑141 (MQ=255)
                                                                  atGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAATCTTCTAATATACGAAAATATGCAAAATTATACATTTTACATATTTTTATTTTAGaaa                                                                         >  2:1080797/1‑141 (MQ=255)
                                                                  atGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTACACAAACTTCTAGTATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGaaa                                                                         >  2:819635/1‑141 (MQ=255)
                                                                   tGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAac                                                                        >  2:635341/1‑141 (MQ=255)
                                                                       ttAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCa                                                                    <  2:184031/141‑1 (MQ=255)
                                                                       ttAGACACACCATATTGTTTTAAAAATTTACAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCa                                                                    <  2:981208/141‑1 (MQ=255)
                                                                          gACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTTACAAACTAATAATATACGAAAATATCCAAAATTATACAATTTACAAATTTTACTTTTAGAAACACCAtaa                                                                 >  2:611512/1‑139 (MQ=255)
                                                                          gACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTTTATAAACCATGa                                                                 >  2:627052/1‑141 (MQ=255)
                                                                          gACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGa                                                                 >  1:582600/1‑141 (MQ=255)
                                                                          gACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAAAAACACCATGa                                                                 >  1:646296/1‑141 (MQ=255)
                                                                            cacacCATATTTTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGaaa                                                               <  2:304512/141‑1 (MQ=255)
                                                                              cacCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTTTTTTTTTTAAACTTCTAATATACGAAGATATAAAAAATTATACATTTTACAAATTTTTTTTTTAGAATAACCTTGAAAtt                                                             >  2:124936/1‑140 (MQ=255)
                                                                               acCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCt                                                            <  2:568269/141‑1 (MQ=255)
                                                                               acCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCt                                                            <  2:54801/141‑1 (MQ=255)
                                                                                ccATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAACCTTCTACTATACAAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATAtt                                                           >  2:1029949/1‑141 (MQ=255)
                                                                                ccATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGACATCtt                                                           >  1:859958/1‑141 (MQ=255)
                                                                                 cATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGCAATCTTa                                                          >  1:451362/1‑141 (MQ=255)
                                                                                 cATATTGTTTTAAAAATTTATAATCCTTATAAGGAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATAAATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTa                                                          >  1:451386/1‑141 (MQ=255)
                                                                                   tatTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTaaa                                                        >  1:244289/1‑141 (MQ=255)
                                                                                    atTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTTTTTTTTTTAAATTTCTAATATACGAAAATATTAAAAATTATATATTTTATATATTTTTTTTTTAGTACAACCTTGATATTTTaaat                                                       >  2:774354/1‑140 (MQ=255)
                                                                                         tttAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACa                                                  <  1:414491/141‑1 (MQ=255)
                                                                                          ttAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGTCCTTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACaa                                                 <  1:1172503/141‑1 (MQ=255)
                                                                                              aaaTTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTTTTTTTTTTTAACTTCTAAAATACGAAAATATTCAAAATTATACATTTGACAAATTTTTCTTTTAGATACACTCTATAATTTTAAAATTTTAAAttt                                             >  2:206405/1‑141 (MQ=255)
                                                                                                 tttATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTTTTTTTTTTCAAAATTCTAATATACGAAAATATGCTAAATTATACATTTTGTAAATTTTTCTTTTTGTTACCCCATGATATATTAATTTTACATTTTTaat                                          >  2:1006118/1‑140 (MQ=255)
                                                                                                 tttATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAg                                          <  1:635341/141‑1 (MQ=255)
                                                                                                          cTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGAtaatttaa                                 <  1:1067937/141‑1 (MQ=255)
                                                                                                           ttATAAGTAAATCGATAATCCTTCTTTTTTCCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGAtaatttaat                                <  1:1170048/141‑1 (MQ=255)
                                                                                                            tataAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTCAAGAtaatttaatt                               <  1:640475/141‑1 (MQ=255)
                                                                                                            tataAGTAAATCGATAATCACTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATATTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGAtaatttaatt                               <  2:435462/141‑1 (MQ=255)
                                                                                                                aGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAGATTACAAATTTATGATAATTTAATTGAtt                           >  1:809080/1‑141 (MQ=255)
                                                                                                                aGTAAATCGATAATCATTCTTTCTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGATAATTTAATTGAtt                           <  2:365901/141‑1 (MQ=255)
                                                                                                                 gTAAATCGATAATCATTGTTTTTCTCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATGACGAAGTTAAGATAATTTAATTGAttt                          <  2:548293/141‑1 (MQ=255)
                                                                                                                            aTCATTATTTATTCCCATTTCACAAACTTGTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGATAATTTAATTGATTTGAACCATCaaa               <  1:1051011/141‑1 (MQ=255)
                                                                                                                               aTTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGATAATTTAATTGATTTGAACCATCaaaaaa            <  1:202555/141‑1 (MQ=255)
                                                                                                                                     ttgagCCATTCCATAAACTCCTAATATACAAAAATATGCAAAATTATACATTTTACAAATTTTCCTTTTAGAAACACCATGAAAACTGAAAATTACAAATTTAAGATAATTTAATTGATTTGAACCATCAAAAAAAGgtgt      <  1:396684/137‑1 (MQ=255)
                                                                                                                                      tttGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGATAATTTAATTGATTTGAACCATCAAAAAAAGGTGTa     <  1:126710/141‑1 (MQ=255)
                                                                                                                                         gCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGATAATTTAATTGATTTGAACCATCAAAAAAAGGTGTAAAc  >  2:829995/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                            
CCTATCGGTTCATCACCTCAAAAAATTCTAAGTCTTCAATTGAAATATGATACATCTTACAAATAAATGTATTAGACACACCATATTGTTTTAAAAATTTATAATCCTTATAAGTAAATCGATAATCATTCTTTTTTGCCATTTCACAAACTTCTAATATACGAAAATATGCAAAATTATACATTTTACAAATTTTTCTTTTAGAAACACCATGAAATCTTAAAATTACAAATTTAAGATAATTTAATTGATTTGAACCATCAAAAAAAGGTGTAAAC  >  USA300TCH1516_ALE/693934‑694211

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