Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,358,482 T→A 18.5% I46F (ATC→TTC)  cadA ← lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,358,4820TA18.5% 314.8 / 61.5 140I46F (ATC→TTC) cadAlysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base T (52/62);  new base A (15/11);  total (67/73)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

Only 100 of 143 total aligned reads displayed.

ACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTAAAA  >  NC_000913/4358383‑4358579
                                                                                                   |                                                                                                 
aCGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATTCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAt                                                                                                   >  2:2651079/1‑100 (MQ=255)
aCGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAt                                                                                                   >  1:59382/1‑100 (MQ=255)
   gCAGGTTCTCGTTGATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCTAGTCAAAAATAACGCCGCATAGACGCGCATTGTTTTCGATCAg                                                                                                <  2:1332478/100‑1 (MQ=255)
      ggTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGttt                                                                                             >  1:2374756/1‑100 (MQ=255)
        ttCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTta                                                                                           >  2:2753714/1‑100 (MQ=255)
         tctcGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGGCGCGCATTGTTTTCGATCAGTTTtaa                                                                                          <  1:1847103/100‑1 (MQ=255)
          ctcGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTtaat                                                                                         <  1:121469/100‑1 (MQ=255)
           tcGTTCATTTTGCTAATTTCTTCGCACAGCTTGATATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTtaata                                                                                        <  1:2380533/100‑1 (MQ=255)
           tcGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTtaata                                                                                        <  1:1048819/100‑1 (MQ=255)
           tcGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTtaata                                                                                        <  1:948517/100‑1 (MQ=255)
            cGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTtaataa                                                                                       >  1:2340990/1‑100 (MQ=255)
             gTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAg                                                                                      >  1:2845492/1‑100 (MQ=255)
             gTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCAAGTCAAAAATAACGCCGAACAGACGCGCATTGTTTTCGATCAGTTTTAATAAg                                                                                      >  2:563500/1‑100 (MQ=255)
              ttCATTTTGCTAATTTCTTCGAACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAgt                                                                                     <  2:2585522/100‑1 (MQ=255)
               tCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAgtc                                                                                    <  2:1246664/100‑1 (MQ=255)
                cATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCCCAGACGCGCATTGTTTTCGATCAGTTTTAATAAgtcg                                                                                   <  2:1614120/100‑1 (MQ=255)
                cATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAgtcg                                                                                   >  1:2763820/1‑100 (MQ=255)
                cATTTTGCTAATTCCTTCGCACAGCTCGAGACTATATTTATCGCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAgtcg                                                                                   >  1:2048184/1‑100 (MQ=255)
                   tttGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCa                                                                                <  1:865507/100‑1 (MQ=255)
                   tttGCTAATTTCTTCACACAGCTCGAGATTATCTCTATCCCAGTAAAAAATAACGTCGAACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGGCGTCa                                                                                <  2:951112/100‑1 (MQ=255)
                    ttGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCAc                                                                               >  2:3004146/1‑100 (MQ=255)
                     tGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACg                                                                              >  2:1618688/1‑100 (MQ=255)
                       cTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGt                                                                            <  1:2856671/100‑1 (MQ=255)
                       cTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGt                                                                            >  1:2439779/1‑100 (MQ=255)
                       cTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGt                                                                            >  2:1845321/1‑100 (MQ=255)
                             tcttcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCg                                                                      <  1:1370483/100‑1 (MQ=255)
                               ttcGCACAGCTCGAGATTATATTTATCCTAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCggg                                                                    <  1:1202636/100‑1 (MQ=255)
                               ttcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCTTTGTTTTCGATCAGTTTTAATAAGTCGTTACGGTCGTTCggg                                                                    <  1:2049074/100‑1 (MQ=255)
                               ttcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCggg                                                                    <  1:2714695/100‑1 (MQ=255)
                                tcGCACAGCTCTAGATTATATTTATCCCAGTCAAAAATATCGCCGTACAGACGCGCATTGTTTTCGATCAGTGTTAATAAGTCGTCACGGTCGTTCGGGt                                                                   <  1:2649959/100‑1 (MQ=255)
                                tcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGt                                                                   <  1:742367/100‑1 (MQ=255)
                                tagCACAGCTCGAGATTATATTTATCCTAGTCAAAAATAACGCCGCACAGACGCGCATAGTTTTCGATCAGTTTTAATAAGTAGTCACGGTCGTTCGGGt                                                                   >  1:2736515/3‑100 (MQ=255)
                                   cacaGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTaaa                                                                <  2:514478/100‑1 (MQ=255)
                                   cacaGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTaaa                                                                >  1:920350/1‑100 (MQ=255)
                                   cacaGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTaaa                                                                >  2:754217/1‑100 (MQ=255)
                                   cacaGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGCTTACGATCATTTTTAATAAGTCGTCACGGTCGTTCGGGTaaa                                                                >  2:2131984/1‑100 (MQ=255)
                                        cTCGAGATTATGTTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATc                                                           >  2:2175401/1‑100 (MQ=255)
                                        cTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATCAATCGTCACGGTCGTTCGGGTAAACAATc                                                           >  2:2116247/1‑100 (MQ=255)
                                        cTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATc                                                           >  1:1947967/1‑100 (MQ=255)
                                        cTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATc                                                           >  1:166004/1‑100 (MQ=255)
                                         tCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCt                                                          >  2:355388/1‑100 (MQ=255)
                                          cGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTg                                                         <  1:1696864/100‑1 (MQ=255)
                                          cGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACACACGCGCATTGTTTTCGATCAATTTTAATAAGTCGTCACGGACATTCGGGTAAACAATCt                                                          >  1:840532/1‑99 (MQ=255)
                                            cgatcttatTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGa                                                       <  1:2252037/94‑1 (MQ=255)
                                            agaTTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGa                                                       <  2:594585/100‑1 (MQ=255)
                                            agaTTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGa                                                       <  2:1652341/100‑1 (MQ=255)
                                            agaTTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATATgca                                                       >  2:1928399/1‑98 (MQ=255)
                                              aTTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAg                                                     >  2:85786/1‑100 (MQ=255)
                                               ttATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGTTCGTACGAGTAAACAATATGGAAGt                                                    >  1:1170146/1‑100 (MQ=255)
                                               ttATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGt                                                    <  2:1347136/100‑1 (MQ=255)
                                                 atatTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTc                                                  >  2:1469900/1‑100 (MQ=255)
                                                  tatTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAgatcg                                                 >  2:2252013/1‑96 (MQ=255)
                                                  tatTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCa                                                 <  2:2277262/100‑1 (MQ=255)
                                                   atTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAATAATCTGGAAGTTCAg                                                <  2:618113/100‑1 (MQ=255)
                                                   atTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCAGGTAAAAAATCTAGAATTTCa                                                 >  1:1488230/1‑99 (MQ=255)
                                                   atTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGATCGTTCGGGTAAACAATCTGGAAGTTCAg                                                >  1:2420911/1‑100 (MQ=255)
                                                     ttATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGAc                                              <  1:2613857/100‑1 (MQ=255)
                                                     ttATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGAc                                              >  2:1275951/1‑100 (MQ=255)
                                                        tCCCAGTCAAAAATAACGCCGCACAGACGCGCATTTTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGt                                            <  1:2132007/99‑1 (MQ=255)
                                                        tCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGtt                                           >  2:1019485/1‑100 (MQ=255)
                                                         cccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTc                                          >  1:899637/1‑100 (MQ=255)
                                                         cccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTc                                          >  1:339695/1‑100 (MQ=255)
                                                          ccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCa                                         <  1:2326891/100‑1 (MQ=255)
                                                          ccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCa                                         >  2:2057603/1‑100 (MQ=255)
                                                            attCAAAAATAACGACGTATAGACGCGTATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAg                                       <  1:752953/98‑1 (MQ=255)
                                                             gTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAgc                                      >  1:2558815/1‑100 (MQ=255)
                                                               cAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAgc                                      >  1:500924/1‑98 (MQ=255)
                                                                 aaaaTAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAg                                       <  2:1510270/95‑1 (MQ=255)
                                                                  aaaTAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGcgcgcg                                 <  2:420246/100‑1 (MQ=255)
                                                                  aaaTAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGcgcgcg                                 <  2:2095989/100‑1 (MQ=255)
                                                                   aaTAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGa                                <  2:2472263/100‑1 (MQ=255)
                                                                     tAACGCCGCACATACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGAt                               <  1:1092069/99‑1 (MQ=255)
                                                                     tAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATg                              <  1:1143666/100‑1 (MQ=255)
                                                                       aCGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGaa                            >  1:1491459/1‑100 (MQ=255)
                                                                        cgccgcACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAg                           >  1:2735684/1‑100 (MQ=255)
                                                                         gccgcACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGATCAGACGTTCAAGCGCGCGATGAAGt                          <  1:1845342/100‑1 (MQ=255)
                                                                           cgcACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTc                        <  2:2288329/100‑1 (MQ=255)
                                                                              acaGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACg                     <  1:570419/100‑1 (MQ=255)
                                                                              acaGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACg                     <  1:1795983/100‑1 (MQ=255)
                                                                                  aCGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATg                 <  2:1694676/100‑1 (MQ=255)
                                                                                   cgcgCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATgg                >  1:1608352/1‑100 (MQ=255)
                                                                                    gcgcATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATggg               <  2:2827009/100‑1 (MQ=255)
                                                                                     cgcATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGCTCAAGCGCGCGATGAAGTTCACGGATggg               >  2:2232540/1‑99 (MQ=255)
                                                                                      gcATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGtt             >  2:1947583/1‑100 (MQ=255)
                                                                                        aTTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGttc            <  2:1037692/99‑1 (MQ=255)
                                                                                        aTTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGttct           <  2:1409086/100‑1 (MQ=255)
                                                                                          ttttttCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGttcttc         <  2:198976/98‑1 (MQ=255)
                                                                                          tGTTTTCGATCAGTTTTAATATGTCGTCAGGGTCGTTCGGGTAAACAATCTTGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGttcttc         <  2:323024/100‑1 (MQ=255)
                                                                                          tGTTTTCGATCAGTTTTAATACGTCGTCACGGTCGTTCGGGTAAACAGTCTGGAAGTTGTGAGGTTCAAGCGCGCGATGAGGTTAACGTATGGGttcttc         <  2:2468225/100‑1 (MQ=255)
                                                                                          tGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCTCGATGAAGTTCACGGATGGGttcttc         <  2:2744987/100‑1 (MQ=255)
                                                                                           gTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGttcttc         <  2:2606331/99‑1 (MQ=255)
                                                                                           gTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTtcttct        <  1:1896397/100‑1 (MQ=255)
                                                                                           gTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTtcttct        >  1:100634/1‑100 (MQ=255)
                                                                                           gTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTtcttct        >  1:1958929/1‑100 (MQ=255)
                                                                                            ttttCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTcttctt       <  1:490015/100‑1 (MQ=255)
                                                                                            ttttCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTcttctt       <  2:22904/100‑1 (MQ=255)
                                                                                             tttCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTcttctt       <  2:1982206/99‑1 (MQ=255)
                                                                                              ttCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAGGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTa     >  1:217198/1‑100 (MQ=255)
                                                                                               tCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTaa    <  2:92292/100‑1 (MQ=255)
                                                                                                 gATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTaaaa  >  1:2411775/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
ACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTAAAA  >  NC_000913/4358383‑4358579

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: