Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,239,937 T→C 8.6% P142P (CCA→CCG)  ldh ← leucine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,239,9370TC 8.6% 184.1 / 5.8 69P142P (CCA→CCG) ldhleucine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (4/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTC  >  NC_002947/5239793‑5240080
                                                                                                                                                |                                                                                                                                               
gcCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGAACGCTACCCGCAACCCTCTCAGGTCATCGCTGCCCAGCCGGCCCGAGGTTGTGGCACGTATGCCCGCGAACACCCCCCTCGGCGCCTGTGGCgg                                                                                                                                            >  1:344824/1‑149 (MQ=255)
 cccACCGCATGCAATTGCAAAGCCACGGCCTAACCATCCTCGACGCTCCCCCGGACCGATACCAGCACGACCCTCCGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGAt                                                                                                                                           <  2:1081604/150‑1 (MQ=255)
  ccACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATg                                                                                                                                          <  1:204985/150‑1 (MQ=255)
   cACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATgg                                                                                                                                         <  1:325504/150‑1 (MQ=255)
   cACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCACATGCGGCGATgg                                                                                                                                         <  2:253504/150‑1 (MQ=255)
        cATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGc                                                                                                                                    >  1:838976/1‑150 (MQ=255)
              aTTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATg                                                                                                                              <  1:437156/150‑1 (MQ=255)
                  cTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgt                                                                                                                          <  1:781912/150‑1 (MQ=255)
                       ccAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTTCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGc                                                                                                                     <  1:685581/150‑1 (MQ=255)
                              aTAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCAtgtg                                                                                                                                                <  1:765925/116‑1 (MQ=255)
                              aTAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCAtgtg                                                                                                                                                >  2:765925/1‑116 (MQ=255)
                                       cctGCCCACCCGCTGGACCGCTACCC‑CGAGCCTCTCAGGTCACCGCTGCCCAGCAGGAACGAGGCTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGCCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTg                                                                                                    <  1:87197/148‑1 (MQ=255)
                                            ccAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgt                                                                                                >  1:1120794/1‑150 (MQ=255)
                                             cAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgtg                                                                                               >  1:381425/1‑150 (MQ=255)
                                               aCCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTggg                                                                                             <  2:439994/150‑1 (MQ=255)
                                                 cccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACAGCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCa                                                                                           >  1:747537/1‑150 (MQ=255)
                                                 cccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCa                                                                                           >  1:572956/1‑150 (MQ=255)
                                                  ccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                              >  1:34961/1‑146 (MQ=255)
                                                  ccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                              <  2:34961/146‑1 (MQ=255)
                                                  ccTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCaa                                                                                          <  1:654459/150‑1 (MQ=255)
                                                         cGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTc                                                                                   >  2:100304/1‑150 (MQ=255)
                                                           cTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCCAACCCTCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTCGGCAATGCAGTCCa                                                                                 <  1:426700/150‑1 (MQ=255)
                                                            tACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgtg                                                                                               >  2:761611/1‑135 (MQ=255)
                                                            tACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgtg                                                                                               <  1:761611/135‑1 (MQ=255)
                                                              cccGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTc                                                                                   >  1:253483/1‑145 (MQ=255)
                                                              cccGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTc                                                                                   <  2:253483/145‑1 (MQ=255)
                                                               ccGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTc                                                                             >  2:909826/1‑150 (MQ=255)
                                                                        tctcAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGAt                                                                                                                               >  2:885642/1‑91 (MQ=255)
                                                                        tctcAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGAt                                                                                                                               <  1:885642/91‑1 (MQ=255)
                                                                         ctcAGGTCATCGCTGCCCAGACGGAACGAGGGTGTGGGACGTATGCCGGCGAACACCCCCAGGGCCGCATGTGGCGAGGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                                   <  2:344824/150‑1 (MQ=255)
                                                                            aGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTg                                                                >  2:28230/1‑150 (MQ=255)
                                                                                           aGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGTAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCg                                                 >  2:1042611/1‑150 (MQ=255)
                                                                                                gAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGc                                                                                            >  1:915878/1‑102 (MQ=255)
                                                                                                gAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGc                                                                                            <  2:915878/102‑1 (MQ=255)
                                                                                                  aCGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGAtgatga                                          >  2:650109/1‑150 (MQ=255)
                                                                                                   cGAGTTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGTGTCGCCCGATGCAGTGTTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAATGTCGAGGTGCCACTGTCCACGGCGATGATGaa                                         <  2:1131250/150‑1 (MQ=255)
                                                                                                   cGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGaa                                         <  2:953408/150‑1 (MQ=255)
                                                                                                   cGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGaa                                         <  1:28230/150‑1 (MQ=255)
                                                                                                     aGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTAGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCACTGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgc                                       >  2:197247/1‑150 (MQ=255)
                                                                                                       gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                                   <  1:217632/120‑1 (MQ=255)
                                                                                                       gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                                   >  2:217632/1‑120 (MQ=255)
                                                                                                       gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcgc                                     >  1:784281/1‑150 (MQ=255)
                                                                                                       gTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcgc                                     >  2:13543/1‑150 (MQ=255)
                                                                                                           tgGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgt                                                                                                <  1:805069/87‑1 (MQ=255)
                                                                                                           tgGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgt                                                                                                >  2:805069/1‑87 (MQ=255)
                                                                                                             gCACGTATGCCCGCGAACACGCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGc                               >  2:100102/1‑150 (MQ=255)
                                                                                                                cGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCGTGCAAg                            <  1:100304/150‑1 (MQ=255)
                                                                                                                       ccGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGa                     >  2:844953/1‑150 (MQ=255)
                                                                                                                               acCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgtg                                                                                               >  1:894617/1‑68 (MQ=255)
                                                                                                                               acCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgtg                                                                                               <  2:894617/68‑1 (MQ=255)
                                                                                                                               acCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCACGTCCAAGGTCGAGCTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCgg             >  2:607324/1‑150 (MQ=255)
                                                                                                                                 ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtg                                                                                                                         <  2:140280/40‑1 (MQ=255)
                                                                                                                                 ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgtg                                                                                                                         >  1:140280/1‑40 (MQ=255)
                                                                                                                                 ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTg                                                                <  1:672094/97‑1 (MQ=255)
                                                                                                                                 ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTg                                                                >  2:672094/1‑97 (MQ=255)
                                                                                                                                 ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACt                                                           <  2:203366/102‑1 (MQ=255)
                                                                                                                                 ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACt                                                           >  1:203366/1‑102 (MQ=255)
                                                                                                                                  cccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGc                                                  >  2:297923/1‑110 (MQ=255)
                                                                                                                                  cccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGc                                                  <  1:297923/110‑1 (MQ=255)
                                                                                                                                  cccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCa          >  1:1104131/1‑150 (MQ=255)
                                                                                                                                   ccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCaa         >  2:23487/1‑150 (MQ=255)
                                                                                                                                    cATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCaaa        >  1:638814/1‑150 (MQ=255)
                                                                                                                                     aTCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAg       <  1:239174/150‑1 (MQ=255)
                                                                                                                                     aTCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCaa                             <  2:647221/128‑1 (MQ=255)
                                                                                                                                     aTCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCaa                             >  1:647221/1‑128 (MQ=255)
                                                                                                                                     aTCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAg       >  2:935458/1‑150 (MQ=255)
                                                                                                                                         ccgcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCt   <  1:100102/150‑1 (MQ=255)
                                                                                                                                          cgcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTc  >  1:462129/1‑150 (MQ=255)
                                                                                                                                          cgcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTc  <  2:838976/150‑1 (MQ=255)
                                                                                                                                          cgcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTc  <  1:533997/150‑1 (MQ=255)
                                                                                                                                          cgcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTc  >  1:889243/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTC  >  NC_002947/5239793‑5240080

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