Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,152,979 G→A G134S (GGT→AGT)  argC → N‑acetyl‑gamma‑glutamylphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,152,9790GA95.7% 625.4 / 6.6 186G134S (GGT→AGT) argCN‑acetyl‑gamma‑glutamylphosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base A (77/101);  total (80/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 194 total aligned reads displayed.

CGGCGCGTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACGGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTG  >  NC_000913/4152888‑4153068
                                                                                           |                                                                                         
cGGCGCGTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACgg                                                                                          >  1:4529169/1‑93 (MQ=255)
 ggCGCGTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGt                                                                                         >  2:4727075/1‑93 (MQ=255)
  gcgcGTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTc                                                                                        <  1:1100612/93‑1 (MQ=255)
   ggcgTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCATTACCCATAATTACCCGGAACTGTTGGAACAGGCAGCCTACAGTCt                                                                                       <  2:1498442/92‑1 (MQ=255)
   cgcgTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCt                                                                                       >  1:1251429/1‑93 (MQ=255)
   cgcgTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCt                                                                                       <  2:1341755/93‑1 (MQ=255)
    gcgTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTg                                                                                      <  2:1587544/93‑1 (MQ=255)
     cgTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTgg                                                                                     <  1:4188190/93‑1 (MQ=255)
      gTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGc                                                                                    >  1:378798/1‑93 (MQ=255)
      gTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGc                                                                                    <  1:587636/93‑1 (MQ=255)
       tttCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCg                                                                                   >  1:2539465/1‑93 (MQ=255)
         tCGTGTTAGCGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGa                                                                                 <  1:3122690/93‑1 (MQ=255)
           gtgtTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgt                                                                               <  1:4039339/93‑1 (MQ=255)
           gtgtTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgt                                                                               >  1:2005204/1‑93 (MQ=255)
            tgtTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtg                                                                              >  2:1563694/1‑93 (MQ=255)
            tgtTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtg                                                                              >  2:3061650/1‑93 (MQ=255)
            tgtTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtg                                                                              <  2:3833769/93‑1 (MQ=255)
             gtTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtgg                                                                             <  1:1461042/93‑1 (MQ=255)
              ttAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtggt                                                                            >  2:1054078/1‑93 (MQ=255)
              ttAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtggg                                                                            >  2:1964295/1‑92 (MQ=255)
              ttAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtggg                                                                            >  1:3095565/1‑92 (MQ=255)
               tAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtggtg                                                                           >  2:3369052/1‑93 (MQ=255)
               tAACGACGCCACCTTCTATGAAAAATAATACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAgtggtg                                                                           >  1:3280567/1‑93 (MQ=255)
                aaCGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGc                                                                          <  1:861353/93‑1 (MQ=255)
                  cgacgCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCgg                                                                        >  1:4621926/1‑93 (MQ=255)
                   gacgCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGt                                                                       <  1:158896/93‑1 (MQ=255)
                   gacgCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGGGGTGCGGt                                                                       >  1:1326988/1‑93 (MQ=255)
                    acgCCCCCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGta                                                                      <  2:1598551/93‑1 (MQ=255)
                       ccACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGGGAGGTGCGGtaata                                                                   >  1:1866945/1‑93 (MQ=255)
                        cACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGtaataa                                                                  >  2:108255/1‑93 (MQ=255)
                         aCCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGACCAGGCAGCCTACAGTCTGCGGGGGGGGTGAGGAAATaaa                                                                 >  2:3677461/1‑93 (MQ=255)
                          ccTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAAt                                                                >  2:1426875/1‑93 (MQ=255)
                           cTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAAtt                                                               >  2:1785348/1‑93 (MQ=255)
                            ttCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTa                                                              <  2:4414249/93‑1 (MQ=255)
                             tCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTaa                                                             <  1:1976339/93‑1 (MQ=255)
                              cTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTaaa                                                            <  1:2731085/93‑1 (MQ=255)
                               tATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACGGTCTGGCGGAGGGGTGCGGTAATAAATTaaaa                                                           >  1:3581376/1‑93 (MQ=255)
                               tATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTaaaa                                                           >  1:3581375/1‑93 (MQ=255)
                               tATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGGAATAAATTaaaa                                                           >  2:3101788/1‑93 (MQ=255)
                                aTGAAAAATATTACGGCTTTACCCATCAATACCGGGACCTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAaag                                                          <  2:121044/93‑1 (MQ=255)
                                  gAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAaagaa                                                        >  2:2294637/1‑93 (MQ=255)
                                   aaaaaTATTACGGCTTTACCCATCAATGCCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAaagaag                                                       <  2:3900963/93‑1 (MQ=255)
                                   aaaaaTATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAaagaag                                                       <  1:1085600/93‑1 (MQ=255)
                                   aaaaaTATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAaagaag                                                       <  2:3796265/93‑1 (MQ=255)
                                    aaaaTATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGc                                                      <  2:2348514/93‑1 (MQ=255)
                                       atatTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGaa                                                   <  1:454777/93‑1 (MQ=255)
                                        tatTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAAt                                                  <  2:1253514/93‑1 (MQ=255)
                                         atTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACGGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAAtt                                                 <  2:4529169/93‑1 (MQ=255)
                                          ttACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCTAAttt                                                <  2:3180752/93‑1 (MQ=255)
                                          ttACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAAttt                                                <  2:4337751/93‑1 (MQ=255)
                                          ttACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAAttt                                                <  2:3924612/93‑1 (MQ=255)
                                           tACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACGGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTg                                               <  2:4037534/93‑1 (MQ=255)
                                            aCGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGa                                              <  1:2761252/93‑1 (MQ=255)
                                             cGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGAt                                             <  1:1677905/93‑1 (MQ=255)
                                              ggCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGAtt                                            <  1:3048477/93‑1 (MQ=255)
                                              ggCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGAtt                                            >  2:2655116/1‑93 (MQ=255)
                                               gCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTg                                           >  1:786555/1‑93 (MQ=255)
                                                 tttACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACGGTCTGGCGGAGAGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCg                                         >  2:4695672/1‑93 (MQ=255)
                                                 tttACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTTCg                                         >  1:2878415/1‑93 (MQ=255)
                                                 tttACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCg                                         >  1:1789094/1‑93 (MQ=255)
                                                 tttACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCg                                         <  1:2469439/93‑1 (MQ=255)
                                                  ttACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCgg                                        >  2:2826137/1‑93 (MQ=255)
                                                  ttACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCgg                                        <  2:1100779/93‑1 (MQ=255)
                                                   tagccATCAATACGCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGt                                       <  2:4238865/90‑1 (MQ=255)
                                                   tACCCATCAATACCCGGAACTGTTGGAACAGCCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGt                                       <  1:3151589/93‑1 (MQ=255)
                                                    aCCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTg                                      <  2:4156524/93‑1 (MQ=255)
                                                          cAATACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGc                                <  2:3747442/93‑1 (MQ=255)
                                                           aaTACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCt                               >  2:2567188/1‑93 (MQ=255)
                                                             tACCCGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGt                             >  1:3965854/1‑93 (MQ=255)
                                                                ccGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTAt                          >  1:3396412/1‑93 (MQ=255)
                                                                ccGGAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTAt                          <  1:2696660/93‑1 (MQ=255)
                                                                 cGGATCTGTTGGAGCAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATc                         <  1:54242/93‑1 (MQ=255)
                                                                  ggAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATcc                        >  2:1646913/1‑93 (MQ=255)
                                                                  ggAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATcc                        >  1:222219/1‑93 (MQ=255)
                                                                   gAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCg                       >  1:2123512/1‑93 (MQ=255)
                                                                   gAACTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCg                       >  2:823380/1‑93 (MQ=255)
                                                                    aaCTGTTGGAACAGGCAGCCTACAGTCTGGCGGCGGGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGGCGGGCTGTTATccta                      >  2:2051769/1‑91 (MQ=255)
                                                                     aCTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGTTGCCGGGCTGTTATCCGAc                     <  2:53232/93‑1 (MQ=255)
                                                                     aCTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGAc                     <  1:3337063/93‑1 (MQ=255)
                                                                     aCTGTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGCAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGAc                     <  1:1153943/93‑1 (MQ=255)
                                                                      cTGTTGGAACAGGCAGCCTACAGTCTGGCGGGGTGGGGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGAcg                    >  2:2876417/1‑93 (MQ=255)
                                                                        gTTGGAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGAcggc                  >  2:353325/1‑93 (MQ=255)
                                                                           ggAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACggcggc               >  1:4005780/1‑93 (MQ=255)
                                                                            gAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGca              >  2:2743601/1‑93 (MQ=255)
                                                                            gAACAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGca              <  2:1756817/93‑1 (MQ=255)
                                                                             aaCAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGcac             <  2:957248/93‑1 (MQ=255)
                                                                             aaCAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGcac             <  1:1912228/93‑1 (MQ=255)
                                                                              aCAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGcaca            <  2:3599607/93‑1 (MQ=255)
                                                                              aCAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGcaca            <  2:2093940/93‑1 (MQ=255)
                                                                               cAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAg           <  2:4466516/93‑1 (MQ=255)
                                                                               cAGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCCGGCTGTTATCCGACGGCGGcacac           >  1:2222460/1‑92 (MQ=255)
                                                                                aGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGc          <  1:1822248/93‑1 (MQ=255)
                                                                                aGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGc          >  1:238031/1‑93 (MQ=255)
                                                                                aGGCAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGc          <  1:3478455/93‑1 (MQ=255)
                                                                                 ggCAGCCTACAGTCTGGCGGAGAGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGTCGGGCTGTTATCCGACGGCGGCACAGcc         >  1:690740/1‑92 (MQ=255)
                                                                                 ggCAGCCTACAGGCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCt         <  2:4298077/93‑1 (MQ=255)
                                                                                   cAGCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTgg       >  1:4448780/1‑93 (MQ=255)
                                                                                     gCCTACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGgcg     <  2:2253236/93‑1 (MQ=255)
                                                                                      ccTACAGTCTGGCGGAGTTGTGCGGTAAGAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGgcgc    <  1:1495986/93‑1 (MQ=255)
                                                                                        tACAGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTg  <  2:2715883/93‑1 (MQ=255)
                                                                                           |                                                                                         
CGGCGCGTTTCGTGTTAACGACGCCACCTTCTATGAAAAATATTACGGCTTTACCCATCAATACCCGGAACTGTTGGAACAGGCAGCCTACGGTCTGGCGGAGTGGTGCGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTG  >  NC_000913/4152888‑4153068

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