Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,039,273 G→C 6.1% P42R (CCT→CGT)  PP_0897 ← class 1 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,039,2730GC 6.1% 259.1 / 7.2 82P42R (CCT→CGT) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (41/36);  new base C (2/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

ACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGG  >  NC_002947/1039126‑1039421
                                                                                                                                                    |                                                                                                                                                    
aCGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGc                                                                                                                                                    <  2:699634/150‑1 (MQ=255)
  gCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGa                                                                                                                                                  >  2:1425593/1‑150 (MQ=255)
  gCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGa                                                                                                                                                  <  2:824583/150‑1 (MQ=255)
  gCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGa                                                                                                                                                  >  1:1026888/1‑150 (MQ=255)
  gCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGa                                                                                                                                                  <  1:1109484/150‑1 (MQ=255)
   cTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGAt                                                                                                                                                 <  2:1794946/150‑1 (MQ=255)
   cTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCTGCACGCGAt                                                                                                                                                 <  2:860967/150‑1 (MQ=255)
    tCATGGTGGCGCCATCCCAGCGGACGTCTAGGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCTGCACGCGAtt                                                                                                                                                <  1:1467004/150‑1 (MQ=255)
    tCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGAtt                                                                                                                                                >  1:453015/1‑150 (MQ=255)
    tCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGAtt                                                                                                                                                >  1:1128945/1‑150 (MQ=255)
     cATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTc                                                                                                                                               <  1:21607/150‑1 (MQ=255)
              gcCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAg                                                                                                                                      >  2:1636822/1‑150 (MQ=255)
               ccATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGa                                                                                                                                     >  1:1126910/1‑150 (MQ=255)
               ccATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGa                                                                                                                                     <  1:773331/150‑1 (MQ=255)
               ccATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGa                                                                                                                                     <  2:1298513/150‑1 (MQ=255)
                 aTCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGAGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATa                                                                                                                                   >  1:118936/1‑150 (MQ=255)
                 aTCCCAGCGCACGTCCATGCCCACACGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATa                                                                                                                                   >  1:1485165/1‑150 (MQ=255)
                  tCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAg                                                                                                                                  >  2:515648/1‑150 (MQ=255)
                    ccAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGc                                                                                                                                >  1:485234/1‑150 (MQ=255)
                         gcACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                                           <  2:485234/150‑1 (MQ=255)
                         gcACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                                           >  1:413983/1‑150 (MQ=255)
                         ccACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                                           <  1:1118652/149‑1 (MQ=255)
                           aCGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGcagcag                                                                                                                                                      >  1:667890/1‑121 (MQ=255)
                           aCGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGcagcag                                                                                                                                                      <  2:667890/121‑1 (MQ=255)
                             gTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGTACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCAt                                                                                                                       >  2:595546/1‑150 (MQ=255)
                             gTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCAt                                                                                                                       >  1:1047596/1‑150 (MQ=255)
                              tCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATg                                                                                                                      <  2:1026888/150‑1 (MQ=255)
                              ccgaagCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATg                                                                                                                      <  2:1151826/145‑1 (MQ=255)
                                 aTGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                                           >  2:1432512/1‑142 (MQ=255)
                                 aTGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                                           <  1:1432512/142‑1 (MQ=255)
                                      cACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGAt                                                                                                              <  2:578583/150‑1 (MQ=255)
                                       aCGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGTCCTGGATg                                                                                                             >  2:1046633/1‑150 (MQ=255)
                                       aCGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGACGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCATGGATg                                                                                                             >  2:717627/1‑150 (MQ=255)
                                           ggATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGAACTACATAGTCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAAt                                                                                                         >  1:593614/1‑150 (MQ=255)
                                            gATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATc                                                                                                        <  2:1128945/150‑1 (MQ=255)
                                                 acacGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACgg                                                                                                   >  2:377667/1‑150 (MQ=255)
                                                           gATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGAGCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATTGCCTGGATGAAATCGCGGGGGTGgtagta                                                                                         >  2:894067/1‑150 (MQ=255)
                                                            aTACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTAc                                                                                        <  2:1332594/150‑1 (MQ=255)
                                                                cGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACgaga                                                                                    >  2:1162651/1‑150 (MQ=255)
                                                                cGGTGTCCTGGCAGAGTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACgaga                                                                                    <  1:1523505/150‑1 (MQ=255)
                                                                   tgtCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGa                                                                                 >  1:1792835/1‑150 (MQ=255)
                                                                        tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCa                                                                                                                                           <  1:88706/87‑1 (MQ=255)
                                                                        tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCa                                                                                                                                           >  2:88706/1‑87 (MQ=255)
                                                                        tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGc                                                                            >  2:24916/1‑150 (MQ=255)
                                                                        tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGAAGGGGTGGTAGGATGAGATGAATTga                                                                            >  1:391018/1‑149 (MQ=255)
                                                                         ggCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCa                                                                           <  1:1425593/150‑1 (MQ=255)
                                                                           cAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATGGGTGGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAgg                                                                         <  2:57703/150‑1 (MQ=255)
                                                                               ttGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTg                                                                             >  2:262354/1‑142 (MQ=255)
                                                                               ttGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTg                                                                             <  1:262354/142‑1 (MQ=255)
                                                                                  gACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCgg                                                                  <  1:1162651/150‑1 (MQ=255)
                                                                                   aCGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGc                                                                 >  2:807432/1‑150 (MQ=255)
                                                                                     gATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                                               <  1:1052651/102‑1 (MQ=255)
                                                                                     gATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAAGCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGa                                                                                                               >  2:1052651/1‑102 (MQ=255)
                                                                                       tGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACg                                                             >  1:1306619/1‑150 (MQ=255)
                                                                                         gCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGct                                                           >  1:804365/1‑150 (MQ=255)
                                                                                           cGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGctct                                                         >  2:839553/1‑150 (MQ=255)
                                                                                           cGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGCGGCGTTGCGAGATGAATTGCAGGGCGTCGGCGACGctct                                                         <  1:1892929/150‑1 (MQ=255)
                                                                                             gTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCGCTGa                                                       >  1:1497433/1‑150 (MQ=255)
                                                                                                gcgcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATc                                                    <  2:413983/150‑1 (MQ=255)
                                                                                                gcgcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATc                                                    <  2:1126910/150‑1 (MQ=255)
                                                                                                 cgcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCa                                                   <  1:807432/150‑1 (MQ=255)
                                                                                                  gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                                           <  1:1100055/77‑1 (MQ=255)
                                                                                                  gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGt                                                                                                                           >  2:1100055/1‑77 (MQ=255)
                                                                                                  gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCa                                                                           <  1:265806/125‑1 (MQ=255)
                                                                                                  gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCa                                                                           >  2:265806/1‑125 (MQ=255)
                                                                                                  gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                                  <  1:24916/150‑1 (MQ=255)
                                                                                                  gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                                  <  2:118936/150‑1 (MQ=255)
                                                                                                          gcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTc                                                                                                                                               <  1:304486/49‑1 (MQ=255)
                                                                                                          gcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTc                                                                                                                                               >  2:304486/1‑49 (MQ=255)
                                                                                                          gcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTctt                                          >  2:1473937/1‑150 (MQ=255)
                                                                                                            gAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgc                                        <  2:1792835/150‑1 (MQ=255)
                                                                                                            gAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgc                                        <  2:453015/150‑1 (MQ=255)
                                                                                                                   tCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCa                                 >  2:1497640/1‑150 (MQ=255)
                                                                                                                   tCAGGATCTAGGCGATAGAATCGCGGGCAGCAGGCGATTACTCACGAAGATACGACTCGTGCATGGCCCGGCTGAAATCGACGGGGTGGTAGTCCGAGATGAATTGCAGGGCGCCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCa                                 >  2:931699/1‑150 (MQ=255)
                                                                                                                     aGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACg                               <  1:377667/150‑1 (MQ=255)
                                                                                                                       gATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAGATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGt                             <  2:1157235/150‑1 (MQ=255)
                                                                                                                           tGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATAACGGTAATg                         >  1:145535/1‑150 (MQ=255)
                                                                                                                                gATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGTGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcg                    >  2:1038540/1‑150 (MQ=255)
                                                                                                                                 gaggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcgc                   >  2:200010/3‑150 (MQ=255)
                                                                                                                                 aTGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcgc                   >  2:1535169/1‑150 (MQ=255)
                                                                                                                                   ggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                                  <  2:1713058/117‑1 (MQ=255)
                                                                                                                                   ggAATCGCCTGCTGCACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                                  >  1:1713058/1‑117 (MQ=255)
                                                                                                                                      aTCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGctcc              <  2:1306619/150‑1 (MQ=255)
                                                                                                                                           cgTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCtt          >  1:1050858/1‑150 (MQ=255)
                                                                                                                                                  cagGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACg   <  1:839553/150‑1 (MQ=255)
                                                                                                                                                   agGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCtt          <  1:171303/142‑1 (MQ=255)
                                                                                                                                                   agGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCtt          >  2:171303/1‑142 (MQ=255)
                                                                                                                                                   agGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACgg  <  2:1246639/150‑1 (MQ=255)
                                                                                                                                                   agGCGATTCCGCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGTGTAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACgg  <  1:595546/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGG  >  NC_002947/1039126‑1039421

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: