Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,047,222 A→G 5.5% D86G (GAC→GGC)  polA → 5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,047,2220AG 5.5% 239.1 / 3.2 73D86G (GAC→GGC) polA5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease
Reads supporting (aligned to +/- strand):  ref base A (31/38);  new base G (2/2);  total (33/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 = 7.65e-01

GATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCA  >  NC_000913/4047079‑4047368
                                                                                                                                               |                                                                                                                                                  
gATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATc                                                                                                                                               <  1:275534/149‑1 (MQ=255)
   gTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGc                                                                                                                                            <  4:467130/149‑1 (MQ=255)
   gTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGc                                                                                                                                            <  3:400753/149‑1 (MQ=255)
             cTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                                  >  4:25086/1‑149 (MQ=255)
             cTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                                  >  3:67512/1‑149 (MQ=255)
             cTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                                  >  2:14065/1‑149 (MQ=255)
              tCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGTCCGGCATTGGCGGGCGAtg                                                                                                                                               <  2:579629/135‑2 (MQ=255)
              tCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGTCCGGCATTGGCGGGCGAtg                                                                                                                                               >  1:579629/1‑134 (MQ=255)
              tCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATc                                                                                                                                 >  4:491198/1‑149 (MQ=255)
              tCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATc                                                                                                                                 >  2:340091/1‑149 (MQ=255)
                aaCATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGa                                                                                                                               >  2:514665/1‑149 (MQ=255)
                aaCATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGa                                                                                                                               >  2:198992/1‑149 (MQ=255)
                 aCATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGaa                                                                                                                              <  2:433313/149‑1 (MQ=255)
                 aCATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGaa                                                                                                                              <  3:50626/149‑1 (MQ=255)
                   aTGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAAcc                                                                                                                            <  1:514665/149‑1 (MQ=255)
                    tgctgcGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAAccc                                                                                                                           <  1:19347/149‑1 (MQ=255)
                      ctgcGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCtt                                                                                                                         >  2:34/1‑149 (MQ=255)
                                tGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATg                                                                                                               >  3:316162/1‑149 (MQ=255)
                                 gATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATgg                                                                                                              <  1:340091/149‑1 (MQ=255)
                                  aTCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGt                                                                                                             >  3:450692/1‑149 (MQ=255)
                                   tCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGaa                                                                                                                              >  1:311128/1‑131 (MQ=255)
                                   tCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGaa                                                                                                                              <  2:311128/131‑1 (MQ=255)
                                   tCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGtt                                                                                                            >  4:262370/1‑149 (MQ=255)
                                    cATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTa                                                                                                           <  2:422139/149‑1 (MQ=255)
                                     aTGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCAcgc                                                                                                                   >  2:164018/1‑140 (MQ=255)
                                     aTGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCAcgc                                                                                                                   <  1:164018/140‑1 (MQ=255)
                                     aTGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTaa                                                                                                          >  3:238254/1‑149 (MQ=255)
                                       gCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAg                                                                                                        >  3:235546/1‑149 (MQ=255)
                                           tataAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTg                                                                                                                        <  2:534095/129‑1 (MQ=255)
                                           tataAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTg                                                                                                                        >  1:534095/1‑129 (MQ=255)
                                           tataAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGAt                                                                                                    >  1:519919/1‑149 (MQ=255)
                                               aaCCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGa                                                                                                >  3:429944/1‑149 (MQ=255)
                                               aaCCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGa                                                                                                >  2:384967/1‑149 (MQ=255)
                                                          gCAGCGGTGGTCTTTGACGCCAAGTGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAg                                                                                                        <  3:51764/130‑1 (MQ=255)
                                                          gCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAg                                                                                                        >  4:51764/1‑130 (MQ=255)
                                                            aGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACACGATGGTTAAAGCGATGGGACTGCCGCTGCTgg                                                                                   <  4:203620/149‑1 (MQ=255)
                                                                       ttGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGc                                                                        <  2:519919/149‑1 (MQ=255)
                                                                         taCGCCAAGGGAAAAACCTTTCTGGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGt                                                                      <  4:92480/148‑1 (MQ=255)
                                                                          aCGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTa                                                                     <  4:316162/149‑1 (MQ=255)
                                                                             ccAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGaa                                                                  <  4:429944/149‑1 (MQ=255)
                                                                                      aaaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGgacgac                                                         <  3:125167/149‑1 (MQ=255)
                                                                                           tttCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTAt                                                    >  4:187038/1‑149 (MQ=255)
                                                                                           tttCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTAt                                                    <  4:67512/149‑1 (MQ=255)
                                                                                           tttCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTAt                                                    <  3:552174/149‑1 (MQ=255)
                                                                                            ttCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATc                                                   <  1:552737/149‑1 (MQ=255)
                                                                                                     aaCTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTg                                          >  2:475173/1‑149 (MQ=255)
                                                                                                       cTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgc                                        >  4:263059/1‑149 (MQ=255)
                                                                                                          tttGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgcgcg                                     >  2:536715/1‑149 (MQ=255)
                                                                                                              aaCATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGaa                                 <  4:103925/149‑1 (MQ=255)
                                                                                                                cATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGc                               <  3:25086/149‑1 (MQ=255)
                                                                                                                  ttACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCg                             >  1:482506/1‑149 (MQ=255)
                                                                                                                   tACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGa                            >  3:498508/1‑149 (MQ=255)
                                                                                                                      aaaTCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGaaaa                         <  1:410612/149‑1 (MQ=255)
                                                                                                                      aaaTCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGaaaa                         >  4:36987/1‑149 (MQ=255)
                                                                                                                       aaTCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGaaaaa                        >  3:398071/1‑149 (MQ=255)
                                                                                                                        aTCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGga                                                             <  2:289660/111‑1 (MQ=255)
                                                                                                                        aTCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGga                                                             >  1:289660/1‑111 (MQ=255)
                                                                                                                           acaTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCg                    <  4:229778/149‑1 (MQ=255)
                                                                                                                           acaTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCg                    <  4:388725/149‑1 (MQ=255)
                                                                                                                            caTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCgg                   <  4:238254/149‑1 (MQ=255)
                                                                                                                            caTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCgg                   <  4:235546/149‑1 (MQ=255)
                                                                                                                             aTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCggg                  >  2:269211/1‑149 (MQ=255)
                                                                                                                             aTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCg                    >  3:45418/1‑147 (MQ=255)
                                                                                                                             aTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAACGAACCCTTTCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCg                    <  4:45418/147‑1 (MQ=255)
                                                                                                                                 cccgccAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTc              <  4:256150/149‑1 (MQ=255)
                                                                                                                                 cccgccAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTc              <  1:14065/149‑1 (MQ=255)
                                                                                                                                   cgccAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCg            <  3:491198/149‑1 (MQ=255)
                                                                                                                                     ccAATGCCGGACGATCTGCGTACACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTActc                                            <  1:284297/115‑1 (MQ=255)
                                                                                                                                     ccAATGCCGGACGATCTGCGTACACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTActc                                            >  2:284297/1‑115 (MQ=255)
                                                                                                                                      cAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTg         <  4:450692/149‑1 (MQ=255)
                                                                                                                                        aTGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCt       <  1:198992/149‑1 (MQ=255)
                                                                                                                                            cGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATc   <  4:50657/149‑1 (MQ=255)
                                                                                                                                             ggACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCa  <  1:378004/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCA  >  NC_000913/4047079‑4047368

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: ‑