Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,682,492 G→A A319V (GCC→GTC)  dctA ← C4‑dicarboxylic acid, orotate and citrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,682,4920GA96.4% 352.4 / ‑0.5 110A319V (GCC→GTC) dctAC4‑dicarboxylic acid, orotate and citrate transporter
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (50/56);  total (53/57)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

AAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGAT  >  NC_000913/3682395‑3682587
                                                                                                 |                                                                                               
aaaGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGAcc                                                                                               <  2:1100982/100‑1 (MQ=255)
 aaGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCt                                                                                              <  2:2682235/100‑1 (MQ=255)
  aGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTg                                                                                             <  2:1109073/100‑1 (MQ=255)
   gCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGc                                                                                                >  2:2414854/1‑96 (MQ=255)
   gCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTgg                                                                                            >  2:2238898/1‑100 (MQ=255)
    ccACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTggg                                                                                           >  1:1501840/1‑100 (MQ=255)
     cACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGc                                                                                          <  2:1705118/100‑1 (MQ=255)
     cACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGc                                                                                          >  1:2440464/1‑100 (MQ=255)
       cgatctgTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGa                                                                                        <  1:601427/94‑1 (MQ=255)
        tACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGAt                                                                                       >  1:2329694/1‑100 (MQ=255)
         aCCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATa                                                                                      <  1:963027/100‑1 (MQ=255)
          cccGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATaa                                                                                     <  1:2344521/100‑1 (MQ=255)
            cGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAac                                                                                   <  1:2833943/100‑1 (MQ=255)
            cGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAac                                                                                   <  1:51602/100‑1 (MQ=255)
             gTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGaagatcg                                                                                  >  2:601419/1‑94 (MQ=255)
               tACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACAcc                                                                                <  2:1979347/100‑1 (MQ=255)
               tACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACAcc                                                                                >  2:1483962/1‑100 (MQ=255)
                aCCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACAccg                                                                               >  2:2012729/1‑100 (MQ=255)
                 ccccTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACAccgc                                                                              <  1:2666221/100‑1 (MQ=255)
                 ccccTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACAccgc                                                                              >  1:1390787/1‑100 (MQ=255)
                 ccccTGCCGCCCCATTAGAAGAAAGCAGCAAAACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACAccgc                                                                              >  1:2054812/1‑100 (MQ=255)
                    cTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACcgccgc                                                                           >  2:2052829/1‑100 (MQ=255)
                    cTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACcgccgc                                                                           >  2:2052850/1‑100 (MQ=255)
                      gccgccCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCa                                                                         >  1:2297551/1‑100 (MQ=255)
                      gccgccCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCa                                                                         >  1:1759409/1‑100 (MQ=255)
                       ccgccCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCAt                                                                        <  2:291385/100‑1 (MQ=255)
                        cgccCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCAtt                                                                       >  1:2517726/1‑100 (MQ=255)
                            ccTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTc                                                                    >  2:2597694/1‑99 (MQ=255)
                            ccTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCa                                                                   >  2:1045861/1‑100 (MQ=255)
                            ccTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCa                                                                   >  2:1556190/1‑100 (MQ=255)
                             cTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAg                                                                  <  2:281991/100‑1 (MQ=255)
                              tttAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAgg                                                                 <  2:1076125/100‑1 (MQ=255)
                               ttAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGt                                                                <  2:1508295/100‑1 (MQ=255)
                                tAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGta                                                               <  2:366839/100‑1 (MQ=255)
                                  gaagaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGtata                                                             >  1:1306504/1‑100 (MQ=255)
                                  gaagaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGtata                                                             >  2:1134634/1‑100 (MQ=255)
                                  gaagaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGtata                                                             >  2:263865/1‑100 (MQ=255)
                                  gaagaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGtata                                                             <  2:352453/100‑1 (MQ=255)
                                  gaagaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGtata                                                             >  2:1107699/1‑100 (MQ=255)
                                   aagaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATGGTCAGGTatat                                                            >  2:2187672/1‑100 (MQ=255)
                                    agaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATc                                                           >  1:1926988/1‑100 (MQ=255)
                                    agaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACATGGGCGATAAACACCGCCGCCATTGTCAGGTATATc                                                           >  2:2699605/1‑100 (MQ=255)
                                    agaaAGCAGCAACAAGATTAACAGCGTGACTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGACGATAAACACCGCCGCCATTGTCAGGTATATc                                                           <  2:1071888/100‑1 (MQ=255)
                                     gaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCg                                                          <  1:2492605/100‑1 (MQ=255)
                                     gaaAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCg                                                          >  1:1270971/1‑100 (MQ=255)
                                            gcaACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGcc                                                   >  1:478900/1‑100 (MQ=255)
                                            gcaACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGcc                                                   >  1:873293/1‑100 (MQ=255)
                                             caACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCAATTGTCAGGTATATCGATGTGCCa                                                  >  1:1015372/1‑100 (MQ=255)
                                                cacGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCa                                               >  2:822930/1‑100 (MQ=255)
                                                   gATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAgg                                            <  1:2207385/100‑1 (MQ=255)
                                                    aTTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGt                                           <  2:2494493/100‑1 (MQ=255)
                                                      tAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTa                                         <  2:1166893/100‑1 (MQ=255)
                                                       aaCAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTaa                                        >  2:2546575/1‑100 (MQ=255)
                                                       aaCAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTaa                                        >  2:912388/1‑100 (MQ=255)
                                                       aaCAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTaa                                        >  1:1270051/1‑100 (MQ=255)
                                                        aCAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTaaa                                       >  1:921135/1‑100 (MQ=255)
                                                         cttccgatcttGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAAc                                      <  1:2394717/91‑1 (MQ=255)
                                                          aGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACg                                     >  2:1755556/1‑100 (MQ=255)
                                                          aGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACg                                     <  2:2795014/100‑1 (MQ=255)
                                                          aGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACg                                     <  1:1483975/100‑1 (MQ=255)
                                                          aGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTCACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAAc                                      <  2:2140370/99‑1 (MQ=255)
                                                           gCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGa                                    >  2:2175384/1‑100 (MQ=255)
                                                                aTTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGc                               <  2:503817/100‑1 (MQ=255)
                                                                aTTTGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGc                               <  2:833472/100‑1 (MQ=255)
                                                                 tttGGTGGACGATATCCATCTGTCTGTTAGTGACCTGCGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGcc                              <  2:2054786/100‑1 (MQ=255)
                                                                 tttGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGcc                              <  1:1245739/100‑1 (MQ=255)
                                                                 tttGGTGGACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGcc                              <  2:2294113/100‑1 (MQ=255)
                                                                   tggtggACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAAcagatcggaag                            >  2:2394686/1‑90 (MQ=255)
                                                                    ggtggACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGt                           <  2:2297522/100‑1 (MQ=255)
                                                                    ggtggACGATATCCATCTGACTGTGAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGt                           <  2:512873/100‑1 (MQ=255)
                                                                      tggACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCg                         >  1:129359/1‑100 (MQ=255)
                                                                      tggACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCg                         >  2:2291474/1‑100 (MQ=255)
                                                                      tggACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCg                         >  2:2291739/1‑100 (MQ=255)
                                                                      tggACGATATCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCg                         <  2:2329664/100‑1 (MQ=255)
                                                                         aCGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGAAGCCTGTCGGGa                      >  2:1233073/1‑100 (MQ=255)
                                                                             tatCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATgtc                  <  1:2052855/100‑3 (MQ=255)
                                                                             tatCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGAc                  <  1:2052876/100‑1 (MQ=255)
                                                                              atCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGAcc                 <  1:1816119/100‑1 (MQ=255)
                                                                              atCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGAcc                 >  2:1729724/1‑100 (MQ=255)
                                                                              atCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGAcc                 <  2:1390774/100‑1 (MQ=255)
                                                                               tCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCa                <  2:1501824/100‑1 (MQ=255)
                                                                               tCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCa                <  1:3045329/100‑1 (MQ=255)
                                                                               tCCATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCa                <  2:1832531/100‑1 (MQ=255)
                                                                                 cATCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACAAGTAGCCTGTCGGGATGACCAg               >  2:611016/1‑99 (MQ=255)
                                                                                  aTCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGcc             >  2:1025139/1‑100 (MQ=255)
                                                                                  aTCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGcc             <  2:2985187/100‑1 (MQ=255)
                                                                                  aTCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGcc             <  1:1867378/100‑1 (MQ=255)
                                                                                  aTCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGcc             <  1:706990/100‑1 (MQ=255)
                                                                                   tCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGc              <  2:1857462/98‑1 (MQ=255)
                                                                                   tCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGccc            <  1:2012754/100‑1 (MQ=255)
                                                                                   tCTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCAAGTCGGGATGACCAGc              >  2:2161849/1‑98 (MQ=255)
                                                                                    cTGACTGTTAGTGACCTGGGCGATAAACACCGCGGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTACACGAGTAGCCTGTCGGGATGACCAGcccc           >  2:1832149/1‑100 (MQ=255)
                                                                                    cTGACTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGcccc           >  1:506732/1‑100 (MQ=255)
                                                                                     tGACTGTGAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGcccga          <  1:2072621/100‑3 (MQ=255)
                                                                                        cTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCccacca       >  1:1301993/1‑100 (MQ=255)
                                                                                        cTGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCccacca       >  2:2037987/1‑100 (MQ=255)
                                                                                        cTGTAAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGGTAAACGAGTAGCCTGTCGGGATGACCAGCCccacca       >  2:874963/1‑100 (MQ=255)
                                                                                         tGTTAGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTAGGCATGACCAGCCccac         >  2:2938634/1‑97 (MQ=255)
                                                                                             aGTGACCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGAt  >  1:2752011/1‑100 (MQ=255)
                                                                                             aGTGACCTGGGCGAACAACACCACCGCCATTGTCATGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCTTATCGCGATGACCAGCCCCaccacc     <  2:736671/97‑1 (MQ=255)
                                                                                                 |                                                                                               
AAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATCAAGGTTAAACGAGTAGCCTGTCGGGATGACCAGCCCCACCACCGAT  >  NC_000913/3682395‑3682587

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