Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,090,860 G→A 55.0% A223A (GCG→GCA hisG → ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,090,8600GA55.0% 31.3 / 167.4 133A223A (GCG→GCAhisGATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base A (42/30);  total (73/60)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 136 total aligned reads displayed.

CGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCT  >  NC_000913/2090763‑2090958
                                                                                                 |                                                                                                  
cGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGgcacg                                                                                                  >  1:1699234/1‑100 (MQ=255)
   aaGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGCATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCCGGGTGTGGTCCAGGCACGCGa                                                                                               >  2:1178239/1‑100 (MQ=255)
     gcctgcctGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAAt                                                                                             <  2:3012026/100‑1 (MQ=255)
      cctgcctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATc                                                                                            >  2:3237442/1‑100 (MQ=255)
      cctgcctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATc                                                                                            >  1:2854743/1‑100 (MQ=255)
         gcctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCaaa                                                                                         >  1:2749360/1‑100 (MQ=255)
         gcctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCaaa                                                                                         <  2:435097/100‑1 (MQ=255)
          cctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGATCTGTATTCAGGGTGTTATCCAGGCACGGGAATCaaaa                                                                                        >  1:608939/1‑100 (MQ=255)
           ctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAAt                                                                                       <  1:885915/100‑1 (MQ=255)
           ctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAAt                                                                                       <  2:813765/100‑1 (MQ=255)
            tgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATa                                                                                      <  1:1468814/100‑1 (MQ=255)
            tgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGCATCAAAAt                                                                                       >  2:5172/1‑99 (MQ=255)
            tgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAAt                                                                                       >  1:310451/1‑99 (MQ=255)
             gATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATAc                                                                                     >  1:147959/1‑100 (MQ=255)
             gATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATAc                                                                                     >  1:3128350/1‑100 (MQ=255)
             gATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATAc                                                                                     >  1:2450626/1‑100 (MQ=255)
               ttCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATAcat                                                                                   >  1:1588957/1‑100 (MQ=255)
                tCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATAcatc                                                                                  <  2:2277462/100‑1 (MQ=255)
                 cAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATAcatca                                                                                 <  2:1129769/100‑1 (MQ=255)
                 cAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATAcatca                                                                                 >  2:3994349/1‑100 (MQ=255)
                  aaCGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATAcatcat                                                                                >  1:112778/1‑100 (MQ=255)
                  aaCGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATAcatcat                                                                                <  1:3852704/100‑1 (MQ=255)
                  aaCGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATAcatcat                                                                                >  1:3792736/1‑100 (MQ=255)
                  aaCGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATAcatcat                                                                                >  1:3791906/1‑100 (MQ=255)
                    cgcgATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCCTATTCAGGGTGTGATCCAGGCACGCGCCTCACAATACATCatga                                                                              >  2:2351876/1‑100 (MQ=255)
                      cgATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCatgatg                                                                            >  2:2636636/1‑100 (MQ=255)
                      cgATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCatgatg                                                                            <  1:1329994/100‑1 (MQ=255)
                      cgATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCatgatg                                                                            <  1:748820/100‑1 (MQ=255)
                      cgATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCatgatg                                                                            >  1:1416573/1‑100 (MQ=255)
                       gATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATgc                                                                           <  2:1168945/100‑1 (MQ=255)
                        aTGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAGCTGCTGACCCGTATTCAGGGTGTGATCTAGGCGCGCGAATCAAAATACATCATGATgc                                                                           <  2:253540/99‑1 (MQ=255)
                        aTGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATgca                                                                          >  2:2221541/1‑100 (MQ=255)
                          ggCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATgcacg                                                                        <  2:3953071/100‑1 (MQ=255)
                               aaTGGAAGAATCCAAACAGCAACTTATCGACAAACTGCTGACCCTTATTCAGGTTTGGATCCAGGCGCTCGAATCAAAATACATCATGATGCACGCacc                                                                    <  2:2036093/99‑1 (MQ=255)
                                aTGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCaccga                                                                  <  1:2636652/100‑1 (MQ=255)
                                  ggAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCaccgacc                                                                <  2:2427852/100‑1 (MQ=255)
                                  ggAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCaccgacc                                                                >  2:1499810/1‑100 (MQ=255)
                                  ggAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCaccgacc                                                                <  1:2431553/100‑1 (MQ=255)
                                      gaaTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAAc                                                            <  2:3051395/100‑1 (MQ=255)
                                       aaTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACg                                                           <  1:3840371/100‑1 (MQ=255)
                                        aTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGt                                                          >  1:2189136/1‑100 (MQ=255)
                                         tacAAACAGCAACAGATCGACAAACTGCTGAACCGTATTCAGTGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTc                                                         >  1:2788103/3‑100 (MQ=255)
                                           cAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTg                                                       <  2:2463567/100‑1 (MQ=255)
                                           cAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTg                                                       >  2:3570867/1‑100 (MQ=255)
                                           cAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTg                                                       <  1:3484973/100‑1 (MQ=255)
                                           cAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTg                                                       >  1:829860/1‑100 (MQ=255)
                                            aaaCAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTgg                                                      >  2:1720355/1‑100 (MQ=255)
                                              aCAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGAt                                                    <  1:1313938/100‑1 (MQ=255)
                                                aGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGa                                                  <  1:1061667/100‑1 (MQ=255)
                                                  cAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAg                                                <  1:813978/100‑1 (MQ=255)
                                                  cAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAg                                                <  1:3997070/100‑1 (MQ=255)
                                                   aaCTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGt                                               <  1:2981700/100‑1 (MQ=255)
                                                       gATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATc                                           >  1:435455/1‑100 (MQ=255)
                                                         tCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGc                                         <  2:1327690/100‑1 (MQ=255)
                                                          cGACAAACTGCTGACCCGTATTCAGGGTGTGTTCCAGGCTCGCGGATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGCTGAAGTCATCg                                          >  2:3965457/1‑98 (MQ=255)
                                                          cGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGcc                                        >  2:170184/1‑100 (MQ=255)
                                                           gACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGccc                                       >  2:156066/1‑100 (MQ=255)
                                                           gACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGccc                                       >  1:716365/1‑100 (MQ=255)
                                                           gACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGccc                                       >  1:2297985/1‑100 (MQ=255)
                                                           gACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGccc                                       <  1:3677233/100‑1 (MQ=255)
                                                              aaaCTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctgc                                    >  1:2177822/1‑100 (MQ=255)
                                                               aaCTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctgct                                   <  1:896943/100‑1 (MQ=255)
                                                               aaCTGCTGACCCGTATTCAGGGTGTGAGCCAGGCACGCGAATCAAAATACGTCATGATGCACGGACCGACCGAACGTCTGGATGAAGTCATCGCCctg                                     <  2:2571267/98‑1 (MQ=255)
                                                                   gctgACCCGTATTCAGGGTGTGATCCAGGCACGTGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCa                               <  1:980024/100‑1 (MQ=255)
                                                                   gctgACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCa                               >  1:211620/1‑100 (MQ=255)
                                                                   gctgACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCa                               <  2:3128331/100‑1 (MQ=255)
                                                                      gACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGt                            <  2:1239623/100‑1 (MQ=255)
                                                                      gACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGt                            >  1:1667485/1‑100 (MQ=255)
                                                                      gACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGt                            >  1:252715/1‑100 (MQ=255)
                                                                       aCCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTg                           <  2:2665838/100‑1 (MQ=255)
                                                                         ccGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGcc                         >  1:818708/1‑100 (MQ=255)
                                                                         ccGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGcc                         <  1:2031972/100‑1 (MQ=255)
                                                                          cGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCg                        >  1:3209965/1‑100 (MQ=255)
                                                                             aTTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAAc                     <  1:477994/100‑1 (MQ=255)
                                                                             aTTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAAc                     <  1:1955052/100‑1 (MQ=255)
                                                                             aTTCAGAGTGAGATCCAGGCGCGCGAATCAAAATACATAATGATGCACGCACCGACCGAACGTCCAGATGAAGCCAACACCCTGCTCCCa                               >  2:1347184/1‑90 (MQ=255)
                                                                              ttCAGGGTGTGATCCAGGCGCGCGAATCAAAATTCATCATGATGCACGCACCGACCGCACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGacag                    >  2:3937264/1‑97 (MQ=255)
                                                                              ttCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACg                    >  1:1771438/1‑100 (MQ=255)
                                                                               tCAGGGTGTGATCCATCCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATCAAGTCATCGCCCTGCTGCCAGGGGCCGAACGc                   >  2:400912/1‑100 (MQ=255)
                                                                               tCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGc                   >  1:556004/1‑100 (MQ=255)
                                                                               tCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGc                   <  1:1383422/100‑1 (MQ=255)
                                                                               tCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGc                   <  2:1746013/100‑1 (MQ=255)
                                                                                cAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGcc                  >  1:3949639/1‑100 (MQ=255)
                                                                                cAGGGTGTGATCCAGGCACGCGAATCAAATTACATCATGTTGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGcc                  <  2:1292879/100‑1 (MQ=255)
                                                                                cAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGcc                  <  2:919529/100‑1 (MQ=255)
                                                                                 aGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGccc                 <  2:692445/100‑1 (MQ=255)
                                                                                  gggTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGccc                 >  2:3874939/1‑99 (MQ=255)
                                                                                  gggTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCa                <  1:1720367/100‑1 (MQ=255)
                                                                                     tgtgATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACt             >  1:3661547/1‑100 (MQ=255)
                                                                                     tgtgATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACt             <  1:1658913/100‑1 (MQ=255)
                                                                                       tgATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTAt           >  1:749598/1‑100 (MQ=255)
                                                                                        gATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTAtt          <  1:2624332/100‑1 (MQ=255)
                                                                                         aTCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTc         >  1:3985743/1‑100 (MQ=255)
                                                                                          tCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCt        <  1:663837/100‑1 (MQ=255)
                                                                                           ccAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTg       <  1:1922957/100‑1 (MQ=255)
                                                                                           ccAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTg       >  2:2775375/1‑100 (MQ=255)
                                                                                           ccAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTg       >  1:3998237/1‑100 (MQ=255)
                                                                                            cAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGc      <  2:1919653/100‑1 (MQ=255)
                                                                                            cAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGc      >  1:2301413/1‑100 (MQ=255)
                                                                                               gcgcgcGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGc   <  2:2313978/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
CGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCT  >  NC_000913/2090763‑2090958

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: