Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,357:1 +T 17.3% intergenic (+151/+21) ydbL → / ← feaR DUF1318 family protein/transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,3571.T17.3% 304.2 / 68.4 98intergenic (+151/+21)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base . (39/42);  new base T (12/5);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

Only 100 of 101 total aligned reads displayed.

TTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTC  >  NC_000913/1446259‑1446455
                                                                                                   |                                                                                                  
tttCTTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                   <  1:2902380/100‑1 (MQ=255)
tttCCTGAAACAAGGTGAATATGCAAAAACTCCTGTCAAATTGCCTTGTGCCCTGGAAAATGCATAGTCAATTATCTGTTTTAACTAATTGGCGTTGCAg                                                                                                   <  2:1068199/100‑1 (MQ=255)
  tCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGTCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTa                                                                                                 <  2:1915523/100‑1 (MQ=255)
  ccgatctaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTa                                                                                                 <  1:1105848/93‑1 (MQ=255)
   ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTAc                                                                                                >  2:127935/1‑100 (MQ=255)
   ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTAc                                                                                                <  2:506456/100‑1 (MQ=255)
    cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                               >  2:2274850/1‑100 (MQ=255)
    cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACa                                                                                               >  2:2657986/1‑100 (MQ=255)
      gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCatgtactg                                                                                              <  2:1434653/100‑8 (MQ=255)
      gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATg                                                                                             <  1:2366278/100‑1 (MQ=255)
       aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCTCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGc                                                                                            <  1:195224/100‑1 (MQ=255)
         aCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTaagatcgg                                                                                          >  2:1105839/1‑93 (MQ=255)
            aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGc                                                                                       <  1:370003/100‑1 (MQ=38)
             ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGcc                                                                                      <  1:2870668/100‑1 (MQ=38)
             ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                      >  1:986798/1‑100 (MQ=255)
             ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                      >  2:2533520/1‑100 (MQ=255)
             ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGAGTTGCAGTACATGCAACGc                                                                                       >  1:508835/1‑99 (MQ=255)
               tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTTCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGcc                                                                                      <  1:2971739/98‑1 (MQ=255)
               tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                    <  2:2771722/100‑1 (MQ=255)
                gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGTGCTGAAAAATGCATAGTCAAATATCGGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAt                                                                                   <  1:1019089/100‑1 (MQ=255)
                 aaaatacAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCaa                                                                                    <  2:1616566/92‑1 (MQ=37)
                 aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAtt                                                                                  <  1:91037/100‑1 (MQ=255)
                   tatTCAAAATCTCCTGTCAAATTGTCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAt                                                                                   <  1:2551335/97‑1 (MQ=255)
                   tatTCAAAAACTCCTGTCAAATTACATAAAAAAAAAAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttag                                                                                <  2:1028722/100‑1 (MQ=25)
                        aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGc                                                                           >  2:225339/1‑100 (MQ=255)
                        aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGc                                                                           >  1:1525350/1‑100 (MQ=255)
                         aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCg                                                                          <  2:3039404/100‑1 (MQ=255)
                             cTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAt                                                                      <  1:1382777/100‑1 (MQ=255)
                             cTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAt                                                                      <  1:1762711/100‑1 (MQ=255)
                               ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                    <  2:2275565/100‑1 (MQ=255)
                                 tGTCAGATTTCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCTACGCTAATTAGTTAGCGGAATTTa                                                                   <  2:2018027/99‑1 (MQ=255)
                                  gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                 >  1:1060354/1‑100 (MQ=255)
                                  gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                 >  2:678298/1‑100 (MQ=255)
                                   tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGt                                                                <  2:1339658/100‑1 (MQ=255)
                                     aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                              >  2:960147/1‑100 (MQ=255)
                                        ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                           >  1:2403014/1‑100 (MQ=255)
                                        ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                           >  1:2403255/1‑100 (MQ=255)
                                              tttGCCCTGAAATATGCATAGTCAACTATCTGTTTTAACTAATTGGCCTTGTAATACATGCAACATCAATTAGTTAGCTTAATTTACGTCGATAGTCGcc                                                     <  2:174552/100‑1 (MQ=255)
                                              tttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGcc                                                     >  1:1448441/1‑100 (MQ=255)
                                                 gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAATGCCAATGAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                  >  2:2133905/1‑100 (MQ=255)
                                                  ctcTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                  <  1:655075/97‑1 (MQ=255)
                                                  cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGc                                                 >  1:546199/1‑100 (MQ=255)
                                                  cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGc                                                 >  1:2804623/1‑100 (MQ=255)
                                                   ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGc                                                 >  2:279596/1‑99 (MQ=255)
                                                   ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                >  2:2255290/1‑100 (MQ=255)
                                                   ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  2:2402986/100‑1 (MQ=255)
                                                   ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  2:1005663/100‑1 (MQ=255)
                                                   ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  1:1608710/100‑1 (MQ=255)
                                                   ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  2:2403227/100‑1 (MQ=255)
                                                    cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                               >  1:171294/1‑100 (MQ=255)
                                                     tGAAAAATGCTTATTCAAATTTTTGTTTTAACTATTTGGCGTTGCAGTACATGCAACGCCAATTAGTGAGCGGAATTTACGTCGATACTCGCCTGGCGt                                               <  1:1974280/99‑1 (MQ=38)
                                                     tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                              >  2:2008007/1‑100 (MQ=255)
                                                     tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                              >  2:1297026/1‑100 (MQ=255)
                                                        aaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTATTTAGCGGAAATTACGACGATACACGCCTGGCGTCa                                             >  2:1440512/1‑98 (MQ=39)
                                                          aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATccc                                         <  1:1549084/100‑1 (MQ=255)
                                                          aaTGCAGAGTCAGATATCTGTTTTAACTAATTCGCGTTGCAGTACAGGCAACTCCAATTAGTTATCGGAATTTACGTCTATTCTCGCCTGGCGtcctcc                                          <  2:2641491/99‑1 (MQ=255)
                                                              cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                     >  2:1374596/1‑100 (MQ=255)
                                                              cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCAAAATTTACGACAATACTCACCTGGCGTCATCCCAAAg                                     >  1:889475/1‑100 (MQ=38)
                                                               aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                    <  2:1211266/100‑1 (MQ=255)
                                                                tAGTCAAATATCTGTTTTAACTAATTGACGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGCTACTCGCCTGGCGTCATCCCAAAGCg                                   >  2:659741/1‑100 (MQ=255)
                                                                 aGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGt                                  <  1:795950/100‑1 (MQ=255)
                                                                  gTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                 >  1:2356389/1‑100 (MQ=255)
                                                                  gTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGAGGAATTTACGTAGAAAATCGCCTGGCGGCATCCAAAAGCGtt                                 >  1:2475091/1‑100 (MQ=38)
                                                                  gTCAAATATCTGTTTGAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                 <  1:2993774/100‑1 (MQ=255)
                                                                   tCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCg                                   >  2:1937827/1‑97 (MQ=255)
                                                                    cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               <  1:788045/100‑1 (MQ=255)
                                                                    cAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               <  1:127936/100‑1 (MQ=255)
                                                                     aaaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCt                              >  1:780681/1‑100 (MQ=255)
                                                                     aaaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCt                              <  1:2214267/100‑1 (MQ=255)
                                                                     aaaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCt                              <  2:1593730/100‑1 (MQ=255)
                                                                      aaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCAAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               <  1:1138863/98‑1 (MQ=255)
                                                                       atatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTa                            <  1:166984/100‑1 (MQ=255)
                                                                       atatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTa                            <  1:1118876/100‑1 (MQ=255)
                                                                             gTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                      <  1:2188399/100‑1 (MQ=255)
                                                                               tttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGt                    <  1:2529077/100‑1 (MQ=255)
                                                                                  cgatctttGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                 <  1:1984547/94‑1 (MQ=37)
                                                                                   aCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                >  1:372620/1‑100 (MQ=255)
                                                                                      aaTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa             >  1:1055325/1‑100 (MQ=38)
                                                                                      aaTTGGCGTTGCAGTACATGCAACGCCAATAAGATAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAGCGTTGaaaaa             >  2:1997907/1‑100 (MQ=255)
                                                                                       aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  1:541342/1‑100 (MQ=38)
                                                                                       aTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            <  1:2617999/100‑1 (MQ=255)
                                                                                        ttGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCCTCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg           >  2:2892689/1‑100 (MQ=38)
                                                                                        ttGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaagatcg           >  2:1984524/1‑95 (MQ=37)
                                                                                        ttGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg           >  2:2881036/1‑100 (MQ=38)
                                                                                        ctGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg           <  1:258831/99‑1 (MQ=38)
                                                                                         tGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa          >  2:1684699/1‑100 (MQ=38)
                                                                                          ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGag         >  2:258828/1‑99 (MQ=38)
                                                                                          ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGACGTCATCACAAAGAGTTGCTTAAATACCGTTGAAAAATGAc         >  1:190397/1‑100 (MQ=37)
                                                                                          ggCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTAGCCTAGCGTCATCCCAAAAAGTTGCTTAAATACAGTTGAAAAATTAc         >  1:659733/1‑100 (MQ=25)
                                                                                          ggCGTTGCATGTACTGCAACGCCAATTAGATAGCGGAATTTACGTCAATACTCGCCTGACGTAATACCAAAGCGTTGCTTAAATACCGTTGAAAAATGa          >  2:1261770/1‑99 (MQ=25)
                                                                                          ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAACTCCCGTTGAAAAATGAc         >  2:600868/1‑100 (MQ=255)
                                                                                          ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  1:2781564/1‑100 (MQ=255)
                                                                                           gCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct        >  2:1553514/1‑100 (MQ=38)
                                                                                            cGTTGCAGTACATGCAACGCAAATTAGTTAGCGGAATTTACGTCGATTCTCGCCTGGCGTCATACCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  2:3025855/1‑98 (MQ=255)
                                                                                             gTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActct      >  2:1544217/1‑100 (MQ=255)
                                                                                              ttGCAGTACATGCAACGACAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc       <  1:2296762/98‑1 (MQ=255)
                                                                                                gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt   <  1:1544235/100‑1 (MQ=255)
                                                                                                gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt   >  1:1499837/1‑100 (MQ=255)
                                                                                                gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt   >  1:877703/1‑100 (MQ=255)
                                                                                                 cAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc  <  2:2458342/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTC  >  NC_000913/1446259‑1446455

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