Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,078,600 T→C 5.4% I95V (ATC→GTC) ‡ putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,078,6000TC 5.4% 311.7 / 5.5 93I95V (ATC→GTC) ‡putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (40/48);  new base C (2/3);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

CCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAGCGCAGGTAGCTCCGGCAGAGTATCGCTGTT  >  NC_000913/1078455‑1078747
                                                                                                                                                 |                                                                                                                                                   
ccACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATcgc                                                                                                                                                  <  4:230562/149‑1 (MQ=255)
 cACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATcgcg                                                                                                                                                 >  3:235253/1‑149 (MQ=255)
    tGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATcgcggcg                                                                                                                                              >  3:123252/1‑149 (MQ=255)
    tGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATcgcggcg                                                                                                                                              >  4:498356/1‑149 (MQ=255)
    tGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGATCACcgcggc                                                                                                                                               <  1:29412/148‑1 (MQ=255)
     ggCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATcgcggcg                                                                                                                                              >  2:827676/1‑148 (MQ=255)
     ggCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATcgcggcg                                                                                                                                              >  2:181484/1‑148 (MQ=255)
     ggCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATcgcggcg                                                                                                                                              >  1:80883/1‑148 (MQ=255)
       cATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATcgcggcgcgg                                                                                                                                           <  4:748574/149‑1 (MQ=255)
            ttttGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGaaa                                                                                                                                       <  2:864130/148‑1 (MQ=255)
            ttttGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAAc                                                                                                                                      <  4:382954/149‑1 (MQ=255)
            ttttGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAAc                                                                                                                                      <  3:38069/149‑1 (MQ=255)
            ttttGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGAGGGaaa                                                                                                                                       <  2:494660/148‑1 (MQ=255)
             tttGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAAc                                                                                                                                      <  2:279984/148‑1 (MQ=255)
                gATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGa                                                                                                                                   <  1:848121/148‑1 (MQ=255)
                gATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGAc                                                                                                                                  <  4:606588/149‑1 (MQ=255)
                   tACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTgg                                                                                                                               <  4:123252/149‑1 (MQ=255)
                   tACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTgg                                                                                                                               <  3:113771/149‑1 (MQ=255)
                    aCGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTggg                                                                                                                              >  3:471320/1‑149 (MQ=255)
                         gTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGATCACcgcggcctat                                                                                                                                           >  1:861160/1‑127 (MQ=255)
                         gTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGATCACcgcggcctat                                                                                                                                           <  2:861160/131‑5 (MQ=255)
                           ttATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGATCACcgcggccta                                                                                                                                            >  2:694464/1‑125 (MQ=255)
                           ttATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGATCACcgcggccta                                                                                                                                            <  1:694464/128‑4 (MQ=255)
                            tATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATAtt                                                                                                                      <  4:560892/149‑1 (MQ=255)
                               cGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGc                                                                                                                   <  4:420109/149‑1 (MQ=255)
                               cGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGc                                                                                                                   <  3:510619/149‑1 (MQ=255)
                                ggCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCt                                                                                                                  >  4:492452/1‑149 (MQ=255)
                                 gCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCTCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAat                                                                                                                         >  3:80217/1‑141 (MQ=255)
                                 gCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAat                                                                                                                         <  4:80217/141‑1 (MQ=255)
                                 gCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCt                                                                                                                  <  1:468980/148‑1 (MQ=255)
                                 gCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCt                                                                                                                  <  2:910753/148‑1 (MQ=255)
                                 gCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCt                                                                                                                  >  2:681324/1‑148 (MQ=255)
                                  ccAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCg                                                                                                                <  3:248269/149‑1 (MQ=255)
                                        gATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCaaa                                                                                                           <  2:732880/148‑1 (MQ=255)
                                        gATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAg                                                                                                          >  3:212303/1‑149 (MQ=255)
                                           aCGCCAGTTTGTGCGCCTGTTCAGCAACTGGCGGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCg                                                                                                       >  4:32499/1‑149 (MQ=255)
                                            cGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCg                                                                                                       <  2:141319/148‑1 (MQ=255)
                                             gCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAg                                                                                                     <  4:235253/149‑1 (MQ=255)
                                               cAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAgg                                                                                                    >  2:826618/1‑148 (MQ=255)
                                                aGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGa                                                                                                   >  2:279966/1‑148 (MQ=255)
                                                  tttGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTgg                                                                                                                               >  4:367182/1‑118 (MQ=255)
                                                  tttGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTgg                                                                                                                               <  3:367182/118‑1 (MQ=255)
                                                  tttGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATg                                                                                                >  4:402008/1‑149 (MQ=255)
                                                     gtgCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGc                                                                                              <  1:470459/148‑1 (MQ=255)
                                                     gtgCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCt                                                                                             <  4:634460/149‑1 (MQ=255)
                                                          ccTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGgtgt                                                                                        <  4:121824/149‑1 (MQ=255)
                                                           cTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGtgtg                                                                                       <  3:771490/149‑1 (MQ=255)
                                                                cAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTcc                                                                                  >  4:698676/1‑149 (MQ=255)
                                                                 aGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTgggg                                                                                                                             >  3:682626/1‑105 (MQ=255)
                                                                 aGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTgggg                                                                                                                             <  4:682626/105‑1 (MQ=255)
                                                                  gCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCct                                                                                 >  1:222267/1‑148 (MQ=255)
                                                                  gCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCct                                                                                 <  2:624574/148‑1 (MQ=255)
                                                                    aaCTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGgtgt                                                                                        <  1:94417/139‑1 (MQ=255)
                                                                    aaCTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGgtgt                                                                                        >  2:94417/1‑139 (MQ=255)
                                                                    aaCTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCctct                                                                               <  1:48877/148‑1 (MQ=255)
                                                                    aaCTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTg                                                                              >  4:714169/1‑149 (MQ=255)
                                                                     aCTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTg                                                                              >  2:612892/1‑148 (MQ=255)
                                                                      cTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGcc                                                                            <  3:402008/149‑1 (MQ=255)
                                                                       tGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCg                                                                           >  3:544685/1‑149 (MQ=255)
                                                                          cTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAg                                                                        >  3:229934/1‑149 (MQ=255)
                                                                           tGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCaa                                                                                                                          <  1:373063/98‑1 (MQ=255)
                                                                           tGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCaa                                                                                                                          >  2:373063/1‑98 (MQ=255)
                                                                             cGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCg                                                                     <  3:498356/149‑1 (MQ=255)
                                                                             cGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCg                                                                     <  4:212303/149‑1 (MQ=255)
                                                                              ggCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCg                                                                     <  1:827676/148‑1 (MQ=255)
                                                                                 aGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGc                                                                 >  3:733467/1‑149 (MQ=255)
                                                                                 aGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGc                                                                 <  4:692351/149‑1 (MQ=255)
                                                                                     gggCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTc                                                              <  1:680035/148‑1 (MQ=255)
                                                                                      ggCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCAt                                                            >  3:315671/1‑149 (MQ=255)
                                                                                          tGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGct                                                        >  4:877562/1‑149 (MQ=255)
                                                                                          tGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGc                                                         >  1:646065/1‑148 (MQ=255)
                                                                                            ttCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGctct                                                      >  3:564429/1‑149 (MQ=255)
                                                                                               cagcagCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCa                                                    <  1:181484/148‑1 (MQ=255)
                                                                                                agcagcATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCAt                                                   >  2:390061/1‑148 (MQ=255)
                                                                                                agcagcATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCAt                                                   >  2:782731/1‑148 (MQ=255)
                                                                                                 gcagcaTAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTg                                                 <  4:4190/149‑1 (MQ=255)
                                                                                                       tAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTgg                                                                                      >  2:622080/1‑106 (MQ=255)
                                                                                                       tAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTgg                                                                                      <  1:622080/106‑1 (MQ=255)
                                                                                                        aGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCgcgc                                          >  4:691748/1‑149 (MQ=255)
                                                                                                         gAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTc                                                                                   >  2:343690/1‑107 (MQ=255)
                                                                                                         gAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTc                                                                                   <  1:343690/107‑1 (MQ=255)
                                                                                                         gAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGcc                                         >  4:265240/1‑149 (MQ=255)
                                                                                                          aaaCCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCa                                        >  4:99157/1‑149 (MQ=255)
                                                                                                              cGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGaa                                     <  1:425160/148‑1 (MQ=255)
                                                                                                                 ttCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGTGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAagca                                 <  3:32499/149‑1 (MQ=255)
                                                                                                                  tCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAagca                                 <  1:279966/148‑1 (MQ=255)
                                                                                                                     gTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAgcgc                             <  3:99157/149‑1 (MQ=255)
                                                                                                                     gTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAgcgc                             >  3:54296/1‑149 (MQ=255)
                                                                                                                                  gATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAGCGCAGGTAGCTCCGGc                <  3:698676/149‑1 (MQ=255)
                                                                                                                                     aGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAGCGCAGGTAGCTCCGGCaga             <  4:733467/149‑1 (MQ=255)
                                                                                                                                                aTCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAGCGCAGGTAGCTCCGGCAGAGTATCGCt     >  4:349374/1‑146 (MQ=255)
                                                                                                                                                aTCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAGCGCAGGTAGCTCCGGCAGAGTATCGCt     <  3:349374/146‑1 (MQ=255)
                                                                                                                                                aTCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAGCGCAGGTAGCTCCGGCAGAGTATCGCTGtt  >  3:137553/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTGCGCCTGTTCAGCAACTGGCTGCGGCAGGCGGGCTTGTTCCAGCAGCATAGAAACCGCTTCGGTTTCCGGGCGGCGATAGGCCGCGGTGATCGCGGCGCGGGAAACCGACTGGGGCAATATTTGCTCGGCAAAGTCGAGGAATGGCTGGTGTGGTTCCTCTGCCGGAGTCGGTGCTTCATCGCTCTCATTGGCCGCGCCAGAAAGCAGCGCAGGTAGCTCCGGCAGAGTATCGCTGTT  >  NC_000913/1078455‑1078747

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: