Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,096,628 A→G 5.9% Q5Q (CAA→CAG hisI → phosphoribosyl‑AMP cyclohydrolase and phosphoribosyl‑ATP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,096,6280AG 5.9% 208.2 / 3.5 68Q5Q (CAA→CAGhisIphosphoribosyl‑AMP cyclohydrolase and phosphoribosyl‑ATP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (26/38);  new base G (2/2);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

TGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTC  >  NC_000913/2096481‑2096775
                                                                                                                                                   |                                                                                                                                                   
tGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGaacaa                                                                                                                                                     >  5:166178/1‑148 (MQ=255)
 ggAACACTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGaacaac                                                                                                                                                    >  3:237902/1‑148 (MQ=255)
  gAACACTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGaacaacaa                                                                                                                                                  <  1:140954/149‑1 (MQ=255)
     cacTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGt                                                                                                                                               <  3:16388/149‑1 (MQ=255)
        ttCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTcgc                                                                                                                                            >  5:169585/1‑149 (MQ=255)
        ttCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTcgc                                                                                                                                            >  4:146039/1‑149 (MQ=255)
          ccTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTcgc                                                                                                                                            <  6:169585/147‑1 (MQ=255)
           cTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGaa                                                                                                                                         >  8:297660/1‑149 (MQ=255)
              gAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAactg                                                                                                                                      >  7:20464/1‑149 (MQ=255)
                aGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAactgga                                                                                                                                    <  4:300974/149‑1 (MQ=255)
                 gCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAAAATTGGTGAATTAAAAGCGTACCTCGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAactggcg                                                                                                                                   >  6:98249/1‑147 (MQ=255)
                   cTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAactggactg                                                                                                                                 <  8:30414/149‑1 (MQ=255)
                   cTTCCGCGATGCCGACGGTGCCGGCGCGCGGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTAACTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAactggactg                                                                                                                                 <  1:332703/149‑1 (MQ=255)
                     tCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAactggactgg                                                                                                                                >  8:226817/1‑148 (MQ=255)
                       cgcgATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGa                                                                                                                              >  7:316207/1‑148 (MQ=255)
                          gATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGaaaaa                                                                                                                          <  3:168437/149‑1 (MQ=255)
                               cGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACTGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGa                                                                                                                     <  1:138197/149‑1 (MQ=255)
                               cGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGACCAACAACGTCGCGAACTGGACTGGGAAAAAACCGa                                                                                                                     <  3:183201/149‑1 (MQ=255)
                                    ttGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGAc                                                                                                                >  5:155104/1‑149 (MQ=255)
                                    ttGACGGCGCGCTGGCAGATTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGAc                                                                                                                <  2:217377/149‑1 (MQ=255)
                                       aCGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTa                                                                                                             <  4:257175/149‑1 (MQ=255)
                                        cGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTAt                                                                                                            <  7:226817/149‑1 (MQ=255)
                                             cgcTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCTAACTGGACTGGGAAAAAACCGACGGACTTATGCCgg                                                                                                       <  4:237902/149‑1 (MQ=255)
                                                  gCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAGCAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGAtt                                                                                                  >  4:269121/1‑149 (MQ=255)
                                                  gCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAGCAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGAtt                                                                                                  >  4:264704/1‑149 (MQ=255)
                                                    aGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATtgt                                                                                                <  4:141766/149‑1 (MQ=255)
                                                              ttCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATTTGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACAcgc                                                                                       <  2:40217/148‑1 (MQ=255)
                                                              ttCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACAcgcg                                                                                      <  3:4070/149‑1 (MQ=255)
                                                              ttCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACAcgcg                                                                                      <  1:177609/149‑1 (MQ=255)
                                                                 cacaAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACCCAGGGCGTGGAGATCAGGAGATGTTAACAGAACAACAACGTCGCGAACTGGACTCGGAAAAAACCGCCGGACTTATGCCGGTGATTGTGCAACACGCGGTa                                                                                   >  4:181763/1‑149 (MQ=255)
                                                                 cacaAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTa                                                                                   >  3:3593/1‑149 (MQ=255)
                                                                   caaacaaaTAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGCACAACAACGTCGCGAACTGGACTGGGAAAAATCCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTAtg                                                                                 >  2:139600/1‑148 (MQ=255)
                                                                    aaacaaaTAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATcc                                                                                >  2:135474/1‑149 (MQ=255)
                                                                      acaaaTAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCgg                                                                              >  5:19055/1‑149 (MQ=255)
                                                                          ataatCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGaa                                                                          <  1:8941/149‑1 (MQ=255)
                                                                          ataatCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGaa                                                                          <  3:146039/149‑1 (MQ=255)
                                                                           taatCAATATTGGTGAATTAAAAGCGTACCTCGCAACACAGGTCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGCGAAAAAACCGACGGCCTTATGCCGGTGAATGTGCAACAAGCGGAAACCGGCGAAg                                                                         >  1:171735/1‑149 (MQ=255)
                                                                              tCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCAGCGAAGTGc                                                                      <  6:108294/149‑1 (MQ=255)
                                                                                 atatTGGTGAATTAACAGCGTACCTGGCAACACAGGGCGTGGAGATCCGGATATGTTAACAGAACAACAACGTCGCGAACTGCACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTaa                                                                   <  6:272393/149‑1 (MQ=255)
                                                                                 atatTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTaa                                                                   <  6:166178/149‑1 (MQ=255)
                                                                                    ttGGTGAATTAAAAGCGTACCTGTCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGc                                                                <  3:193834/149‑1 (MQ=255)
                                                                                    ttGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAGCAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATg                                                                 <  3:269121/148‑1 (MQ=255)
                                                                                    ttGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAGCAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATg                                                                 <  3:264704/148‑1 (MQ=255)
                                                                                       gTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTgg                                                             <  7:325152/149‑1 (MQ=255)
                                                                                       gTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTgg                                                             >  2:297713/1‑149 (MQ=255)
                                                                                        tGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGCCGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTggg                                                            >  2:166466/1‑149 (MQ=255)
                                                                                            ttAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCta                                                         >  4:282062/1‑148 (MQ=255)
                                                                                              aaaaGCGTACCTGGCAACACAGGGCGTGGAGACCAGGATATGTTAACAGAACAACAACGTCGCGCACTGGACTGGGAAAAAACCGACGGACGTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCtata                                                       <  5:243048/148‑1 (MQ=255)
                                                                                               aaaGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTatat                                                      <  2:96845/148‑1 (MQ=255)
                                                                                                  gCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAAc                                                  <  4:168233/149‑1 (MQ=255)
                                                                                                       ccTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGa                                             <  2:106151/149‑1 (MQ=255)
                                                                                                        cTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCtata                                                       >  4:63988/1‑138 (MQ=255)
                                                                                                        cTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCtata                                                       <  3:63988/138‑1 (MQ=255)
                                                                                                        cTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCgg                                              >  2:114093/1‑147 (MQ=255)
                                                                                                        cTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGaa                                            <  6:155104/149‑1 (MQ=255)
                                                                                                         tGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTg                                                                       >  3:260946/1‑121 (MQ=255)
                                                                                                         tGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTg                                                                       <  4:260946/121‑1 (MQ=255)
                                                                                                                acaGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGa                                    >  4:242404/1‑149 (MQ=255)
                                                                                                                 cgatctgTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTa                                      <  3:52062/140‑1 (MQ=255)
                                                                                                                     gCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACaaaa                               <  6:211320/149‑1 (MQ=255)
                                                                                                                        tGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTa                                      >  4:52062/1‑139 (MQ=255)
                                                                                                                         ggAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCt                           <  8:316207/149‑1 (MQ=255)
                                                                                                                          gagaTCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCt                           <  8:20464/148‑1 (MQ=255)
                                                                                                                          gagaTCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTc                          <  6:19055/149‑1 (MQ=255)
                                                                                                                            gaTCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGt                                                                                    >  7:110561/1‑89 (MQ=255)
                                                                                                                            gaTCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGt                                                                                    <  8:110561/89‑1 (MQ=255)
                                                                                                                            gaTCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCg                         >  5:99573/1‑148 (MQ=255)
                                                                                                                                 ggATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCg                   <  7:297660/149‑1 (MQ=255)
                                                                                                                                 ggATATGTTAACACAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCTGTATCCGGCGAAGTGCTAATGCTGTGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCg                   <  7:256428/149‑1 (MQ=255)
                                                                                                                                       gTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACACAACCCTCGAAAGCGGCAAAg             >  4:275202/1‑149 (MQ=255)
                                                                                                                                             aGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACct       >  1:222156/1‑149 (MQ=255)
                                                                                                                                               aacaacaacGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACcttc     >  2:319438/1‑149 (MQ=255)
                                                                                                                                                 caacaacGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACcttctt   >  3:139390/1‑149 (MQ=255)
                                                                                                                                                  aacaacGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCttcttc  >  8:141721/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCAGGATATGTTAACAGAACAACAACGTCGCGAACTGGACTGGGAAAAAACCGACGGACTTATGCCGGTGATTGTGCAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTC  >  NC_000913/2096481‑2096775

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