Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,926,157:1 +T 8.0% coding (378/1629 nt) peaA ← quinohemoprotein amine dehydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,926,1571.T 8.0% 330.6 / 21.0 87coding (378/1629 nt)peaAquinohemoprotein amine dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base . (36/44);  new base T (3/4);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTG  >  NC_002947/3926010‑3926306
                                                                                                                                                    |                                                                                                                                                     
aaGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCAc                                                                                                                                                     <  2:1927457/150‑1 (MQ=255)
        aGTGGACCAGGTGTTCCCAGCCTTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtgga                                                                                                                                             <  2:1344806/150‑4 (MQ=255)
        aGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtgga                                                                                                                                             <  1:585993/150‑4 (MQ=255)
        aGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCAATGGCGTGCGCTACGCCAtgga                                                                                                                                             <  2:2213502/150‑4 (MQ=255)
           ggACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGttt                                                                                                                                          <  2:1632614/150‑1 (MQ=255)
              ccAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGc                                                                                                                                       <  2:1071224/150‑1 (MQ=255)
              ccAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGc                                                                                                                                       <  1:669585/150‑1 (MQ=255)
                aGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTg                                                                                                                                     <  1:2230038/150‑1 (MQ=255)
                 ggTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCg                                                                                                                                              >  2:718044/1‑140 (MQ=255)
                 ggTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCg                                                                                                                                              <  1:718044/140‑1 (MQ=255)
                 ggTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTgg                                                                                                                                    <  2:1103841/150‑1 (MQ=255)
                      tCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTGGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGcc                                                                                                                               <  1:1355925/150‑1 (MQ=255)
                        ccACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAg                                                                                                                             >  1:1576496/1‑150 (MQ=255)
                        ccACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAg                                                                                                                             >  1:779498/1‑150 (MQ=255)
                        ccACTCCTTGGCCGGGCGCCGTTGCAAGGAGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAg                                                                                                                             >  1:487920/1‑150 (MQ=255)
                            tCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTTCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGccc                                                                                                                         >  2:694871/1‑150 (MQ=255)
                              cTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTg                                                                                                                       >  1:773162/1‑150 (MQ=255)
                              cTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTg                                                                                                                       >  2:1184798/1‑150 (MQ=255)
                              cTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTg                                                                                                                       <  2:2005213/150‑1 (MQ=255)
                               ttcgCCGGGCGCCGTTCCAAGGCGACACGCGCACCCGAGTGGCACCGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACACTATTGCGCACATACGCCATGGCGTAGCGCACGCCATCGGTTTCGCTCGGCGCCAGGCCCTGt                                                                                                                      <  2:338538/147‑1 (MQ=255)
                               ttGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGgc                                                                                                                                  >  1:2174775/1‑138 (MQ=255)
                               ttGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGgc                                                                                                                                  <  2:2174775/138‑1 (MQ=255)
                               ttGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGt                                                                                                                      >  2:2064419/1‑150 (MQ=255)
                               ttGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGt                                                                                                                      <  2:779498/150‑1 (MQ=255)
                                 ggCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGttt                                                                                                                    >  1:1600893/1‑150 (MQ=255)
                                    cGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTc                                                                                                                 <  1:2214221/150‑1 (MQ=255)
                                     gggCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCg                                                                                                                <  2:2119186/150‑1 (MQ=255)
                                       gcgcCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGc                                                                                                              <  1:2105085/150‑1 (MQ=255)
                                         gcCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCCTCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCa                                                                                                            >  2:172536/1‑150 (MQ=255)
                                         gcCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCa                                                                                                            >  1:364738/1‑150 (MQ=255)
                                         gcCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCa                                                                                                            <  1:1184798/150‑1 (MQ=255)
                                           cGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACATGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAgg                                                                                                          >  1:789792/1‑150 (MQ=255)
                                           cGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAgg                                                                                                          >  2:1446704/1‑150 (MQ=255)
                                           cGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAgg                                                                                                          >  2:1465623/1‑150 (MQ=255)
                                            gTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGt                                                                                                         <  2:773162/150‑1 (MQ=255)
                                             ttGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTAGCCGCACCGGCTGTCGCAATGCTCGACAGTATTGAGGTGACGTTCCAGTGCTTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTa                                                                                                        <  2:880616/150‑1 (MQ=255)
                                              tGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACACTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAt                                                                                                       <  2:1344968/150‑1 (MQ=255)
                                                 aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTg                                                                                                    >  2:885308/1‑150 (MQ=255)
                                                  aGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGa                                                                                                   <  2:1857020/150‑1 (MQ=255)
                                                       acacGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGc                                                                                              <  1:1446704/150‑1 (MQ=255)
                                                         acGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCTAGGTATTTGACCAGCGt                                                                                            <  1:694871/150‑1 (MQ=255)
                                                           gggCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGc                                                                                          <  2:487920/150‑1 (MQ=255)
                                                              cACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGcgac                                                                                       >  2:18335/1‑150 (MQ=255)
                                                                 ccGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGCGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTc                                                                                                                 <  2:348354/121‑1 (MQ=255)
                                                                 ccGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGCGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTc                                                                                                                 >  1:348354/1‑121 (MQ=255)
                                                                               ccACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATc                                                                      >  2:1963655/1‑150 (MQ=255)
                                                                               ccACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATc                                                                      >  2:174643/1‑150 (MQ=255)
                                                                                cacaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCt                                                                          >  1:87680/1‑145 (MQ=255)
                                                                                cacaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCt                                                                          <  2:87680/145‑1 (MQ=255)
                                                                                cacaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCt                                                                     >  2:309877/1‑150 (MQ=255)
                                                                                cacaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCt                                                                     >  2:447306/1‑150 (MQ=255)
                                                                                 acaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTg                                                                    >  2:1177192/1‑150 (MQ=255)
                                                                                   aGGTCTCGCTGAGCTGGGTGTCGACGTGCTCGACAGTATTGAGCCGACGCTCCCTGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTTTGCCAGGTATTTGACCAGCGTGCGACGAGCGTCATCGCTGGTGTGaa                                                                  >  1:273362/1‑150 (MQ=255)
                                                                                   aGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTc                                                                                                                 <  2:845853/103‑1 (MQ=255)
                                                                                   aGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTc                                                                                                                 >  1:845853/1‑103 (MQ=255)
                                                                                   aGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGaa                                                                  <  1:1177192/150‑1 (MQ=255)
                                                                                        tcGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCa                                                             >  2:2353700/1‑150 (MQ=255)
                                                                                         cGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCAt                                                            <  1:1080362/150‑1 (MQ=255)
                                                                                           cTGAGCTGGGTGTCGAAGTGCTCGATAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGCCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTCCGACGATCGTCGTCGCTGATCTGAAGGCCATGc                                                          <  2:2081832/150‑1 (MQ=255)
                                                                                           cTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGc                                                          <  1:1465623/150‑1 (MQ=255)
                                                                                              aGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATc                                                       >  1:2015213/1‑150 (MQ=255)
                                                                                                cTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCAc                                                     <  1:309877/150‑1 (MQ=255)
                                                                                                 tgtgtgcCGAAGTGCTCGACAGAATTGAGCCGACGCTCCATTGCGTAGCTCGCGCCATCGGTTTCTATGGGCGCCAGGCCCTGTCTGTCGTCCAGCTCTTTGACCGGCGTATGACGATCATCGTCGCTGATCTGAAGGCCATGCATCAcc                                                    >  2:1471132/4‑150 (MQ=255)
                                                                                                   ggTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTg                                                  <  1:1922082/150‑1 (MQ=255)
                                                                                                          aaGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGAc                                                                                                  <  2:701031/95‑1 (MQ=255)
                                                                                                          aaGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGAc                                                                                                  >  1:701031/1‑95 (MQ=255)
                                                                                                           aGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCggg                                          <  1:2064419/150‑1 (MQ=255)
                                                                                                             tGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCg                                        >  1:833637/1‑150 (MQ=255)
                                                                                                               cTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCGCCCGCATGCAGGTGAt                                      <  1:174643/150‑1 (MQ=255)
                                                                                                               cTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGAt                                      <  1:2089003/150‑1 (MQ=255)
                                                                                                                tCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATg                                     <  1:2103063/150‑1 (MQ=255)
                                                                                                                 cGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGc                                    >  2:2176788/1‑150 (MQ=255)
                                                                                                                 cGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGc                                    >  1:1369844/1‑150 (MQ=255)
                                                                                                                  gACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGAt                                      >  2:1427346/1‑147 (MQ=255)
                                                                                                                  gACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGAt                                      <  1:1427346/147‑1 (MQ=255)
                                                                                                                  gACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCt                                   >  1:2025598/1‑150 (MQ=255)
                                                                                                                      gTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCtcatc                               <  2:1369844/150‑1 (MQ=255)
                                                                                                                       tATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCtcatca                              <  1:2176788/150‑1 (MQ=255)
                                                                                                                            aGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATcagccag                         <  2:789792/150‑1 (MQ=255)
                                                                                                                                   gCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCg                  <  1:110894/150‑1 (MQ=255)
                                                                                                                                      cgaTGGCGTGCGCTACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGGCAGGGATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgtg               >  2:973463/3‑150 (MQ=255)
                                                                                                                                      cgaTGGCGTGCGCTACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGAtc                                                                                   >  1:1415866/3‑82 (MQ=255)
                                                                                                                                      cgaTGGCGTGCGCTACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGAtc                                                                                   <  2:1415866/80‑1 (MQ=255)
                                                                                                                                       gaTGGCGTGCGCTACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgtgt              >  1:1414405/2‑150 (MQ=255)
                                                                                                                                       cATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCAt                                               <  2:1095090/117‑1 (MQ=255)
                                                                                                                                       cATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCAt                                               >  1:1095090/1‑117 (MQ=255)
                                                                                                                                              tAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTgcgc       <  2:364738/150‑1 (MQ=255)
                                                                                                                                                  gcACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCtggt   >  2:816137/1‑150 (MQ=255)
                                                                                                                                                   cACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCtggtg  <  1:885308/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTG  >  NC_002947/3926010‑3926306

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: