Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,362 Δ1 bp 17.5% intergenic (+156/+16) 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,3620A.17.5% 287.9 / 58.4 103intergenic (+156/+16)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base A (44/40);  new base . (10/8);  total (54/49)
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.93e-01

Only 100 of 105 total aligned reads displayed.

TGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAA  >  NC_000913/1446264‑1446461
                                                                                                   |                                                                                                   
tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACan                                                                                                    <  1:2328079/100‑2 (MQ=255)
tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACAt                                                                                                    <  2:3175660/100‑1 (MQ=255)
tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTATGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACAt                                                                                                    >  1:2587718/1‑100 (MQ=255)
  aaaCAAGGTGAATATTCAAAATCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACAt                                                                                                    <  1:3105909/98‑1 (MQ=255)
    aCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCaa                                                                                                <  1:1035547/100‑1 (MQ=255)
     cAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAAc                                                                                               <  2:25771/100‑1 (MQ=39)
     cAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAAc                                                                                               >  2:352090/1‑100 (MQ=255)
      aaGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACg                                                                                              <  1:59030/100‑1 (MQ=38)
       aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGc                                                                                             <  1:2345498/100‑1 (MQ=38)
        ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                            <  1:2191964/100‑1 (MQ=38)
        ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGcc                                                                                            <  2:3679079/100‑1 (MQ=255)
        ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTACCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                            <  2:2659699/100‑1 (MQ=38)
         gTGAATATTCAAAAACTCCTGTCATATTGCATTTTGCCCTGATAAATTCATATTCATATATCTTTTTTAACTACTTGGCGTTGCA‑GTACATGCatc                                                                                               <  2:3090401/96‑3 (MQ=255)
           gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAt                                                                                         <  1:2640291/100‑1 (MQ=255)
           gAATATAAAAAAATTCTTGTCAAATTGCTTTTTGACCTAAAAAATGCATCCTCAAATATCTCTTTAAACTAATTGGCGTTGCA‑GTACTTGCAACGCCa                                                                                           <  2:1088691/98‑1 (MQ=255)
            aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtt                                                                                        <  2:34420/100‑1 (MQ=255)
                ttCAAAAACTCCTGTCAAATTGCCTTTTTCCCTGAAAAATTCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGGAC‑TGCAACGCCAAttagt                                                                                     <  2:3790533/99‑1 (MQ=25)
                ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                    >  2:2565789/1‑100 (MQ=255)
                ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                    >  2:3362742/1‑100 (MQ=255)
                  cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                                  >  1:3657348/1‑100 (MQ=255)
                   aaaagCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttagttaa                                                                                  <  2:1341723/99‑2 (MQ=25)
                   aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGc                                                                                 >  1:151965/1‑100 (MQ=255)
                    aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCg                                                                                <  1:3198414/100‑1 (MQ=255)
                     aaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCCGTTTTAACTAATTGGCGTTGCA‑ATACATGCAAAGCCCATTAATTAGCgg                                                                               >  2:2159595/1‑100 (MQ=255)
                          ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                          <  2:399221/100‑1 (MQ=255)
                              tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGt                                                                      >  1:1835554/1‑100 (MQ=255)
                               cAAATTTCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                     <  2:2937453/100‑1 (MQ=255)
                               cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                     >  1:260704/1‑100 (MQ=255)
                                 aaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGa                                                                   >  1:2564904/1‑100 (MQ=255)
                                  aTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGAt                                                                  <  2:1182944/100‑1 (MQ=255)
                                   ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                                 >  2:742685/1‑100 (MQ=255)
                                        ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                            >  2:456132/1‑100 (MQ=255)
                                        ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                            <  1:3362768/100‑1 (MQ=255)
                                        ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGc                                                            <  1:1906781/100‑1 (MQ=255)
                                         tttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGcc                                                           >  1:2814346/1‑100 (MQ=255)
                                            gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCt                                                          >  2:3380917/1‑98 (MQ=255)
                                             cccTGAAAAATGCATAGTAAAATATCTGTTTTAACTAATTGGCGTAGCA‑GTACATGCAACGCCTATTAGTTAGCGGAATTTACGTCGAGACTCTcc                                                           >  2:1265809/1‑96 (MQ=255)
                                              ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                      <  2:1954369/100‑1 (MQ=255)
                                              ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                      <  1:2307984/100‑1 (MQ=255)
                                              ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                      >  1:1909312/1‑100 (MQ=255)
                                              ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                      <  1:1692292/100‑1 (MQ=255)
                                              ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                      >  2:1768541/1‑100 (MQ=255)
                                                tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                    >  2:102822/1‑100 (MQ=255)
                                                 gAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                    >  2:2418760/1‑99 (MQ=255)
                                                   aaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                                 <  1:2549519/100‑1 (MQ=255)
                                                   aaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTAGCA‑GTACATGCAACGCCAATCAGTTAGCGGAATTTACGTCGATACTCGCCTGGGGt                                                     >  2:1880771/1‑96 (MQ=255)
                                                     aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATccc                                               <  2:679305/100‑1 (MQ=255)
                                                         cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                           >  1:2899238/1‑100 (MQ=255)
                                                          aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                          >  2:1485099/1‑100 (MQ=255)
                                                          aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                          <  1:2061105/100‑1 (MQ=255)
                                                           tAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTATCGGCATTTACGTCTATACTCCCCTGTCGTTATCCCAAAgtg                                         >  2:1630087/1‑98 (MQ=37)
                                                           tAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCg                                         <  1:59482/100‑1 (MQ=255)
                                                             gTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                       <  1:3065087/100‑1 (MQ=255)
                                                              tCAAATATCTGTTTTGACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                      <  1:3764253/100‑1 (MQ=255)
                                                               cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                     >  1:3703894/1‑100 (MQ=255)
                                                               cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                     <  2:1567191/100‑1 (MQ=255)
                                                               cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACAGGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                     <  1:1270057/100‑1 (MQ=255)
                                                                 aaTATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                                   <  2:2564887/100‑1 (MQ=255)
                                                                   tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATAAAACGCCTATTAGTTAGCAGTATTTACATCGATACCCCCCAGGCGTCAACCCAAAGCGt                                        >  2:1687317/1‑93 (MQ=255)
                                                                     tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGAAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                               <  1:2565806/100‑1 (MQ=255)
                                                                       tGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATa                              >  2:479484/1‑99 (MQ=255)
                                                                        gTTTTAACTGATGGGCGTTGCA‑GTAAATGCAACGCCCATGAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGGCATCCCAAAGCGTTGCTTAAATa                              <  2:3414192/98‑1 (MQ=255)
                                                                        gTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                            >  1:3441897/1‑100 (MQ=255)
                                                                           ttAACTAATTGGCTTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                         <  1:3065446/100‑1 (MQ=255)
                                                                           ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                         >  2:3404551/1‑100 (MQ=255)
                                                                            tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                        >  1:2191767/1‑100 (MQ=255)
                                                                             aaCTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                       >  2:2380562/1‑100 (MQ=255)
                                                                                tAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaa                    >  1:3654457/1‑100 (MQ=255)
                                                                                tAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaa                    <  2:2530711/100‑1 (MQ=255)
                                                                                tAATTGGCGTTGCA‑GTACATGAAACCCCAATTAGTTCGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                     <  2:2517715/99‑1 (MQ=255)
                                                                                 aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                   >  1:306018/1‑100 (MQ=255)
                                                                                 aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                   >  2:3079679/1‑100 (MQ=255)
                                                                                 aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                   <  1:2879625/100‑1 (MQ=255)
                                                                                  aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  1:1885326/1‑100 (MQ=38)
                                                                                  aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  2:69705/1‑100 (MQ=38)
                                                                                  aTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  2:2025366/1‑100 (MQ=255)
                                                                                  aTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                  >  2:1958267/1‑100 (MQ=255)
                                                                                  aTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGACGATACTCGCCTGGCGTCATCCCAAAGCGTTGATTACATACAGTTGAAAAAt                  >  1:414771/1‑100 (MQ=255)
                                                                                   ttGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                 <  2:3098103/100‑1 (MQ=255)
                                                                                   ttGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg                 <  1:2173268/100‑1 (MQ=255)
                                                                                    tGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                >  1:2090322/1‑100 (MQ=38)
                                                                                     ggCGTTGCATGTAC‑TGCAACGCCTATTAGTTAGCGGAATTTACTTCGATACTCGCCTTGCGTCATCCCAAAGCTTTGCTTAAATACCGATGAAAAATGAc               >  2:1669320/1‑100 (MQ=37)
                                                                                     ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTATATACCGTTGAAAAATGAc               >  2:2163036/1‑100 (MQ=38)
                                                                                     ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:3733127/1‑100 (MQ=38)
                                                                                     ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  2:1760565/1‑100 (MQ=38)
                                                                                     ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:1347847/1‑100 (MQ=38)
                                                                                     ggCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:2437951/1‑100 (MQ=255)
                                                                                     ggCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:728621/1‑100 (MQ=255)
                                                                                     ggCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc               >  1:2436925/1‑100 (MQ=255)
                                                                                      gCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGCAATTTACGTCGATACTCGCCTGGCGTCATCCCAAACCGTTGCTTAAATACCGTTGAAAAATGAc               <  2:1270824/99‑1 (MQ=255)
                                                                                          tGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg          <  2:2845547/100‑1 (MQ=255)
                                                                                          tGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg          >  1:3161749/1‑100 (MQ=255)
                                                                                           gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt         <  1:3404577/100‑1 (MQ=255)
                                                                                            cA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc        <  1:3351679/100‑1 (MQ=255)
                                                                                               gTACATGCAACGCCAATTAGTTAGCGGAATTTATGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAg      >  2:640665/1‑100 (MQ=255)
                                                                                                tACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGa     <  1:3489611/100‑1 (MQ=255)
                                                                                                tACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGa     <  1:352092/100‑1 (MQ=255)
                                                                                                tACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGa     >  1:2734190/1‑100 (MQ=255)
                                                                                                 aCATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaa    <  2:1842279/100‑1 (MQ=255)
                                                                                                   aTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGATTAAAAAAAAAAAAAAAATGACTCTGGTCAGaaaa  <  2:2372949/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAA  >  NC_000913/1446264‑1446461

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