Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,603 T→A 15.9% intergenic (+36/‑763) hdeD → / → gadE acid‑resistance membrane protein/gad regulon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,657,6030TA15.9% 236.1 / 28.3 93intergenic (+36/‑763)hdeD/gadEacid‑resistance membrane protein/gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (44/34);  new base A (4/11);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 4.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCT  >  NC_000913/3657504‑3657702
                                                                                                   |                                                                                                   
tGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAt                                                                                                     <  2:155855/100‑1 (MQ=255)
 gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                    <  2:1759980/100‑1 (MQ=255)
 gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                    <  2:566878/100‑1 (MQ=255)
 gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAta                                                                                                    >  1:155857/1‑99 (MQ=255)
  aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                    <  2:1118194/99‑1 (MQ=255)
   aCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                  <  2:1482743/100‑1 (MQ=255)
     tGATATTTAGGGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGAc                                                                                                <  2:2047722/100‑1 (MQ=255)
     tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGAc                                                                                                >  1:1800333/1‑100 (MQ=255)
     tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGAc                                                                                                >  1:1563931/1‑100 (MQ=255)
      gATAATTAGCGCCACCAGCCGGTTCAGCTTCGACAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACGCGGTGCATAGGCACGAGATGTGACt                                                                                               >  2:646383/1‑100 (MQ=255)
        tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                             <  1:2188321/100‑1 (MQ=255)
        tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                             <  2:245511/100‑1 (MQ=255)
        tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                             <  1:1426167/100‑1 (MQ=39)
        tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                             <  2:529012/100‑1 (MQ=39)
         atTTAGCTTCGCCAGCCTATTTAGCTTCGCTTGTTTTTTTGTTAAGCAGCAATAATTATCCCGGTTGTTGCCCGGATCATAGTCACATCAAGTGATTATg                                                                                            <  2:298734/100‑1 (MQ=255)
         atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATg                                                                                            <  1:639450/100‑1 (MQ=39)
         atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATg                                                                                            <  1:553187/100‑1 (MQ=39)
          tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                           <  1:437228/100‑1 (MQ=38)
           ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGAt                                                                                          <  1:687596/100‑1 (MQ=255)
           ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGAt                                                                                          <  2:651387/100‑1 (MQ=255)
             aGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATcc                                                                                        <  1:1946645/100‑1 (MQ=255)
             aGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATcc                                                                                        <  2:2371991/100‑1 (MQ=255)
              gcgcCGCCTGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATcc                                                                                        <  1:1251778/99‑1 (MQ=38)
              gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                       >  2:606432/1‑100 (MQ=255)
              gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCg                                                                                       <  1:1240109/100‑1 (MQ=38)
               cgccgcCAGCCTGTTCCGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGa                                                                                      >  2:698665/1‑100 (MQ=255)
               cgccgcCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCg                                                                                       >  1:2160794/1‑99 (MQ=255)
                gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                       >  2:753183/1‑98 (MQ=255)
                 ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                    <  1:1123195/100‑1 (MQ=255)
                 ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                    <  1:1123181/100‑1 (MQ=255)
                   gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGCTTa                                                                                   >  2:467068/1‑99 (MQ=255)
                      aGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAc                                                                               <  2:572262/100‑1 (MQ=255)
                      aGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATAc                                                                               >  1:1032495/1‑100 (MQ=255)
                       gCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATAct                                                                              >  1:338344/1‑100 (MQ=255)
                        ccTGTTCAGCTTCTCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActc                                                                             <  1:1490962/100‑1 (MQ=255)
                        ccTGTTCAGCTTCACCCGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTCGTCACCCGGACCATAGTCACTTGATGTGACTATGATCCGATTAATAc                                                                               >  2:1148935/1‑98 (MQ=255)
                          tGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctc                                                                           <  2:1645373/100‑1 (MQ=255)
                          tGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActc                                                                             <  2:995730/98‑1 (MQ=255)
                             tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTcc                                                                        <  2:714177/100‑1 (MQ=255)
                                gCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCt                                                                     >  1:297717/1‑100 (MQ=255)
                                gCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCt                                                                     >  2:114351/1‑100 (MQ=255)
                                gCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCt                                                                     >  1:298986/1‑100 (MQ=255)
                                   tCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTAcn                                                                  <  2:830005/100‑2 (MQ=255)
                                   tCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACg                                                                  <  2:338340/100‑1 (MQ=255)
                                     gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGCTACGATGAATACTCTCTCCGCTACGc                                                                 >  2:598617/1‑99 (MQ=255)
                                     gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                >  1:353430/1‑100 (MQ=255)
                                     gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                >  1:1580631/1‑100 (MQ=255)
                                     gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCAGCTACGCa                                                                >  2:576666/1‑100 (MQ=255)
                                     gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTGAATTGATGTGACTATGATCCGATTAATACTCTCTCTGCTACg                                                                  >  2:1030630/1‑98 (MQ=255)
                                      ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAg                                                               >  1:878101/1‑100 (MQ=255)
                                      ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAg                                                               >  1:2495822/1‑100 (MQ=255)
                                         gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                            <  1:1167753/100‑1 (MQ=255)
                                          tttGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgtt                                                           <  2:316672/100‑1 (MQ=255)
                                            tGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                         <  1:1112770/100‑1 (MQ=255)
                                             gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                        <  1:395121/100‑1 (MQ=255)
                                             gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                        <  2:468981/100‑1 (MQ=255)
                                              ttCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAg                                                       >  2:2149483/1‑100 (MQ=255)
                                              ttCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAg                                                       >  2:1901299/1‑100 (MQ=255)
                                               tCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGa                                                      >  2:1019429/1‑100 (MQ=255)
                                                cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                     >  1:2242351/1‑100 (MQ=255)
                                                cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                     >  2:1542298/1‑100 (MQ=255)
                                                cGTTAAGCAGCAATAATTACCCCGGTTGTCACACGGATAATAGTCACTTGATGTGACTATGATCCGATTAATACTATCTCCGCTACGAAGTGTTGTAGAt                                                     >  2:1542269/1‑100 (MQ=255)
                                                 gTTAAGCAGCAATGATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAGCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                                    <  2:68441/100‑1 (MQ=255)
                                                 gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                                    >  1:2048126/1‑100 (MQ=255)
                                                 gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                                    >  1:377648/1‑100 (MQ=255)
                                                 gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                                    <  2:323252/100‑1 (MQ=255)
                                                 gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAGCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                                    >  2:1010493/1‑100 (MQ=255)
                                                     agcagcAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                >  2:1696877/1‑100 (MQ=255)
                                                     agcagcAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                >  2:779610/1‑100 (MQ=255)
                                                     agcagcAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                >  2:328040/1‑100 (MQ=255)
                                                     accagcaATAATTACCCCGGTTGTCACCCGGGTCATAGTCACATGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                >  1:1734050/3‑100 (MQ=255)
                                                         gcaATACTTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGAGACTATTATCCGACTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgagc                                            >  2:444516/1‑100 (MQ=255)
                                                         gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTTTGGTAGATCAGGTgcgc                                            >  1:2041626/1‑100 (MQ=255)
                                                          caATAATTACACCGGTTGTCACCCGGATCATAGTCACTTGATGTGAGTATGATCCGTTTAATACTCTCTCCGCTACGAAGTTTTGTAGAGCAATTGCGCa                                           >  2:2459578/1‑100 (MQ=255)
                                                               atTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATc                                      >  1:2384334/1‑100 (MQ=255)
                                                               atTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATc                                      >  1:520799/1‑100 (MQ=255)
                                                                   ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCACACTATCATTg                                  <  2:894502/100‑1 (MQ=255)
                                                                    cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGa                                 <  1:2047820/100‑1 (MQ=255)
                                                                     gcGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTATGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                <  2:364829/99‑1 (MQ=255)
                                                                      cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaaa                               <  1:101809/100‑1 (MQ=255)
                                                                             cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                        <  1:1435245/100‑1 (MQ=255)
                                                                              aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                       <  1:767798/100‑1 (MQ=255)
                                                                              aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                       <  1:435183/100‑1 (MQ=255)
                                                                               cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCTGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                      >  2:1810809/1‑100 (MQ=255)
                                                                                 cGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGCAGATCAATTGCGCACTATCATTGAAATAATTACCTGc                    >  1:495381/1‑100 (MQ=37)
                                                                                  ggATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                   <  1:1667395/100‑1 (MQ=255)
                                                                                  ggAGCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATAGCGCACTATCATTGAAATAATTACCTGCt                   >  1:1652673/1‑100 (MQ=38)
                                                                                   gATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                  >  1:681570/1‑100 (MQ=38)
                                                                                     tCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGt                <  2:2495807/100‑1 (MQ=255)
                                                                                      cATAGTCGCTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg               >  1:1340392/1‑100 (MQ=255)
                                                                                       aTAGTCACATCAAGTGACTATCATCCGATTAATACTCTCTCCACTAAGCATTGCTGTAGATCAATTGCGCACCATCACTGAAATAATTAGCTGCTACTGa              >  2:627950/1‑100 (MQ=255)
                                                                                        tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat             <  2:1259238/100‑1 (MQ=255)
                                                                                           tctcTTGATGTGTCTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCTCACTATCATTGAAATAATTACCTGCTAGTGattat          <  2:2178269/97‑1 (MQ=255)
                                                                                           tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGattat          <  2:1832140/100‑1 (MQ=255)
                                                                                             aCTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTAttt        >  2:2426630/1‑100 (MQ=255)
                                                                                                tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCaa     >  1:1496272/1‑100 (MQ=255)
                                                                                                tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCaa     >  2:1508931/1‑100 (MQ=255)
                                                                                                 gATGTGACTATGATCCGATTTATACTCTCTCCGCTACGCAGTGTTGTAGTTCAAGTGCGCACTATCGTTGAAATAATTACCTGCTAGTGATTATTTCAAc    <  2:1669308/100‑1 (MQ=255)
                                                                                                  aTGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGTTCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAcc   <  2:955868/100‑1 (MQ=255)
                                                                                                   tgtgACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCt  >  2:57961/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCT  >  NC_000913/3657504‑3657702

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: