Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,331,289 C→T A81V (GCG→GTG)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,331,2890CT98.4% 214.3 / NA 63A81V (GCG→GTG) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/31);  minor base A (1/0);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTG  >  NC_000913/1331143‑1331435
                                                                                                                                                  |                                                                                                                                                     
tCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCAc                                                                                                                                                   <  2:844323/151‑1 (MQ=255)
   caTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTag                                                                                                                                                >  1:552616/1‑150 (MQ=255)
   caTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTAt                                                                                                                                                <  2:904900/151‑1 (MQ=255)
   caTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTAt                                                                                                                                                <  2:789773/151‑1 (MQ=255)
   caTCCGCGATTTGCCGACCAGTGGCTCAGCTGCAAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTTAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTAt                                                                                                                                                <  2:529695/151‑1 (MQ=255)
   caTACGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCCCTAt                                                                                                                                                >  2:144964/1‑151 (MQ=255)
     tCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGa                                                                                                                                              <  1:687448/151‑1 (MQ=255)
         cgATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAgtg                                                                                                                                          <  2:483978/151‑1 (MQ=255)
           aTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAgtgt                                                                                                                                         <  2:601703/150‑1 (MQ=255)
           aTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAgtgt                                                                                                                                         <  2:712584/150‑1 (MQ=255)
                ccGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCt                                                                                                                                    <  2:211743/150‑1 (MQ=255)
                ccGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCt                                                                                                                                    <  2:21003/150‑1 (MQ=255)
                    ccAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTaa                                                                                                                               >  2:773836/1‑151 (MQ=255)
                       gTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAaga                                                                                                                            >  1:579751/1‑151 (MQ=255)
                        tGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAagag                                                                                                                           <  2:254510/151‑1 (MQ=255)
                        tGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAagag                                                                                                                           <  1:152467/151‑1 (MQ=255)
                        tGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAaga                                                                                                                            <  1:816597/150‑1 (MQ=255)
                        tGGCTCAGCTGCCAAAAAGAATGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAagag                                                                                                                           <  1:301356/151‑1 (MQ=255)
                          gCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAgaga                                                                                                                          <  1:421314/150‑1 (MQ=255)
                                  nccAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑gtcgt                                                                                                                >  2:147284/2‑149 (MQ=255)
                                    cAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGtct                                                                                                              >  2:138920/1‑149 (MQ=255)
                                      aaaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTctga                                                                                                          >  1:170214/1‑151 (MQ=255)
                                      aaaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTctga                                                                                                          <  2:444893/151‑1 (MQ=255)
                                                  aCTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAaacaa                                                                                                >  2:915510/1‑149 (MQ=255)
                                                  aCTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACt                                                                                              <  1:889349/151‑1 (MQ=255)
                                                           ccaccaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGaaaa                                                                                     >  1:11642/1‑151 (MQ=255)
                                                            caccaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGaaaaa                                                                                    >  1:532169/1‑151 (MQ=255)
                                                              ccaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGc                                                                                  <  1:868361/151‑1 (MQ=255)
                                                                    cGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACca                                                                            >  2:277553/1‑151 (MQ=255)
                                                                      gCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCTTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACcaca                                                                          >  1:96840/1‑151 (MQ=255)
                                                                      gCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGGCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGCAACAACTGGCTGAAAAAGCCGACcaca                                                                          >  1:48080/1‑151 (MQ=255)
                                                                             aGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctat                                                                     >  2:687313/1‑149 (MQ=255)
                                                                               ccTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctatct                                                                   <  1:151036/149‑1 (MQ=255)
                                                                               ccTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctatct                                                                   <  2:233879/149‑1 (MQ=255)
                                                                                cTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctatct                                                                   <  2:335865/151‑1 (MQ=255)
                                                                                           gAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGGGGCACAATTGGGAGGTGCACTATGGAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTATGAACTGAAACAACTGGCTGAAAAAGCGGACAACATCTATCTCGCAACCGACCTTg                                                     >  2:481904/1‑151 (MQ=255)
                                                                                               gTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACcg                                                 >  2:696495/1‑151 (MQ=255)
                                                                                                    gcTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAg                                            >  1:1038572/1‑151 (MQ=255)
                                                                                                     ctctCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAgg                                           <  1:38628/151‑1 (MQ=255)
                                                                                                       ctcGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAgggg                                         >  2:538543/1‑151 (MQ=255)
                                                                                                       ctcGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAgggg                                         >  2:592290/1‑151 (MQ=255)
                                                                                                       ctcGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGACTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAgggg                                         >  2:150042/1‑151 (MQ=255)
                                                                                                       ctcGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGAAGGAGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAgggg                                         >  2:686879/1‑151 (MQ=255)
                                                                                                            cAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGcc                                    <  2:119926/151‑1 (MQ=255)
                                                                                                            cAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGcc                                    <  1:554474/151‑1 (MQ=255)
                                                                                                                cGTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTg                                >  2:369542/1‑151 (MQ=255)
                                                                                                                 gTATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTg                                >  1:445033/1‑150 (MQ=255)
                                                                                                                 gTATGGGGGTTGACCCGTGGCACAATTGGGAAGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTg                                >  1:332271/1‑150 (MQ=255)
                                                                                                                  tATGGGGGTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCa                              >  2:748918/1‑151 (MQ=255)
                                                                                                                      ggggTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATgg                           >  2:801644/1‑150 (MQ=255)
                                                                                                                         gTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCAc                        >  1:1045939/1‑150 (MQ=255)
                                                                                                                         gTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCAc                        <  2:296721/150‑1 (MQ=255)
                                                                                                                         gTTGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCAc                        <  2:216425/150‑1 (MQ=255)
                                                                                                                         gTTGACCCGTGGCACAATTGGGAGGGGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCAc                        <  2:767663/150‑1 (MQ=255)
                                                                                                                          ttGACCCGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCAcc                       <  2:204945/150‑1 (MQ=255)
                                                                                                                              cccGTGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCgg                  >  2:964649/1‑151 (MQ=255)
                                                                                                                                  tGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGCCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAg              <  2:486879/151‑1 (MQ=255)
                                                                                                                                  tGGCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAg              <  1:1030193/151‑1 (MQ=255)
                                                                                                                                    gCACAATTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTg            >  1:162759/1‑151 (MQ=255)
                                                                                                                                         aTTGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTgg       >  2:69067/1‑151 (MQ=255)
                                                                                                                                          ttGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTggg      <  2:53454/151‑1 (MQ=255)
                                                                                                                                          ttGGGAGGTGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTggg      <  1:79742/151‑1 (MQ=255)
                                                                                                                                              gAGGTGCACTATGAAGTGTCGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCCCATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGtg  >  2:161745/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAG‑‑‑GTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTG  >  NC_000913/1331143‑1331435

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑