Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,031 T→A 100% D410V (GAT→GTT)  astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,0310TA96.1% 288.6 / ‑0.1 92D410V (GAT→GTT) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base A (42/46);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGCC  >  NC_000913/1826934‑1827126
                                                                                                 |                                                                                               
ctcGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCg                                                                                               >  2:2386861/1‑100 (MQ=255)
     gcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCCGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgc                                                                                             >  2:2411692/1‑97 (MQ=255)
     gcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcagc                                                                                          <  2:485397/100‑1 (MQ=255)
     gcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGacagc                                                                                          <  2:786338/100‑1 (MQ=255)
       tGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGcagcag                                                                                        <  1:757174/100‑1 (MQ=255)
          aaCGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTaa                                                                                     >  1:182955/1‑100 (MQ=255)
           aCGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAg                                                                                    <  2:2026834/100‑1 (MQ=255)
            cGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAgg                                                                                   <  1:489660/100‑1 (MQ=255)
             ggATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGc                                                                                  <  1:2132393/100‑1 (MQ=255)
              gATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCg                                                                                 <  1:540760/100‑1 (MQ=255)
               aTAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAgg                                                                                   >  2:879289/1‑97 (MQ=255)
               aTAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGa                                                                                <  1:2259451/100‑1 (MQ=255)
                tAAACCGAACCGAGATTCAGTAATTCGCCCAATACATCAAGTCCTCCGCCCCTTCCGCACAACAAGTACGAGTCCGCCAGAACGGCAGCAGCAAAGc                                                                                  >  2:334434/1‑97 (MQ=255)
                 aaaCCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGAAGCAATTGCTGGTCGGCCCGAAGGGCAGCAGTAAGACGATc                                                                              >  2:2248883/1‑100 (MQ=255)
                  aaccgaaccgAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcg                                                                             >  1:51998/1‑100 (MQ=255)
                   accgaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcgc                                                                            <  2:2172447/100‑1 (MQ=255)
                     cgaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCgg                                                                          <  1:1424627/100‑1 (MQ=255)
                      gaaccgaGGTTCAGTGACTGGCGCAATCCATCCGGTGCTGCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgt                                                                         <  2:154313/100‑1 (MQ=255)
                      gaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgt                                                                         <  1:370007/100‑1 (MQ=255)
                        accgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGGCAGAACGGCAGCAGTAAGGCGATCGCGgtag                                                                       <  2:1967218/100‑1 (MQ=255)
                        accgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgtag                                                                       <  1:1425961/100‑1 (MQ=255)
                         ccgaGATTCAGTAATTGGCTCACTACATCCAGTGCTTCCCGCCCTTAGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCACTAAGGCGATcgcg                                                                           >  2:1159654/1‑95 (MQ=255)
                             gaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACg                                                                  >  1:1991872/1‑100 (MQ=255)
                              aTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGa                                                                 <  1:1574499/100‑1 (MQ=255)
                                 cAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATcc                                                              <  1:483749/100‑1 (MQ=255)
                                 cAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATcc                                                              <  1:498560/100‑1 (MQ=255)
                                  aGTAATTGGCTCAATACATCCAGGGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATcc                                                              <  1:50863/99‑1 (MQ=255)
                                   gTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAc                                                            >  1:291242/1‑100 (MQ=255)
                                    tAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAc                                                            >  2:1217324/1‑99 (MQ=255)
                                          gCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCa                                                     <  2:457773/100‑1 (MQ=255)
                                           cTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCAt                                                    >  1:46093/1‑100 (MQ=255)
                                             cAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTg                                                  <  2:1369262/100‑1 (MQ=255)
                                             cAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTg                                                  <  2:844739/100‑1 (MQ=255)
                                               aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAg                                                <  1:868809/100‑1 (MQ=255)
                                               aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAg                                                <  1:2504014/100‑1 (MQ=255)
                                               aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAg                                                <  1:467493/100‑1 (MQ=255)
                                                tACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgc                                               <  1:1941701/100‑1 (MQ=255)
                                                  cATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcgc                                             >  1:458641/1‑100 (MQ=255)
                                                  cATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcgc                                             >  1:1778151/1‑100 (MQ=255)
                                                  cATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcga                                             <  1:2071235/100‑2 (MQ=255)
                                                   aTCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCa                                            <  1:2056837/100‑1 (MQ=255)
                                                    tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                           <  2:458637/100‑1 (MQ=255)
                                                    tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                           >  1:522510/1‑100 (MQ=255)
                                                    tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                           >  2:311017/1‑100 (MQ=255)
                                                     ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                          >  2:70561/1‑100 (MQ=255)
                                                      cAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTa                                         >  2:747681/1‑100 (MQ=255)
                                                      cAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTa                                         <  1:1307697/100‑1 (MQ=255)
                                                       aGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTa                                         >  2:1531136/1‑99 (MQ=255)
                                                         tGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAAc                                      <  2:402179/100‑1 (MQ=255)
                                                         tGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAAc                                      <  1:2536003/100‑1 (MQ=255)
                                                          gCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACa                                     >  2:131751/1‑100 (MQ=255)
                                                            ttCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGc                                   >  2:2363554/1‑100 (MQ=255)
                                                            ttCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGc                                   <  1:980623/100‑1 (MQ=255)
                                                            ttCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGc                                   <  1:2379142/100‑1 (MQ=255)
                                                             tCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCg                                  <  2:1756624/100‑1 (MQ=255)
                                                             tCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCg                                  >  1:95408/1‑100 (MQ=255)
                                                             tCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCg                                  >  1:88413/1‑100 (MQ=255)
                                                              cccgcaaatCGCGAAGCAATTGCGGGTCCGCCAGATAGGCAGCAGTAAGAAGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCg                                  <  1:1660693/91‑1 (MQ=255)
                                                              cccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGt                                 <  1:1248763/100‑1 (MQ=255)
                                                               ccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACAGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGt                                 >  1:1237262/1‑99 (MQ=255)
                                                                cGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTAt                               >  2:1298174/1‑100 (MQ=255)
                                                                    cTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGtt                           <  1:2156860/100‑1 (MQ=255)
                                                                    cTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGtt                           >  1:1533788/1‑100 (MQ=255)
                                                                      tCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                         <  2:937525/100‑1 (MQ=255)
                                                                      tCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                         >  1:8571/1‑100 (MQ=255)
                                                                       cgcgCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcat                        >  1:1793/1‑100 (MQ=255)
                                                                         cgcAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                      <  1:2386876/100‑1 (MQ=255)
                                                                         cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                      >  1:1328899/1‑100 (MQ=255)
                                                                         cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                      <  1:933122/100‑1 (MQ=255)
                                                                           cagcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatcatc                    <  1:85088/100‑1 (MQ=255)
                                                                           cagcaATTGCGGGTCGGCCAGAACGGCAGCAGCAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatcatc                    >  2:71635/1‑100 (MQ=255)
                                                                            agcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatcatca                   <  1:1896846/100‑1 (MQ=255)
                                                                             gcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAc                  >  1:1012359/1‑100 (MQ=255)
                                                                               aaTTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccg                >  2:1919486/1‑100 (MQ=255)
                                                                               aaTTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccg                >  1:1296116/1‑100 (MQ=255)
                                                                                aTTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgc               >  2:1134941/1‑100 (MQ=255)
                                                                                aTTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgc               >  2:1735562/1‑100 (MQ=255)
                                                                                 ttGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgcc              >  2:1359862/1‑100 (MQ=255)
                                                                                  tGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTACCGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgccg             <  2:2487409/100‑1 (MQ=255)
                                                                                  tGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgccg             >  2:352540/1‑100 (MQ=255)
                                                                                  tGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgccg             >  2:1326744/1‑100 (MQ=255)
                                                                                    cGGGTCGACCAGAACGGCAGCAGTAAGTCGATCGCTGTATTAACTATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCAGCACCGCCgg            <  2:1213740/99‑1 (MQ=255)
                                                                                     gggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAt          <  1:333557/100‑1 (MQ=255)
                                                                                     gggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAt          >  1:489492/1‑100 (MQ=255)
                                                                                       gTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTc        >  2:409872/1‑100 (MQ=255)
                                                                                        tCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCa       <  1:2434280/100‑1 (MQ=255)
                                                                                        tCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCa       <  1:1776821/100‑1 (MQ=255)
                                                                                          ggCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACGCAGTGATGGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGACTCAcc     >  2:892914/1‑100 (MQ=255)
                                                                                          ggCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCAcc     <  1:848968/100‑1 (MQ=255)
                                                                                            ccAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGc   >  2:375286/1‑100 (MQ=255)
                                                                                            ccAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGc   >  1:133221/1‑100 (MQ=255)
                                                                                             cAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGcc  >  2:2105827/1‑100 (MQ=255)
                                                                                                 |                                                                                               
CTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGCC  >  NC_000913/1826934‑1827126

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