Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,396 A→T 16.7% I184N (ATT→AAT)  trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,3960AT16.7% 461.7 / 75.6 191I184N (ATT→AAT) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base A (89/70);  new base T (24/8);  total (113/78)
Fisher's exact test for biased strand distribution p-value = 5.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 195 total aligned reads displayed.

TGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAA  >  NC_000913/1322297‑1322494
                                                                                                   |                                                                                                  
tGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCaa                                                                                                    <  2:955097/100‑1 (MQ=255)
tGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCaa                                                                                                    <  2:1796208/100‑1 (MQ=255)
 gACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAt                                                                                                   <  2:4151571/100‑1 (MQ=255)
  acttAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAt                                                                                                   <  1:1949719/96‑1 (MQ=255)
   cGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATc                                                                                                  >  2:2402325/1‑99 (MQ=255)
    gTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCAc                                                                                                >  2:4403358/1‑100 (MQ=255)
      aGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCa                                                                                              >  2:1717517/1‑100 (MQ=255)
       gttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCAt                                                                                             >  2:1885490/1‑100 (MQ=255)
        ttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATc                                                                                            <  2:2941920/100‑1 (MQ=255)
          cgttcAGGCGAGCAGTGAGACGTTTTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCa                                                                                           <  2:1873338/99‑1 (MQ=255)
          cgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAATCAGGCTGGCCTTAATACGGTTGCTTTTTTTCTGATGGTCAAGCACCATCAg                                                                                          <  1:844749/100‑1 (MQ=255)
          cgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAg                                                                                          <  2:2297107/100‑1 (MQ=255)
           gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGc                                                                                         >  1:950481/1‑100 (MQ=255)
           gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGc                                                                                         >  1:792456/1‑100 (MQ=255)
           gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGc                                                                                         <  1:3363235/100‑1 (MQ=255)
            ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaa                                                                                        <  1:1115269/100‑3 (MQ=255)
            ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCg                                                                                        <  2:3579317/100‑1 (MQ=255)
            ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCg                                                                                        >  1:2271951/1‑100 (MQ=255)
             tcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaaa                                                                                       <  1:3650129/100‑4 (MQ=255)
              cAGGCGAGCAGTGAGACGCTGTTTTTCTTCCTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaaaa                                                                                      <  2:4077226/100‑5 (MQ=255)
                ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                    <  1:1364227/100‑1 (MQ=255)
                ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                    >  2:2839665/1‑100 (MQ=255)
                    aGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCg                                                                                <  1:2704891/100‑1 (MQ=255)
                    aGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCg                                                                                <  2:1639069/100‑1 (MQ=255)
                    aGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCg                                                                                <  2:1684034/100‑1 (MQ=255)
                      cAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgag                                                                              <  2:2862382/100‑1 (MQ=255)
                       aGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                             <  1:3177949/100‑1 (MQ=255)
                       aGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                             <  2:1580963/100‑1 (MQ=255)
                        gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTTAATCACCATCAGCGTTTTAGCGAGAt                                                                            >  1:1451642/1‑100 (MQ=255)
                        gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGAt                                                                            >  2:3145422/1‑100 (MQ=255)
                         tGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATa                                                                           <  1:3056118/100‑1 (MQ=255)
                         tGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATa                                                                           <  1:4062719/100‑1 (MQ=255)
                              cGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCGGAATACGGGTGCTTTTTTTCTGATGTTCAATTACCATCAGCGTTTCAGCGAGATaa                                                                          >  1:4153824/1‑96 (MQ=255)
                               gttgttTTTCTTCTTCATTCGGATCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAc                                                                      <  2:2568268/99‑1 (MQ=255)
                               gttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAcn                                                                     <  2:3141383/100‑2 (MQ=255)
                                   tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAACACAGaa                                                                  >  2:325913/1‑99 (MQ=255)
                                   tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGaaa                                                                 <  1:3457735/100‑1 (MQ=255)
                                        cttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                            <  2:701124/100‑1 (MQ=255)
                                        cttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                            >  1:1575078/1‑100 (MQ=255)
                                          tcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGc                                                          >  2:3983476/1‑100 (MQ=255)
                                           cttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCa                                                         >  2:1272588/1‑100 (MQ=255)
                                            ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAg                                                        >  1:1868589/1‑100 (MQ=255)
                                            ttcATTCGGACCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTCTCAGCGAGATAAAAACAGAAATCAGGGCAg                                                        <  1:2888394/100‑1 (MQ=255)
                                             tcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGt                                                       >  1:2988493/1‑100 (MQ=255)
                                               aTTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGtt                                                      >  2:588329/1‑99 (MQ=255)
                                               aTTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTCTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTa                                                     >  1:2353612/1‑100 (MQ=255)
                                                 tCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTTTTTTTTTCAATCACCATAAGCGTTTCAGCGACCTAAAAACAGAAATCAGGGCAGTTAtt                                                   >  1:1683646/1‑100 (MQ=255)
                                                 tCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtt                                                   <  2:2717838/100‑1 (MQ=255)
                                                   ggAGCAAACAGGCTGGCCATAATACGGGTGCTTTTTTTCTGATGGTCCATCGCCAACAGCGTCTCAGCGAGAGATAAACAGAGATCAGGGCAGTTAtttt                                                 >  1:3128214/1‑100 (MQ=255)
                                                    gAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTTTTTTTTTTTATTACCAACAGCGGTTTAACGGGATTACGAAAGAAATGAGGGCAGGTAttt                                                  >  1:2556566/1‑98 (MQ=255)
                                                     aGCAAACAGGCTGGCCTTAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTcc                                               <  1:2009129/100‑1 (MQ=255)
                                                       cAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                             <  2:1227010/100‑1 (MQ=255)
                                                       cAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                             <  1:2999913/100‑1 (MQ=255)
                                                       cAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                             >  2:945125/1‑100 (MQ=255)
                                                        aaaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCt                                            >  2:1981675/1‑100 (MQ=255)
                                                        aaaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCt                                            <  1:63546/100‑1 (MQ=255)
                                                             ggCTGGCCTGAATACGGGTGCTTTTTTTTTTTTTGTTTATCACCATCAGCGTTTCAGCGAGATAAATACAGAAATCAGGGCAGTTATTTTCCGCTGACAg                                       >  1:402401/1‑100 (MQ=255)
                                                              gCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGAc                                         <  1:3566257/97‑1 (MQ=255)
                                                                 ggCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGc                                   >  2:4336214/1‑100 (MQ=255)
                                                                 ggCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGc                                   <  2:359063/100‑1 (MQ=255)
                                                                   ccTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                 >  2:1088969/1‑100 (MQ=255)
                                                                   ccTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                 >  2:716805/1‑100 (MQ=255)
                                                                         tACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAAt                            >  2:2078151/1‑99 (MQ=255)
                                                                          aCGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCt                          <  2:3352400/100‑1 (MQ=255)
                                                                           cGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCtt                         >  1:2875388/1‑100 (MQ=255)
                                                                             ggTGCTTTTTTTTTTTTTTTTAATCACCATCAGCGTGTCAGCGATATAAGAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGGGGTAAATCTTc                        >  1:4298214/1‑99 (MQ=255)
                                                                             ggTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCa                       <  2:900902/100‑1 (MQ=255)
                                                                             ggTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCa                       >  2:459993/1‑100 (MQ=255)
                                                                             ggTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCa                       <  2:3721046/100‑1 (MQ=255)
                                                                               tGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCaan                     <  1:3141555/100‑2 (MQ=255)
                                                                                gctgaaacgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                    >  2:2202813/10‑100 (MQ=37)
                                                                                 cTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCATCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATc                   >  1:2157291/1‑100 (MQ=255)
                                                                                  tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATcc                  >  1:3171442/1‑100 (MQ=255)
                                                                                    tttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCg                <  2:2703655/100‑1 (MQ=255)
                                                                                     ttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc               <  2:2136495/100‑1 (MQ=255)
                                                                                     ttgtCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc               <  2:3935692/97‑1 (MQ=255)
                                                                                     aacgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc               >  1:4095116/5‑100 (MQ=255)
                                                                                      tttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGcc              >  2:1590943/1‑100 (MQ=255)
                                                                                       ttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGCAATCAGGGCCGTTATTATCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca             >  1:1929975/1‑100 (MQ=255)
                                                                                       ttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca             >  2:2616257/1‑100 (MQ=255)
                                                                                       cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGACATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGcc              >  1:2411329/3‑99 (MQ=255)
                                                                                       cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGAGAGTTGCGGTACCCATTCAAATCCCg                >  1:1471534/3‑97 (MQ=255)
                                                                                       cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGcc              >  1:132572/3‑99 (MQ=255)
                                                                                       cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca             >  2:3407708/3‑100 (MQ=255)
                                                                                       cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca             >  2:3523595/3‑100 (MQ=255)
                                                                                       cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca             >  1:666099/3‑100 (MQ=255)
                                                                                        gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAGATCTGCAAATCCCGCcac            >  1:4400468/2‑100 (MQ=255)
                                                                                        gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac            >  1:3411988/2‑100 (MQ=255)
                                                                                        gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac            >  2:3030250/2‑100 (MQ=255)
                                                                                        gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac            >  2:2585797/2‑100 (MQ=255)
                                                                                        gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac            >  1:1112344/2‑100 (MQ=255)
                                                                                         cTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac            >  2:628096/1‑99 (MQ=255)
                                                                                          tGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACaa          >  2:4309218/1‑100 (MQ=255)
                                                                                            aTGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAgg        <  2:3422415/100‑1 (MQ=255)
                                                                                             tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt       <  2:3437848/100‑1 (MQ=255)
                                                                                             tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt       >  1:925401/1‑100 (MQ=255)
                                                                                              ggTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTc      >  2:1161236/1‑100 (MQ=255)
                                                                                               gTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAggaca     >  2:3610021/1‑97 (MQ=255)
                                                                                                tCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt    <  1:716803/100‑1 (MQ=255)
                                                                                                 cAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGACAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATa   >  1:707848/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAA  >  NC_000913/1322297‑1322494

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