Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,074,948 C→A 7.5% G1312V (GGT→GTT)  putA ← fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,074,9480CA 7.5% 310.2 / 9.6 93G1312V (GGT→GTT) putAfused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base C (42/44);  new base A (2/5);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATC  >  NC_000913/1074800‑1075096
                                                                                                                                                    |                                                                                                                                                    
aCGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAcc                                                                                                                                                     <  2:799843/150‑1 (MQ=255)
  gTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGc                                                                                                                                                   <  2:1451249/150‑1 (MQ=255)
    tGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCca                                                                                                                                                 >  1:1754471/1‑150 (MQ=255)
      cgcgATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagc                                                                                                                                               <  1:333537/150‑1 (MQ=255)
      cgcgATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAAGGCcagc                                                                                                                                               >  2:196020/1‑150 (MQ=255)
        cgATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagctt                                                                                                                                             <  2:1109342/150‑3 (MQ=255)
        cgATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagctt                                                                                                                                             <  1:9378/150‑3 (MQ=255)
         gATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGAGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggc                                                                                                                                            <  2:1374963/150‑1 (MQ=255)
          aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAAGGCcagcttaa                                                                                                                                           <  1:196020/150‑5 (MQ=255)
          aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagcttaa                                                                                                                                           <  1:551319/150‑5 (MQ=255)
          aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagcttaa                                                                                                                                           <  1:1493127/150‑5 (MQ=255)
             aGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggcagcg                                                                                                                                        <  1:1159657/150‑1 (MQ=255)
              ggTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggcagcgg                                                                                                                                       <  1:1852503/150‑1 (MQ=255)
                   ttttCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTc                                                                                                                                  <  1:210679/150‑1 (MQ=255)
                     ttCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTcac                                                                                                                                <  1:25357/150‑1 (MQ=255)
                       caatcaaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCacac                                                                                                                              >  1:1028456/1‑150 (MQ=255)
                       caatcaaACGGTTGAAAACGTGCAGCGGAGTATTATTGAGGGGATTCGCGCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACGTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCacac                                                                                                                              >  2:1642567/1‑150 (MQ=255)
                        aatcaaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACt                                                                                                                             <  2:906320/150‑1 (MQ=255)
                        aatcaaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACt                                                                                                                             <  1:212826/150‑1 (MQ=255)
                         atcaaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTc                                                                                                                            <  2:1378167/150‑1 (MQ=255)
                         atcaaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTc                                                                                                                            >  1:311964/1‑150 (MQ=255)
                            aaaCGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGc                                                                                                                         >  1:1206563/1‑150 (MQ=255)
                            aaaCGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGc                                                                                                                         >  1:353193/1‑150 (MQ=255)
                            aaaCGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGc                                                                                                                         >  2:928872/1‑150 (MQ=255)
                                ggTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAAc                                                                                                                     >  1:1209851/1‑150 (MQ=255)
                                 gTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACg                                                                                                                    >  2:1030937/1‑150 (MQ=255)
                                 gTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACg                                                                                                                    >  2:888588/1‑150 (MQ=255)
                                     aaaaCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCg                                                                                                                <  2:933554/150‑1 (MQ=255)
                                            gCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACa                                                                                                         >  1:830797/1‑150 (MQ=255)
                                               gCGGATTATTATTGAGGGGATTTGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCAGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGcc                                                                                                      >  2:1455137/1‑150 (MQ=255)
                                               gCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGat                                                                                                              >  2:588900/1‑142 (MQ=255)
                                               gCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGat                                                                                                              <  1:588900/142‑1 (MQ=255)
                                                 ggATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGt                                                                                                    <  2:587070/150‑1 (MQ=255)
                                                    ttattattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                                                                                 <  2:1961645/150‑1 (MQ=255)
                                                    ttattattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                                                                                 <  2:1206563/150‑1 (MQ=255)
                                                    ttattattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                                                                                 <  1:363334/150‑1 (MQ=255)
                                                     tattattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCa                                                                                                <  1:1375699/150‑1 (MQ=255)
                                                      attattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAg                                                                                               <  1:1171035/150‑1 (MQ=255)
                                                         attGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAg                                                                                            <  2:1033547/150‑1 (MQ=255)
                                                          ttGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAgg                                                                                           <  2:980270/150‑1 (MQ=255)
                                                            gAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGat                                                                                         <  2:1637308/150‑1 (MQ=255)
                                                             aGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                                                                        >  2:1782578/1‑150 (MQ=255)
                                                               gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATAtt                                                                                      <  2:1515867/150‑1 (MQ=255)
                                                                 gATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGc                                                                                    <  1:194103/150‑1 (MQ=255)
                                                                 gATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGc                                                                                    >  1:609473/1‑150 (MQ=255)
                                                                 gATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGc                                                                                    >  2:213960/1‑150 (MQ=255)
                                                                 gATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGc                                                                                    >  1:847950/1‑150 (MQ=255)
                                                                  aTTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCt                                                                                   >  1:957033/1‑150 (MQ=255)
                                                                      gACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCg                                                                               >  1:402951/1‑150 (MQ=255)
                                                                        ccACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGc                                                                                                       >  2:1138494/1‑124 (MQ=255)
                                                                        ccACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGc                                                                                                       <  1:1138494/124‑1 (MQ=255)
                                                                           cGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACg                                                                          >  1:813687/1‑150 (MQ=255)
                                                                               aaaGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCa                                                                      >  2:973116/1‑150 (MQ=255)
                                                                                aaGATGCCGGAGGAGGTTGTAACATCCTCCGTCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCTGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaa                                                                     <  1:840342/150‑1 (MQ=255)
                                                                                  gtagCCGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaaaa                                                                   >  1:689424/4‑150 (MQ=255)
                                                                                   aTGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAc                                                                  <  1:259231/150‑1 (MQ=255)
                                                                                    tGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGt                                                                                                    >  2:872752/1‑115 (MQ=255)
                                                                                    tGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGt                                                                                                    <  1:872752/115‑1 (MQ=255)
                                                                                    agCCGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAcc                                                                 >  1:1671537/2‑150 (MQ=255)
                                                                                      ccGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCt                                                               >  1:363226/1‑150 (MQ=255)
                                                                                            gaggTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCg                                                         <  2:169982/150‑1 (MQ=255)
                                                                                                  gTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAc                                                                  >  2:6385/1‑135 (MQ=255)
                                                                                                  gTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAc                                                                  <  1:6385/135‑1 (MQ=255)
                                                                                                    aaCATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCGCCGAAACAAtt                                                 >  1:148192/1‑150 (MQ=255)
                                                                                                       aTCCTCCTGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgtg                                              <  2:1754471/150‑1 (MQ=255)
                                                                                                       aTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgtg                                              <  2:813687/150‑1 (MQ=255)
                                                                                                        tcctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGc                                             <  1:1030937/150‑1 (MQ=255)
                                                                                                        tcctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGc                                             <  1:15660/150‑1 (MQ=255)
                                                                                                          ctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCa                                           <  2:311964/150‑1 (MQ=255)
                                                                                                            ccGGCTACCAGTATAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTTCTTTCGCCACGGGCAAACCCCTGCACCGAAACAATTGTGCCaac                                         >  2:1951777/1‑148 (MQ=255)
                                                                                                             cGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATcc                                        >  2:553891/1‑150 (MQ=255)
                                                                                                             cGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATcc                                        <  1:1875074/150‑1 (MQ=255)
                                                                                                              ggCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATccc                                       >  2:604559/1‑150 (MQ=255)
                                                                                                                  aCCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCcgcgc                                   >  1:834309/1‑150 (MQ=255)
                                                                                                                     tGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGc                                >  1:230367/1‑150 (MQ=255)
                                                                                                                           aCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGc                          >  1:1610310/1‑150 (MQ=255)
                                                                                                                            ccTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCt                         <  1:1100596/150‑1 (MQ=255)
                                                                                                                            ccTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCt                         >  1:28368/1‑150 (MQ=255)
                                                                                                                              tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTc                       <  2:1610310/150‑1 (MQ=255)
                                                                                                                              tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTc                       <  2:1400455/150‑1 (MQ=255)
                                                                                                                                taGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTcac                     <  2:609473/150‑1 (MQ=255)
                                                                                                                                  gTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTcacac                   <  2:334207/150‑1 (MQ=255)
                                                                                                                                     aTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAAt                >  1:230278/1‑150 (MQ=255)
                                                                                                                                         cgCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATgcgc            >  2:244877/2‑150 (MQ=255)
                                                                                                                                         aGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATgcgc            >  1:1756252/1‑150 (MQ=255)
                                                                                                                                          gCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGcgcg           >  1:183867/1‑150 (MQ=255)
                                                                                                                                          gCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGcgcg           <  1:82294/150‑1 (MQ=255)
                                                                                                                                           cTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGa          <  2:1626498/150‑1 (MQ=255)
                                                                                                                                           cTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGa          >  2:1579467/1‑150 (MQ=255)
                                                                                                                                              gCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGc       >  1:687400/1‑150 (MQ=255)
                                                                                                                                                gTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTg     <  1:888588/150‑1 (MQ=255)
                                                                                                                                                   aCCGCCAGCGGCAGCGGTATTCACACTCAGCGCACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGaa         >  2:1267441/1‑143 (MQ=255)
                                                                                                                                                   aCCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGaa         <  1:1267441/143‑1 (MQ=255)
                                                                                                                                                   aCCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATc  >  1:883832/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATC  >  NC_000913/1074800‑1075096

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