Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,232,174 G→C 5.1% A111P (GCC→CCC)  cdd → cytidine/deoxycytidine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,232,1740GC 5.1% 265.8 / 5.5 79A111P (GCC→CCC) cddcytidine/deoxycytidine deaminase
Reads supporting (aligned to +/- strand):  ref base G (43/30);  new base C (2/2);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

TGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGCCGCGAGGCACACGCGCTGCGTGAC  >  NC_000913/2232030‑2232320
                                                                                                                                                |                                                                                                                                                  
tgtgCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTg                                                                                                                                                >  5:293304/1‑149 (MQ=255)
   gcgcGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGtt                                                                                                                                             >  7:325868/1‑149 (MQ=255)
    cgcgTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGtt                                                                                                                                             <  5:324706/148‑1 (MQ=255)
     gcgTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGa                                                                                                                                           <  6:293304/149‑1 (MQ=255)
       gTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgt                                                                                                                                         <  5:367416/149‑1 (MQ=255)
       gTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATAGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCCCGCCTGGTTGAgt                                                                                                                                         >  6:298034/1‑149 (MQ=255)
          cGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgtggt                                                                                                                                      <  5:184951/149‑1 (MQ=255)
          cGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgtggt                                                                                                                                      <  1:75730/149‑1 (MQ=255)
             cATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCACACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGa                                                                                                                                    >  4:109925/1‑148 (MQ=255)
             cATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgtggt                                                                                                                                      <  6:60675/146‑1 (MQ=255)
             cATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgtggt                                                                                                                                      >  5:60675/1‑146 (MQ=255)
              aTTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGAGCGATCAGCAACGCCTGGTTGAGTGGTGaa                                                                                                                                   >  5:13655/1‑148 (MQ=255)
                 gTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGaaaaa                                                                                                                                <  3:231736/148‑1 (MQ=255)
                  tCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAgc                                                                                                                              >  1:341190/1‑149 (MQ=255)
                    gAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgtgtt                                                                                                                                      >  1:298702/1‑137 (MQ=255)
                    gAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAgtgtt                                                                                                                                      <  2:298702/139‑3 (MQ=255)
                    gAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAgcgc                                                                                                                            <  7:331226/149‑1 (MQ=255)
                    gAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAgcgc                                                                                                                            <  5:33837/149‑1 (MQ=255)
                     aaTTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCt                                                                                                                           >  5:265519/1‑149 (MQ=255)
                     aaTTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCt                                                                                                                           >  8:53671/1‑149 (MQ=255)
                      aTTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCtt                                                                                                                          <  3:109925/149‑1 (MQ=255)
                       ttttAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTgg                                                                                                                                               <  4:152460/127‑1 (MQ=255)
                       ttttAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTgg                                                                                                                                               >  3:152460/1‑127 (MQ=255)
                         ttAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTg                                                                                                                         >  8:6750/1‑147 (MQ=255)
                          tAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAg                                                                                                                      >  1:6071/1‑149 (MQ=255)
                            aTGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGcc                                                                                                                    <  1:391462/149‑1 (MQ=255)
                                 ggCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCa                                                                                                                <  1:318933/148‑1 (MQ=255)
                                  gcgcAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCAcc                                                                                                              >  2:276328/1‑149 (MQ=255)
                                        ttGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAAc                                                                                                        <  3:29874/149‑1 (MQ=255)
                                           cgcgcgGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTac                                                                                                     <  8:325868/149‑1 (MQ=255)
                                            gcgcgGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTaca                                                                                                    >  8:183697/1‑149 (MQ=255)
                                               cgGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGc                                                                                                 <  7:6750/149‑1 (MQ=255)
                                                      aGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTtgtg                                                                                           <  5:364272/148‑1 (MQ=255)
                                                      aGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGACTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACCCGCATTGTgg                                                                                          >  1:126098/1‑149 (MQ=255)
                                                       gCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGt                                                                                         <  7:53671/149‑1 (MQ=255)
                                                          gAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCAc                                                                                      >  6:148310/1‑149 (MQ=255)
                                                              cTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGcc                                                                                  <  7:220740/149‑1 (MQ=255)
                                                              cTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGcc                                                                                  <  1:276328/149‑1 (MQ=255)
                                                               tGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCg                                                                                 >  2:358538/1‑149 (MQ=255)
                                                               tGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCg                                                                                 >  7:151830/1‑149 (MQ=255)
                                                                ggTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGt                                                                                >  6:271135/1‑149 (MQ=255)
                                                                ggTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGt                                                                                >  1:232717/1‑149 (MQ=255)
                                                                ggTATTTCGGTGCAAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCAACGCCTGGTGGCGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGt                                                                                >  6:363029/1‑149 (MQ=255)
                                                                         gTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGtt                                                                           <  2:248387/145‑1 (MQ=255)
                                                                         gTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGtt                                                                           >  1:248387/1‑145 (MQ=255)
                                                                             cAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGTTCAGCCACCCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCa                                                                              >  2:367952/1‑138 (MQ=255)
                                                                             cAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGTTCAGCCACCCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCa                                                                              >  2:362559/1‑138 (MQ=255)
                                                                             cAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGTTCAGCCACCCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCa                                                                              <  1:367952/138‑1 (MQ=255)
                                                                             cAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGTTCAGCCACCCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCa                                                                              <  1:362559/138‑1 (MQ=255)
                                                                              aaTATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTatg                                                                       >  2:274919/1‑144 (MQ=255)
                                                                              aaTATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTatg                                                                       <  1:274919/144‑1 (MQ=255)
                                                                                    gAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAAtgaactgaa                                                            >  1:157120/1‑149 (MQ=255)
                                                                                               tGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTgg                                                 <  2:232717/149‑1 (MQ=255)
                                                                                               tGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTg                                                  >  7:202691/1‑148 (MQ=255)
                                                                                               tGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACAGAACAGCGGTCTg                                                  >  5:145713/1‑148 (MQ=255)
                                                                                                gCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGa                                                >  5:114205/1‑149 (MQ=255)
                                                                                                  gACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACCCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATc                                              >  8:389282/1‑149 (MQ=255)
                                                                                                   tcAATGCAGCAAACCGTTCATGCCGAACAAAGCGCTATCAGCCACTCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCt                                             <  5:298034/148‑1 (MQ=255)
                                                                                                       tGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACTCCTGGTTGAGTGGTTAAAAAGCTCTTGCAGCCAAAACCGTTAACTACAAGCCTTCTGGTCACATCCGTCACTTTAAGAATGAACTGAACAGCCGTCTCGATCTGCg                                          >  8:372204/1‑148 (MQ=255)
                                                                                                       tGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGcc                                                                                                                    <  8:343763/74‑1 (MQ=255)
                                                                                                       tGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGcc                                                                                                                    >  7:343763/1‑74 (MQ=255)
                                                                                                         cagcaAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTAt                                       <  7:372204/149‑1 (MQ=255)
                                                                                                          agcaAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTAtt                                      >  5:175569/1‑149 (MQ=255)
                                                                                                           gcaAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTc                                     <  8:151830/149‑1 (MQ=255)
                                                                                                            caAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCa                                    >  4:118069/1‑149 (MQ=255)
                                                                                                               aCCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCt                                 <  5:148310/149‑1 (MQ=255)
                                                                                                                 cGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGc                               >  6:234381/1‑149 (MQ=255)
                                                                                                                 cGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGc                               >  8:71549/1‑149 (MQ=255)
                                                                                                                  gTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGcc                              >  1:340313/1‑149 (MQ=255)
                                                                                                                      aTGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGc                          >  6:223724/1‑149 (MQ=255)
                                                                                                                        gcagAACAAAGAGCGATCCGCCACGCCTGGTTGAGTGGTGAAAAAGCGCATGCAGCAATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTGTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATATCCCGGCCcg                        >  2:131457/4‑149 (MQ=255)
                                                                                                                           gAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGACAGGAACTGCGTATTCATCTGCCGTGCcgc                       >  8:310937/1‑147 (MQ=255)
                                                                                                                              cAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGCCGCGAgg                   >  3:6771/1‑148 (MQ=255)
                                                                                                                                 aGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGCCGCGAGGcaca               >  1:316423/1‑149 (MQ=255)
                                                                                                                                    gcgATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTa                                        <  1:346523/121‑1 (MQ=255)
                                                                                                                                            ccACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGCCGCGAGGCACACGCGCTGCGt     >  3:217438/1‑148 (MQ=255)
                                                                                                                                            ccACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGCCGCGAGGCACACGCGCTGCGt     >  2:350266/1‑148 (MQ=255)
                                                                                                                                              aCGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGCCGCGAGGCACACGCGCTGCGTGAc  >  7:161604/1‑149 (MQ=255)
                                                                                                                                              aCGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGCCGCGAGGCACACGCGCTGCGTGAc  <  6:114205/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCTGGTTGAGTGGTGAAAAAGCGCTTGCAGCCATCACCGTTAACTACACGCCTTGTGGTCACTGCCGTCAGTTTATGAATGAACTGAACAGCGGTCTGGATCTGCGTATTCATCTGCCGGGCCGCGAGGCACACGCGCTGCGTGAC  >  NC_000913/2232030‑2232320

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