Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,791,137 A→G V70A (GTA→GCA)  tdh ← L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,791,1370AG98.6% 729.1 / ‑3.8 209V70A (GTA→GCA) tdhL‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (98/108);  total (100/109)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

ACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTACCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAGTCCCGC  >  NC_000913/3791040‑3791233
                                                                                                 |                                                                                                
accaccACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAAtacctgcc                                                                                                <  1:2988603/100‑1 (MQ=255)
 ccaccaCGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGACCGCCGATCTTGCAGCCTTTCACTTCCTGACCAATACCTGCCa                                                                                               >  2:695166/1‑100 (MQ=255)
   accacGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACt                                                                                             >  2:5378026/1‑100 (MQ=255)
    ccacGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCAcag                                                                                            >  2:1459385/1‑98 (MQ=255)
      acGACAGTTGCGGAAATGACCGCAGGTGATATGGCCTTCGCCAGAAACGCGATTGCGGAGCTTGAAGAATTTCACTTCCTGATCAATACCTGCAACTTCa                                                                                          >  2:3894252/1‑100 (MQ=255)
        gACAGTTGCGGCAATGACCACAGGTGTTATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGGCTTTCACTTCCTGACCAATTCCTGCCACTGCAcg                                                                                        <  2:1894542/100‑2 (MQ=255)
          cAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGa                                                                                      >  1:1602040/1‑100 (MQ=255)
           aGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGAc                                                                                     <  1:1951535/100‑1 (MQ=255)
            gTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACa                                                                                    >  1:881879/1‑100 (MQ=255)
            gTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACa                                                                                    <  1:4364774/100‑1 (MQ=255)
             ttGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACat                                                                                   <  2:2664598/100‑1 (MQ=255)
             ttGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACat                                                                                   <  2:3814912/100‑1 (MQ=255)
                cGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATAtt                                                                                <  1:5383528/100‑1 (MQ=255)
                 ggCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTc                                                                               >  2:2569948/1‑100 (MQ=255)
                    aaTGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATg                                                                            >  1:1659619/1‑100 (MQ=255)
                    aaTGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATg                                                                            >  2:1211147/1‑100 (MQ=255)
                     aTGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATgg                                                                           <  1:2489620/100‑1 (MQ=255)
                     aTGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTAACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATgg                                                                           >  1:5427780/1‑100 (MQ=255)
                      tGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGc                                                                          <  2:2158246/100‑1 (MQ=255)
                      tGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGc                                                                          <  2:4819176/100‑1 (MQ=255)
                      tGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGc                                                                          <  1:3278005/100‑1 (MQ=255)
                       gACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGcc                                                                         >  1:4487779/1‑100 (MQ=255)
                       gACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGcc                                                                         >  2:740573/1‑100 (MQ=255)
                        aCCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGccc                                                                        >  2:3924476/1‑100 (MQ=255)
                        aCCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATCGccc                                                                        >  2:3591341/1‑100 (MQ=255)
                         ccACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCa                                                                       >  1:4998520/1‑100 (MQ=255)
                            caGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGa                                                                    >  1:1986049/1‑100 (MQ=255)
                              ggTGATATGGCCTTCGCCAGAAAGGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACTTGCGACTTCACCGACATATTAATGGCCCACGAcc                                                                  >  2:1240683/1‑100 (MQ=255)
                               gTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCa                                                                 >  1:5028167/1‑100 (MQ=255)
                                tGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCAt                                                                <  2:5164061/100‑1 (MQ=255)
                                 gATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATc                                                               <  2:1582550/100‑1 (MQ=255)
                                  atatGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCg                                                              <  2:2464050/100‑1 (MQ=255)
                                  atatGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCg                                                              <  1:2255143/100‑1 (MQ=255)
                                    ctGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGc                                                            <  1:214778/99‑1 (MQ=255)
                                      ggCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCAc                                                          >  1:1255150/1‑100 (MQ=255)
                                       gCCTTCGCCGGAAACCCGATCGCCGTGCTTGACGCCTTTGACGTCCGTACCAGTCCCTGCTACTTCGCCGACATAGTCATGGCGTACGACCATCGGCAcc                                                         <  2:3244218/100‑1 (MQ=255)
                                       gCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGACTTTCACTTCCTGACCAAGACCTGCCAGTTCACCGACATGTTCATGGCCCACGACCATCGGCAcc                                                         >  1:5174155/1‑100 (MQ=255)
                                         cTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCgg                                                       <  1:5093259/100‑1 (MQ=255)
                                         cTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCgg                                                       <  1:5276665/100‑1 (MQ=255)
                                         cTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCgg                                                       >  1:2879964/1‑100 (MQ=255)
                                            cGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGAt                                                    >  1:109963/1‑100 (MQ=255)
                                             gCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATTGGCACCGGGATg                                                   >  1:1850473/1‑100 (MQ=255)
                                              ccAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCTCTTCGCCGACATATTCATGGCCCACGACCGTCGGCACCGGGATgg                                                  <  1:4835028/100‑1 (MQ=255)
                                              ccAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATgg                                                  >  1:5532279/1‑100 (MQ=255)
                                               cAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGt                                                 <  1:2595071/100‑1 (MQ=255)
                                               cAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGt                                                 >  1:964734/1‑100 (MQ=255)
                                                 gAAACGCGATCGCCGATCTTGAAGCCTTTCACTCCCTGACCAATACCTACCACTTCACCCACATATTCATGGCCCACGACCATCGGCACCGGGAGGGatt                                               >  2:4555295/1‑100 (MQ=255)
                                                  aaaCGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGtttt                                              <  1:4803568/100‑1 (MQ=255)
                                                    aCGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGtttttt                                            >  1:2634061/1‑100 (MQ=255)
                                                      gcgATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGc                                          >  1:1188644/1‑100 (MQ=255)
                                                      gcgATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGc                                          >  1:1314711/1‑100 (MQ=255)
                                                       cgATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCg                                         <  1:35084/100‑1 (MQ=255)
                                                       cgATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCg                                         <  1:2910001/100‑1 (MQ=255)
                                                       cgATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCg                                         <  1:5109551/100‑1 (MQ=255)
                                                        gATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGa                                        >  1:3424544/1‑100 (MQ=255)
                                                        gATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGG‑TTTTTGCGAc                                       <  1:1477297/100‑1 (MQ=255)
                                                         aTCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGAc                                       >  1:2412914/1‑100 (MQ=255)
                                                         aTCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGAc                                       <  1:2694187/100‑1 (MQ=255)
                                                          tCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGAcc                                      <  2:127539/100‑1 (MQ=255)
                                                          tCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGAcc                                      <  2:2837433/100‑1 (MQ=255)
                                                           cGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCa                                     <  1:59394/100‑1 (MQ=255)
                                                           cGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCa                                     >  2:24569/1‑100 (MQ=255)
                                                           cGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCa                                     >  2:2748640/1‑100 (MQ=255)
                                                            gCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCAc                                    >  2:5476307/1‑100 (MQ=255)
                                                            gCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCAc                                    >  1:3919328/1‑100 (MQ=255)
                                                             ccGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACt                                   >  2:4303045/1‑100 (MQ=255)
                                                             ccGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACt                                   <  2:4870232/100‑1 (MQ=255)
                                                              cGATCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTc                                  >  1:67087/1‑100 (MQ=255)
                                                                 tCTTGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATc                               <  2:1697292/100‑1 (MQ=255)
                                                                    tGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCa                            >  1:748537/1‑100 (MQ=255)
                                                                    tGAAGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCa                            >  1:1237765/1‑100 (MQ=255)
                                                                      aaGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGt                          <  2:4471130/100‑1 (MQ=255)
                                                                      aaGCCTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGt                          <  2:3726024/100‑1 (MQ=255)
                                                                          cTTTCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTtata                      >  2:4277095/1‑100 (MQ=255)
                                                                           tttCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCCTCCCAGTTATAg                     >  1:4622888/1‑100 (MQ=255)
                                                                           tttCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAg                     <  1:3779066/100‑1 (MQ=255)
                                                                            ttCACTTCCTGACCAATACCTACCACTTCACCGACATATTCATGGCCCACGACCATCGGCACTGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGa                    >  2:3999544/1‑100 (MQ=255)
                                                                             tCACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGAt                   <  1:3258194/100‑1 (MQ=255)
                                                                              cACTTCCTGACCAATCCCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGAtg                  <  2:3069829/100‑1 (MQ=255)
                                                                              cACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGAtg                  <  2:3422802/100‑1 (MQ=255)
                                                                              cACTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGAtg                  <  2:2008675/100‑1 (MQ=255)
                                                                                cTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGAtgtg                <  2:3334202/100‑1 (MQ=255)
                                                                                cTTCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGAtgtg                >  2:4620371/1‑100 (MQ=255)
                                                                                  tCCTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTTaa              >  2:2140143/1‑100 (MQ=255)
                                                                                   ccTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAAc             <  1:4316805/100‑1 (MQ=255)
                                                                                    cTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACg            >  1:3084153/1‑100 (MQ=255)
                                                                                    cTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACg            >  2:760689/1‑100 (MQ=255)
                                                                                    cTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACg            >  1:4667337/1‑100 (MQ=255)
                                                                                    cTGACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACg            <  2:4825085/100‑1 (MQ=255)
                                                                                      gACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACGACTCATCCCAGTTATAGATGTGAACGTc          <  2:931844/100‑1 (MQ=255)
                                                                                      gACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTc          >  2:4813760/1‑100 (MQ=255)
                                                                                      gACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTc          <  1:3617005/100‑1 (MQ=255)
                                                                                      gACCAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTc          <  1:98159/100‑1 (MQ=255)
                                                                                        ccAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAg        <  1:4011994/100‑1 (MQ=255)
                                                                                        ccAATACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAg        >  2:1310055/1‑100 (MQ=255)
                                                                                           aTACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTTCGACCACTCATCCCAGTTATAGATGTGAACGTCAGTcc     >  2:4407377/1‑100 (MQ=255)
                                                                                           aTACCTGCCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAGTcc     <  2:4036351/100‑1 (MQ=255)
                                                                                             acctgccACTTCACCGACGTATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAGTCCCg   <  2:1789106/100‑1 (MQ=255)
                                                                                              cctgccACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAGTCCCGc  >  1:789778/1‑100 (MQ=255)
                                                                                              cctgccACTTCACCGACATATTCATGGCACACGACCATCGGCACCGGGATGGTTTTTTGCGACCATTCATACCAGTTATAGATGTGAACGTCAGTCCCGc  >  2:3230098/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTACCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAGTCCCGC  >  NC_000913/3791040‑3791233

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: ‑