Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,441,367 A→G 7.5% N169D (AAT→GAT)  topA → DNA topoisomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,441,3670AG 7.5% 189.6 / 4.1 66N169D (AAT→GAT) topADNA topoisomerase I
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base G (3/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GCGAAGGGGAAGCCATTGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGTTGAGCGTGAGCGTGAGATCCG  >  NC_002947/2441221‑2441515
                                                                                                                                                  |                                                                                                                                                    
gcgAAGGGGAAGCCATTGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAAt                                                                                                                                                    >  7:45842/1‑149 (MQ=255)
  gAAGGGGAAGCCATTGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGc                                                                                                                                                  <  6:24018/149‑1 (MQ=255)
     ggggAAGCCATTGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCca                                                                                                                                               <  5:92135/149‑1 (MQ=255)
        gAAGCCATTGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCcag                                                                                                                                              >  2:58041/1‑147 (MQ=255)
             cATTGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGgcgc                                                                                                                                       >  8:31749/1‑149 (MQ=255)
              aTTGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGGCATCGACCGGGTCAATGCCCAGCAGGcgcg                                                                                                                                      >  2:128286/1‑149 (MQ=255)
               ttGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCgcgg                                                                                                                                     <  6:63782/149‑1 (MQ=255)
                  ccTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGt                                                                                                                                  >  1:146843/1‑149 (MQ=255)
                      gCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCACGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTcc                                                                                                                              <  8:53143/149‑1 (MQ=255)
                        aCCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTc                                                                                                                            <  4:39424/149‑1 (MQ=255)
                         ccTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTc                                                                                                                            >  5:297112/1‑148 (MQ=255)
                             cgcgAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGAcc                                                                                                                        <  7:31517/148‑1 (MQ=255)
                                gAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGt                                                                                                                    >  5:239593/1‑149 (MQ=255)
                                 aGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTg                                                                                                                   <  3:39884/149‑1 (MQ=255)
                                   gCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGt                                                                                                                 >  5:10918/1‑149 (MQ=255)
                                      aTCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCgg                                                                                                              >  5:194626/1‑149 (MQ=255)
                                         gtgggTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTa                                                                                                           <  6:271166/147‑1 (MQ=255)
                                             gTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATg                                                                                                       <  7:31749/149‑1 (MQ=255)
                                                     acCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGcc                                                                                               <  2:146843/149‑1 (MQ=255)
                                                     acCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGcc                                                                                               >  8:269754/1‑149 (MQ=255)
                                                              aaGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCtgtg                                                                                      >  1:275971/1‑149 (MQ=255)
                                                                gcgcGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTggg                                                                                    <  6:194626/149‑1 (MQ=255)
                                                                     tggtgTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGcc                                                                                  >  5:29557/1‑146 (MQ=255)
                                                                      ggtgTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGa                                                                              <  7:269754/149‑1 (MQ=255)
                                                                       gtgtTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGGCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGa                                                                              <  6:160131/148‑1 (MQ=255)
                                                                        tgtTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACcgc                                                                                                                      <  3:7033/107‑1 (MQ=255)
                                                                        tgtTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACcgc                                                                                                                      >  4:7033/1‑107 (MQ=255)
                                                                          ttCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTc                                                                                                  <  6:162707/125‑1 (MQ=255)
                                                                          ttCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTc                                                                                                  >  5:162707/1‑125 (MQ=255)
                                                                          ttCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGc                                                                          >  1:26813/1‑149 (MQ=255)
                                                                             aaCGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGcc                                                                                                                                                 <  8:50855/75‑1 (MQ=255)
                                                                             aaCGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGcc                                                                                                                                                 >  7:50855/1‑75 (MQ=255)
                                                                             aaCGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCg                                                                       >  8:5283/1‑149 (MQ=255)
                                                                             aaCGAAATTACCAAGAAGGACATTAAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGACCAGCAGCCGCGTCGTTTCCTCGAACGCTTGGTCGGTTACATGGTTTCGCAACTGCTGTGCGCCAAGATTGccct                                                                       >  7:217917/1‑148 (MQ=255)
                                                                                   aTTACCGAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCt                                                                 <  1:14391/149‑1 (MQ=255)
                                                                                   aTTACCAAGAATGCCATTCAGGAATCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGAGCGGCTTTTCCTCGACCGCGTGGTCGTTTAAATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCt                                                                 <  8:217917/149‑1 (MQ=255)
                                                                                             agGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGGCATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTgg                                                        <  1:128286/148‑1 (MQ=255)
                                                                                                  aTTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCgtg                                                   >  7:218016/1‑148 (MQ=255)
                                                                                                  aTTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCgtg                                                   >  7:129983/1‑148 (MQ=255)
                                                                                                   ttCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTa                                                 <  7:5283/149‑1 (MQ=255)
                                                                                                   ttCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTa                                                 <  6:29557/149‑1 (MQ=255)
                                                                                                           gCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGt                                         <  8:128298/149‑1 (MQ=255)
                                                                                                           gCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGt                                         <  8:132545/149‑1 (MQ=255)
                                                                                                             cTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTATGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCgt                                                    >  3:293777/1‑136 (MQ=255)
                                                                                                             cTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTATGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCgt                                                    <  4:293777/136‑1 (MQ=255)
                                                                                                              ttCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCt                                                                                                                             <  5:211772/62‑1 (MQ=255)
                                                                                                              ttCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCt                                                                                                                             >  6:211772/1‑62 (MQ=255)
                                                                                                                   gCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCgg                                    <  2:254/146‑1 (MQ=255)
                                                                                                                         aGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCg                                           <  5:181597/133‑1 (MQ=255)
                                                                                                                         aGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCg                                           >  6:181597/1‑133 (MQ=255)
                                                                                                                           gCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGAt                                                                             <  7:23893/97‑1 (MQ=255)
                                                                                                                           gCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGAt                                                                             >  8:23893/1‑97 (MQ=255)
                                                                                                                           gCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGACCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCtggtg                         <  6:239593/149‑1 (MQ=255)
                                                                                                                              tgCTGGGCATTGACCCGGTCGATGTCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGACAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGtt                      >  7:287920/2‑149 (MQ=255)
                                                                                                                              aGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGtt                      <  2:26813/149‑1 (MQ=255)
                                                                                                                               gCTGGGCATTGACCCGGTCGATGTCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTTGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCtggtg                         <  8:24923/145‑1 (MQ=255)
                                                                                                                               gCTGGGCATTGACCCGGTCGATGTCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCtggtg                         >  7:24923/1‑145 (MQ=255)
                                                                                                                                cTGGGCATTGACCCGGTCGATGTCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGtt                      >  5:260368/1‑147 (MQ=255)
                                                                                                                                cTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGTtga                    <  8:129983/149‑1 (MQ=255)
                                                                                                                                cTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTATACAGTCGGTGGCGGTGAAGCTGGTGGTtga                    <  2:275971/149‑1 (MQ=255)
                                                                                                                                  ggACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGTtgagc                  >  7:63040/1‑149 (MQ=255)
                                                                                                                                   gACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGTtgagcg                 >  4:308849/1‑149 (MQ=255)
                                                                                                                                         gACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGTtgagcgtgagc            <  6:10918/148‑1 (MQ=255)
                                                                                                                                            ccggTCGATGTCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGTTGAgcgtgagcgtga        <  6:260368/147‑1 (MQ=255)
                                                                                                                                              ggTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCg                                                            <  2:305846/95‑1 (MQ=255)
                                                                                                                                              ggTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCg                                                            >  1:305846/1‑95 (MQ=255)
                                                                                                                                                  aaTGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGTTGAGCGTGAGCGTGAGATCcg  <  1:58041/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGAAGGGGAAGCCATTGCCTGGCACCTGCGCGAGGCCATCGGGGGTGACGACACCCGCTACAAGCGCGTGGTGTTCAACGAAATTACCAAGAAGGCCATTCAGGAAGCCTTCTCGCAGCCAGGCGAGCTGGACATCGACCGGGTCAATGCCCAGCAGGCGCGGCGTTTCCTCGACCGCGTGGTCGGTTACATGGTTTCGCCACTGCTGTGGGCCAAGATTGCCCGTGGTCTGTCCGCTGGCCGTGTACAGTCGGTGGCGGTGAAGCTGGTGGTTGAGCGTGAGCGTGAGATCCG  >  NC_002947/2441221‑2441515

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