Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,245,336 G→C 7.6% A66A (GCC→GCG)  pyrC ← dihydroorotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,245,3360GC 7.6% 202.8 / 7.4 66A66A (GCC→GCG) pyrCdihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base C (3/2);  total (34/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 = 7.80e-01

GTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAG  >  NC_002947/1245190‑1245485
                                                                                                                                                  |                                                                                                                                                     
gTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTgccggcc                                                                                                                                                    <  2:640290/150‑1 (MQ=255)
   aCGCCGGAATCGGAGTTGGTGGTGGCCCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcggg                                                                                                                                                 >  1:1404907/1‑150 (MQ=255)
     gCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcg                                                                                                                                               <  1:620633/150‑1 (MQ=255)
     gCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGCCTGCCGGCcgggcg                                                                                                                                               >  2:283602/1‑150 (MQ=255)
         gAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCa                                                                                                                                           <  1:522271/150‑1 (MQ=255)
           aTCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGc                                                                                                                                         <  1:1295381/150‑1 (MQ=255)
           aTCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCCCGGCCGGCAGc                                                                                                                                         <  2:1073179/150‑1 (MQ=255)
                 gTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAAt                                                                                                                                   >  2:123408/1‑150 (MQ=255)
                     gtggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcgc                                                                                                                               <  2:1404907/150‑1 (MQ=255)
                     gtggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcgc                                                                                                                               >  1:1079640/1‑150 (MQ=255)
                       ggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAACCCGCGGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTc                                                                                                                             <  1:524797/150‑1 (MQ=255)
                                 gccggGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGgcg                                                                                                                   >  1:1060949/1‑150 (MQ=255)
                                   cggGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGcc                                                                                                                 >  2:1309839/1‑150 (MQ=255)
                                    gggTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCg                                                                                                                >  1:926869/1‑150 (MQ=255)
                                          aaCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCa                                                                                                                                           <  2:316247/117‑1 (MQ=255)
                                          aaCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCa                                                                                                                                           >  1:316247/1‑117 (MQ=255)
                                          aaCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCg                                                                                                          <  1:831692/150‑1 (MQ=255)
                                            cTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGCTGGACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGt                                                                                                        <  1:89219/150‑1 (MQ=255)
                                            cTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGt                                                                                                        <  2:926869/150‑1 (MQ=255)
                                            cTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGt                                                                                                        <  1:84200/150‑1 (MQ=255)
                                            cTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGCGGCCACAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGt                                                                                                        >  1:374368/1‑150 (MQ=255)
                                            cTTGGAGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGt                                                                                                        <  1:1309839/150‑1 (MQ=255)
                                                gcggcgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACa                                                                                                    <  1:403703/150‑1 (MQ=255)
                                                   gcgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCa                                                                                                                                           >  1:222934/1‑108 (MQ=255)
                                                   gcgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCa                                                                                                                                           <  2:222934/108‑1 (MQ=255)
                                                      tACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgc                                                                                                >  1:47497/1‑148 (MQ=255)
                                                      tACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgc                                                                                                <  2:47497/148‑1 (MQ=255)
                                                           gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACggg                                                                                         >  1:297197/1‑150 (MQ=255)
                                                              aCCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggcgg                                                                                      >  2:880320/1‑150 (MQ=255)
                                                                       cTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGtt                                                                             <  1:880320/150‑1 (MQ=255)
                                                                       cTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGtt                                                                             >  1:725830/1‑150 (MQ=255)
                                                                        ttGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTg                                                                            >  1:1228671/1‑150 (MQ=255)
                                                                         tGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTgg                                                                           <  1:659639/150‑1 (MQ=255)
                                                                          ggCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTggg                                                                          <  2:528288/150‑1 (MQ=255)
                                                                           gCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGc                                                                         >  1:1405654/1‑150 (MQ=255)
                                                                               cgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatga                                                                     <  2:562844/150‑1 (MQ=255)
                                                                                       tCCTCGGGGCTGGTGCGGTAGGTGAGATACAGCACCATCAGCGGCTTGAAGCGGCTGCCGGCCGGGCGGGCAGCGAGAATGCGCTCACGGTAGGCGCCGGCGTCAGTGACATTGCGCACGGGCGGCACGACGTTGGGCATGATGATGGCa                                                             >  2:63110/1‑150 (MQ=255)
                                                                                         cTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACg                                                           >  1:295552/1‑150 (MQ=255)
                                                                                          tCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatga                                                                  >  1:1141410/1‑142 (MQ=255)
                                                                                          tCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatga                                                                  <  2:1141410/142‑1 (MQ=255)
                                                                                               gCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGaa                                                     <  1:1375561/150‑1 (MQ=255)
                                                                                                cTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAg                                                    <  2:543023/150‑1 (MQ=255)
                                                                                                      cggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcgc                                              <  1:504695/150‑1 (MQ=255)
                                                                                                         tcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGcc                                           >  2:501280/1‑150 (MQ=255)
                                                                                                                 gTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                                 <  1:379590/56‑1 (MQ=255)
                                                                                                                 gTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                                 >  2:379590/1‑56 (MQ=255)
                                                                                                                 gTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGcc                                   >  1:508844/1‑150 (MQ=255)
                                                                                                                 gTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGcc                                   >  2:591618/1‑150 (MQ=255)
                                                                                                                   aCAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgc                                                                                                <  2:530842/87‑1 (MQ=255)
                                                                                                                   aCAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgc                                                                                                >  1:530842/1‑87 (MQ=255)
                                                                                                                    cAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAAc                                >  2:1335764/1‑150 (MQ=255)
                                                                                                                     aGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTa                                                                                                                       >  2:492923/1‑62 (MQ=255)
                                                                                                                     aGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTa                                                                                                                       <  1:492923/62‑1 (MQ=255)
                                                                                                                     aGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACg                               <  2:725830/150‑1 (MQ=255)
                                                                                                                      gCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACgg                              >  2:587676/1‑150 (MQ=255)
                                                                                                                        aCCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtg                            >  1:1412741/1‑150 (MQ=255)
                                                                                                                          cATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                          <  2:1405654/150‑1 (MQ=255)
                                                                                                                           aTCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgg                         <  1:123408/150‑1 (MQ=255)
                                                                                                                              aGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTTGCCAACGGTGTGGGGc                      <  1:283602/150‑1 (MQ=255)
                                                                                                                               gCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCa                     >  2:914664/1‑150 (MQ=255)
                                                                                                                                 ggCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAgg                   >  1:789545/1‑150 (MQ=255)
                                                                                                                                  gCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGTTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGa                  <  2:824160/150‑1 (MQ=255)
                                                                                                                                       ttctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGc                                                                         >  2:389536/5‑90 (MQ=38)
                                                                                                                                       ttctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGc                                                                         <  1:389536/86‑1 (MQ=38)
                                                                                                                                          tggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCa          >  1:1228004/2‑150 (MQ=255)
                                                                                                                                          tggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCCCACGGTCGGCGCCGGCCACGGTCCCATTCCGCAGGAGGAGCACTGGGGTGGTCGTGAGGAGGTCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCa          >  2:161721/2‑150 (MQ=255)
                                                                                                                                           ggCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCAt         >  2:1109147/1‑150 (MQ=255)
                                                                                                                                                  ggccggGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAg  >  2:40272/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAG  >  NC_002947/1245190‑1245485

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