Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 648,883 C→T 5.3% A266A (GCG→GCA)  accC ← acetyl‑CoA carboxylase biotin carboxylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947648,8830CT 5.3% 222.9 / 2.4 75A266A (GCG→GCA) accCacetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base T (2/2);  total (39/36)
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.13e-01

GCATCTCCTTGACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGTTGCCCTGGCCGTCGGACAGTAC  >  NC_002947/648735‑649032
                                                                                                                                                    |                                                                                                                                                     
gCATCTCCTTGACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACg                                                                                                                                                      >  2:933139/1‑150 (MQ=255)
       cTTGACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCTCTGCCCCCCTGTGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGAc                                                                                                                                               <  2:456665/150‑1 (MQ=255)
       cTTGACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGACCTCGAAAGTACCGGCACCACCGTAGCCGATCTCGATGCACGCATCGAc                                                                                                                                               <  2:1485726/150‑1 (MQ=255)
       cTTGACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGAc                                                                                                                                               <  2:230479/150‑1 (MQ=255)
         tGACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTAAACGAATCACGGTAGCCGATCTCGATGCGTGCATCGagat                                                                                                                                             <  1:1206269/150‑4 (MQ=255)
          gACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCa                                                                                                                                            <  1:454293/150‑1 (MQ=255)
          gACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGCCAGCACGGTAGCCGATCTCGATGCGTGCATCGagatc                                                                                                                                            >  2:1206269/1‑146 (MQ=255)
                   gATACCAGTGACAATATCCGACTCTGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGaa                                                                                                                                   <  2:785145/150‑1 (MQ=255)
                    aTACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAac                                                                                                                                  >  1:1000041/1‑150 (MQ=255)
                      aCCAGTGACCATCACCGACCACGGGTGCTCAGCCTGCATGCGGGTTTTCATCTCTATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAacac                                                                                                                                <  2:1000041/150‑1 (MQ=255)
                           tGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCt                                                                                                                           >  1:1326329/1‑150 (MQ=255)
                           tGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCt                                                                                                                           >  1:914257/1‑150 (MQ=255)
                                aTCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACgg                                                                                                                      <  1:838982/150‑1 (MQ=255)
                                 tctcCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCAt                                                                                                                                                   >  2:1543307/1‑120 (MQ=255)
                                 tctcCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCAt                                                                                                                                                   <  1:1543307/120‑1 (MQ=255)
                                             ggTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGa                                                                                                         >  2:546803/1‑150 (MQ=255)
                                                 cTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGcc                                                                                                     >  1:1220513/1‑150 (MQ=255)
                                                    aaCCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCa                                                                                                                                            <  1:1164434/108‑1 (MQ=255)
                                                    aaCCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCa                                                                                                                                            >  2:1164434/1‑108 (MQ=255)
                                                     aCCTGCACGCGGGTGTTCATCTCGATGAAGTGGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCggg                                                                                                 >  2:1010562/1‑150 (MQ=255)
                                                             gcgGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGc                                                                                         <  1:1010562/150‑1 (MQ=255)
                                                                agtgtTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGTAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGcttg                                                                                      <  2:341957/149‑2 (MQ=255)
                                                                   gtTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTg                                                                                             <  2:1012016/140‑1 (MQ=255)
                                                                   gtTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTg                                                                                             >  1:1012016/1‑140 (MQ=255)
                                                                    ttCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGAc                                                                                                                                               >  1:24493/1‑89 (MQ=255)
                                                                    ttCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGAc                                                                                                                                               <  2:24493/89‑1 (MQ=255)
                                                                     tCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCa                                                                                 <  1:570021/150‑1 (MQ=255)
                                                                     tCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCa                                                                                 <  2:1326329/150‑1 (MQ=255)
                                                                      cATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAg                                                                                >  1:1556084/1‑150 (MQ=255)
                                                                       aTCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCt                                                                                                                                       >  1:1016592/1‑94 (MQ=255)
                                                                       aTCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCt                                                                                                                                       <  2:1016592/94‑1 (MQ=255)
                                                                        tctcGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATc                                                                                                                                                  >  2:583245/1‑82 (MQ=255)
                                                                        tctcGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATc                                                                                                                                                  <  1:583245/82‑1 (MQ=255)
                                                                        tctcGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTa                                                                              >  2:618556/1‑150 (MQ=255)
                                                                         ctcGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCCCGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTATCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTAc                                                                             >  1:1341193/1‑150 (MQ=255)
                                                                         ctcGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAaca                                                                                                                                 <  1:326041/98‑1 (MQ=255)
                                                                         ctcGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAaca                                                                                                                                 >  2:326041/1‑98 (MQ=255)
                                                                            gATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCa                                                                                 <  2:280117/143‑1 (MQ=255)
                                                                            gATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCa                                                                                 >  1:280117/1‑143 (MQ=255)
                                                                            gATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACttt                                                                          >  1:118199/1‑150 (MQ=255)
                                                                               gAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTg                                                                       <  1:618556/150‑1 (MQ=255)
                                                                               gAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTg                                                                       <  1:546803/150‑1 (MQ=255)
                                                                                     gAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAaca                                                                                                                                 >  2:588946/1‑86 (MQ=255)
                                                                                     gAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAaca                                                                                                                                 <  1:588946/86‑1 (MQ=255)
                                                                                          ggCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCttcttc                                                                                   >  2:853075/1‑127 (MQ=255)
                                                                                          ggCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCttcttc                                                                                   <  1:853075/127‑1 (MQ=255)
                                                                                           gCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATTTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTg                                                           >  1:1054652/1‑150 (MQ=255)
                                                                                           gCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTg                                                           >  2:256867/1‑150 (MQ=255)
                                                                                              gTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCAt                                                                                                            <  1:1012021/98‑1 (MQ=255)
                                                                                              gTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCAt                                                                                                            >  2:1012021/1‑98 (MQ=255)
                                                                                               ttttCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGc                                                       <  1:824042/150‑1 (MQ=255)
                                                                                                    gTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGc                                                                                                                    >  2:1299980/1‑84 (MQ=255)
                                                                                                    gTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGc                                                                                                                    <  1:1299980/84‑1 (MQ=255)
                                                                                                    gTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCa                                                  >  2:530248/1‑150 (MQ=255)
                                                                                                             cTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTc                                         >  1:1457719/1‑150 (MQ=255)
                                                                                                              tCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGc                                                                                               >  1:881370/1‑95 (MQ=255)
                                                                                                              tCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGc                                                                                               <  2:881370/95‑1 (MQ=255)
                                                                                                              tCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTc                                                                         <  2:300221/117‑1 (MQ=255)
                                                                                                              tCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTc                                                                         >  1:300221/1‑117 (MQ=255)
                                                                                                              tCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCg                                        >  1:1560474/1‑150 (MQ=255)
                                                                                                              gcGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCg                                        >  2:1124681/2‑150 (MQ=255)
                                                                                                                 aaaGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGccc                                     <  1:1124681/150‑1 (MQ=255)
                                                                                                                  aaGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCa                                    <  2:1556084/150‑1 (MQ=255)
                                                                                                                  aaGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCa                                    >  1:91719/1‑150 (MQ=255)
                                                                                                                   aGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAg                                   <  2:914257/150‑1 (MQ=255)
                                                                                                                         ggCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAAt                             <  1:933139/150‑1 (MQ=255)
                                                                                                                             ccacGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTg                                                                                                                                     <  1:1207607/42‑1 (MQ=255)
                                                                                                                             ccacGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTg                                                                                                                                     >  2:1207607/1‑42 (MQ=255)
                                                                                                                              cacGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGt                        >  1:727945/1‑150 (MQ=255)
                                                                                                                                      ccGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGTTGCCCTgg                <  2:1457719/150‑1 (MQ=255)
                                                                                                                                           gtcGATGCGTGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGTTGCCCt                  <  1:1398144/142‑1 (MQ=255)
                                                                                                                                           gtcGATGCGTGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGTTGCCCt                  >  2:1398144/2‑143 (MQ=255)
                                                                                                                                                 gCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGTTGCCCTGGCCGTCGGACAg     >  2:1619386/1‑150 (MQ=255)
                                                                                                                                                   aCGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCa                                                  >  2:919049/1‑103 (MQ=255)
                                                                                                                                                   aCGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGATCGCTGCGGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGTTGCCCTGGCCGTCGGACAGTa   >  1:1350483/1‑150 (MQ=255)
                                                                                                                                                    cGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGTTGCCCTGGCCGTCGGACAGTAc  <  2:727945/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCATCTCCTTGACGATGTCGATACCAGTGACCATCTCCGACACCGGGTGCTCAACCTGCACGCGGGTGTTCATCTCGATGAAGTAGAAGCGGCCGTTTTCGTACAGGAACTCGAAAGTACCGGCACCACGGTAGCCGATCTCGATGCACGCATCGACGCAGCGCTTGAACACTTCCTGACGGGCTTTCTCATCGATGCCCGGGGCTGGTGCTTCTTCCAGTACTTTCTGGTGACGGCGCTGCAGCGAGCAGTCGCGGTCGCCCAGGTGAATGGCGTTGCCCTGGCCGTCGGACAGTAC  >  NC_002947/648735‑649032

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