Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,926,153 Δ1 bp 7.9% coding (382/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,1530A. 7.9% 280.0 / 16.6 89coding (382/1629 nt)peaAquinohemoprotein amine dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (37/45);  new base . (3/4);  total (40/49)
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.33e-01

GGCCAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGC  >  NC_002947/3926006‑3926301
                                                                                                                                                   |                                                                                                                                                     
ggCCAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAgc                                                                                                                                                     <  1:1789380/150‑1 (MQ=255)
 gCCAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAgcg                                                                                                                                                    >  2:1365243/1‑150 (MQ=255)
    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                                                                                                                                  <  1:669585/150‑1 (MQ=255)
                  ccAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGc                                                                                                                                  <  2:1071224/150‑1 (MQ=255)
                    aGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTg                                                                                                                                <  1:2230038/150‑1 (MQ=255)
                     ggTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCg                                                                                                                                         <  1:718044/140‑1 (MQ=255)
                     ggTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCg                                                                                                                                         >  2:718044/1‑140 (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                                                                                                                  <  2:2005213/150‑1 (MQ=255)
                                  cTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTg                                                                                                                  >  2:1184798/1‑150 (MQ=255)
                                  cTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTg                                                                                                                  >  1:773162/1‑150 (MQ=255)
                                   ttcgCCGGGCGCCGTTCCAAGGCGACACGCGCACCCGAGTGGCACCGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACACTATTGCGCACATACGCCATGGCGTAGCGCACGCCATCGGTTTCGCTCGGCGCCAGGCCCTGt                                                                                                                 <  2:338538/147‑1 (MQ=255)
                                   ttGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGgc                                                                                                                             <  2:2174775/138‑1 (MQ=255)
                                   ttGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGgc                                                                                                                             >  1:2174775/1‑138 (MQ=255)
                                   ttGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGt                                                                                                                 <  2:779498/150‑1 (MQ=255)
                                   ttGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGt                                                                                                                 >  2:2064419/1‑150 (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                                                                                                            >  1:348354/1‑121 (MQ=255)
                                                                     ccGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGCGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTc                                                                                                            <  2:348354/121‑1 (MQ=255)
                                                                                   ccACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATc                                                                 >  2:1963655/1‑150 (MQ=255)
                                                                                   ccACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATc                                                                 >  2:174643/1‑150 (MQ=255)
                                                                                    cacaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCt                                                                     <  2:87680/145‑1 (MQ=255)
                                                                                    cacaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCt                                                                     >  1:87680/1‑145 (MQ=255)
                                                                                    cacaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCt                                                                >  2:447306/1‑150 (MQ=255)
                                                                                    cacaGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCt                                                                >  2:309877/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                                                                                             >  1:701031/1‑95 (MQ=255)
                                                                                                              aaGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGAc                                                                                             <  2:701031/95‑1 (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                                                                              <  2:1415866/80‑1 (MQ=255)
                                                                                                                                          cgaTGGCGTGCGCTACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGAtc                                                                              >  1:1415866/3‑82 (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)
                                                                                                                                                   |                                                                                                                                                     
GGCCAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGC  >  NC_002947/3926006‑3926301

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