Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,245,334 G→C 10.7% P65R (CCG→CGG)  pyrC ← dihydroorotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,245,3340GC10.7% 158.2 / 9.0 56P65R (CCG→CGG) pyrCdihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base C (2/4);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

ACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGAT  >  NC_002947/1245187‑1245481
                                                                                                                                                    |                                                                                                                                                    
aCGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCACCAGCGGCTCGAAGCGGCTGccggccgg                                                                                                                                                     >  6:240088/1‑149 (MQ=255)
  gCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggc                                                                                                                                                   >  2:201815/1‑149 (MQ=255)
   ccGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcg                                                                                                                                                  <  4:316240/149‑1 (MQ=255)
         tCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGTCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAg                                                                                                                                             >  1:272751/1‑149 (MQ=255)
         tCGGAGTTGGTGGGGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCCCGGCCGGCAGc                                                                                                                                            <  6:330082/149‑1 (MQ=255)
          cGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCGGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCCCGGCCGGCAGcc                                                                                                                                           <  1:318029/149‑1 (MQ=255)
             aGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcggg                                                                                                                                                >  1:317661/1‑141 (MQ=255)
             aGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcggg                                                                                                                                                <  2:317661/141‑1 (MQ=255)
              gTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGaa                                                                                                                                       <  5:60191/149‑1 (MQ=255)
                  gtggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcg                                                                                                                                   >  2:132832/1‑149 (MQ=255)
                       ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCgggcgggc                                                                                                                                               >  5:109428/1‑132 (MQ=255)
                       ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCgggcgggc                                                                                                                                               <  6:109428/132‑1 (MQ=255)
                       ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGc                                                                                                                                            >  3:296795/1‑135 (MQ=255)
                       ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGc                                                                                                                                            <  4:296795/135‑1 (MQ=255)
                        gcgcCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACg                                                                                                                             >  5:114607/1‑149 (MQ=255)
                             ggccggGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAgg                                                                                                                        >  3:311144/1‑149 (MQ=255)
                                     aCAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTa                                                                                                                          <  4:21152/139‑1 (MQ=255)
                                     aCAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTa                                                                                                                          >  3:21152/1‑139 (MQ=255)
                                       aaCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTc                                                                                                              >  6:21731/1‑149 (MQ=255)
                                       aaCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTc                                                                                                              >  8:293017/1‑149 (MQ=255)
                                            ggcggcGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAAt                                                                                                                                      >  7:113109/1‑120 (MQ=255)
                                            ggcggcGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAAt                                                                                                                                      <  8:113109/120‑1 (MQ=255)
                                               ggcgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACAt                                                                                                      >  5:167251/1‑149 (MQ=255)
                                                   tACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                                    >  7:296341/1‑115 (MQ=255)
                                                   tACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                                    <  8:296341/115‑1 (MQ=255)
                                                    acacGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgcg                                                                                                  >  1:5155/1‑148 (MQ=255)
                                                          aaCCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggc                                                                                          <  1:132832/149‑1 (MQ=255)
                                                                   gCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAg                                                                                  >  2:323013/1‑149 (MQ=255)
                                                                    ccTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAgg                                                                                 <  2:324116/149‑1 (MQ=255)
                                                                     cTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGt                                                                                >  7:244215/1‑149 (MQ=255)
                                                                      ttGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGtt                                                                               <  7:293017/149‑1 (MQ=255)
                                                                         gCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTggg                                                                            >  3:177226/1‑149 (MQ=255)
                                                                         gCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTgg                                                                             <  5:21731/148‑1 (MQ=255)
                                                                            gcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACa                                                                                                       >  6:235693/1‑120 (MQ=255)
                                                                            gcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACa                                                                                                       <  5:235693/120‑1 (MQ=255)
                                                                            gcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCat                                                                         <  2:5155/149‑1 (MQ=255)
                                                                             cgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatg                                                                        <  2:272751/149‑1 (MQ=255)
                                                                             cgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCACCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatg                                                                        <  5:240088/149‑1 (MQ=255)
                                                                              gcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTCCGCACGGGCGGCACGAGGTTGGGCatga                                                                       <  1:201815/149‑1 (MQ=255)
                                                                                 gACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGGGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatga                                                                    >  3:28904/1‑149 (MQ=255)
                                                                                     tCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTg                                                                              <  2:24013/135‑1 (MQ=255)
                                                                                     tCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTg                                                                              >  1:24013/1‑135 (MQ=255)
                                                                                      ccTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCa                                                               >  3:24832/1‑149 (MQ=255)
                                                                                          ggggCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCccggc                                                                                                                                                   >  8:133289/1‑62 (MQ=255)
                                                                                          ggggCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCccggc                                                                                                                                                   <  7:133289/62‑1 (MQ=255)
                                                                                                  tGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTa                                                   <  8:244215/149‑1 (MQ=255)
                                                                                                  tGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTa                                                   >  1:345293/1‑149 (MQ=255)
                                                                                                    cggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcg                                                 >  5:205077/1‑149 (MQ=255)
                                                                                                     ggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcgc                                                <  4:311144/149‑1 (MQ=255)
                                                                                                           tGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACg                                          >  8:244628/1‑149 (MQ=255)
                                                                                                            gAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatg                                                                     <  5:173811/121‑1 (MQ=255)
                                                                                                            gAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatg                                                                     >  6:173811/1‑121 (MQ=255)
                                                                                                                  cAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGTCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGcc                                     <  1:94162/148‑1 (MQ=255)
                                                                                                                     cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatg                                                                     <  4:68788/112‑1 (MQ=255)
                                                                                                                     cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatg                                                                     >  3:68788/1‑112 (MQ=255)
                                                                                                                     cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACgg                                >  4:260132/1‑149 (MQ=255)
                                                                                                                      aCCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCAGGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgt                               <  6:114607/149‑1 (MQ=255)
                                                                                                                       ccATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtg                              <  4:24832/149‑1 (MQ=255)
                                                                                                                          tCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgg                           <  6:205077/149‑1 (MQ=255)
                                                                                                                           cAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                            <  7:201787/147‑1 (MQ=255)
                                                                                                                           cAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                            >  8:201787/1‑147 (MQ=255)
                                                                                                                                    attctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg                 >  1:232799/6‑149 (MQ=255)
                                                                                                                                    attctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg                 <  2:232799/144‑1 (MQ=255)
                                                                                                                                         ggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACATCGCCAACGGTGTggg                          <  1:322852/135‑1 (MQ=255)
                                                                                                                                         ggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACATCGCCAACGGTGTggg                          >  2:322852/1‑135 (MQ=255)
                                                                                                                                                 gccggGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAg    <  3:260132/149‑1 (MQ=255)
                                                                                                                                                    gggcgggcAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAgat  <  4:105866/148‑1 (MQ=255)
                                                                                                                                                    gggcgggcAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAgat  >  3:105866/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGAT  >  NC_002947/1245187‑1245481

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: