Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,719,306:1 +T 5.7% coding (390/585 nt) nfuA → Fe/S biogenesis protein NfuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,719,3061.T 5.7% 352.7 / 14.2 88coding (390/585 nt)nfuAFe/S biogenesis protein NfuA
Reads supporting (aligned to +/- strand):  ref base . (39/44);  new base T (3/2);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGG  >  NC_002947/2719160‑2719451
                                                                                                                                                   |                                                                                                                                                 
gACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTATCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATggcggc                                                                                                                                                <  2:1268483/150‑1 (MQ=255)
 aCTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGcc                                                                                                                                               <  1:1286634/150‑1 (MQ=255)
     cGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGt                                                                                                                                           >  1:1043222/1‑150 (MQ=255)
       ccACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGa                                                                                                                                         <  2:264443/150‑1 (MQ=255)
        cACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAg                                                                                                                                        >  2:742543/1‑150 (MQ=255)
        cACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAg                                                                                                                                        >  2:201768/1‑150 (MQ=255)
            gaccgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtg                                                                                                                                    <  1:1016305/150‑1 (MQ=255)
             accgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAg                                                                                                                                             <  1:93516/140‑1 (MQ=255)
             accgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAg                                                                                                                                             >  2:93516/1‑140 (MQ=255)
             accgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtgg                                                                                                                                   <  2:108857/150‑1 (MQ=255)
             accgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtgg                                                                                                                                   >  2:809443/1‑150 (MQ=255)
               cgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtggtg                                                                                                                                 <  2:848156/150‑1 (MQ=255)
                  aTGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGAc                                                                                                                              >  1:621512/1‑150 (MQ=255)
                   tGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACg                                                                                                                             >  1:1333806/1‑150 (MQ=255)
                       tGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGt                                                                                                                         >  1:206148/1‑150 (MQ=255)
                       tGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGt                                                                                                                         >  2:272065/1‑150 (MQ=255)
                           cAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGAt                                                                                                                     <  1:1266313/150‑1 (MQ=255)
                            aGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATg                                                                                                                    >  2:779556/1‑150 (MQ=255)
                              cTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtggtg                                                                                                                                 <  2:728777/135‑1 (MQ=255)
                              cTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtggtg                                                                                                                                 >  1:728777/1‑135 (MQ=255)
                               tGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACg                                                                                                                 <  1:116748/150‑1 (MQ=255)
                                    aTCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGTCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCAtt                                                                                                            >  2:782537/1‑150 (MQ=255)
                                               gatctccAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCa                                                                                                 <  1:661227/145‑1 (MQ=255)
                                                 aCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGt                                                                                               <  1:809443/150‑1 (MQ=255)
                                                    ccAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCg                                                                                            <  2:1252935/150‑1 (MQ=255)
                                                     cAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCgg                                                                                           >  2:972843/1‑150 (MQ=255)
                                                     cAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCgg                                                                                           >  2:495267/1‑150 (MQ=255)
                                                      aaGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGt                                                                                                                            >  1:1014044/1‑116 (MQ=255)
                                                      aaGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGt                                                                                                                            <  2:1014044/116‑1 (MQ=255)
                                                      aaGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggt                                                                                          <  1:830507/150‑1 (MQ=255)
                                                       aGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtg                                                                                         <  1:495267/150‑1 (MQ=255)
                                                       aGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtg                                                                                         <  2:1004274/150‑1 (MQ=255)
                                                          tggcGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGgtggtg                                                                                      >  1:1421214/4‑150 (MQ=255)
                                                             cGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCt                                                                                   >  2:377241/1‑150 (MQ=255)
                                                               aTGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGc                                                                                 <  2:222268/150‑1 (MQ=255)
                                                                   tCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAgg                                                                             <  2:1327011/150‑1 (MQ=255)
                                                                        gAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGAt                                                                                                                     >  1:731011/1‑105 (MQ=255)
                                                                        gAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGAt                                                                                                                     <  2:731011/105‑1 (MQ=255)
                                                                        gAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGc                                                                        >  2:1353887/1‑150 (MQ=255)
                                                                         aaGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGCCCCAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCg                                                                       <  1:305356/150‑1 (MQ=255)
                                                                             cAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCa                                                                   <  2:1421214/150‑1 (MQ=255)
                                                                                cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGc                                                                >  1:2189/1‑150 (MQ=255)
                                                                                cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGc                                                                >  2:16913/1‑150 (MQ=255)
                                                                                cccGATCAACGAACGTATCAACTACTACCTAAAGACCGAGATCAACCTCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCCGGTGGACGTGCTCGATGACGGCATTGCCGTACTGCAGTTCTGTGGTCGCTGCCAGGGTTTCGGCCAGGc                                                                >  2:582464/1‑150 (MQ=255)
                                                                                   gATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGa                                                             >  1:328801/1‑150 (MQ=255)
                                                                                       aacgaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTc                                                         >  2:624388/1‑150 (MQ=255)
                                                                                         cgaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAc                                                       <  1:1019085/150‑1 (MQ=255)
                                                                                          gaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCAcc                                                      <  1:377241/150‑1 (MQ=255)
                                                                                              gTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTga                                                  >  1:978733/1‑150 (MQ=255)
                                                                                               tATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgat                                                 <  2:978733/150‑2 (MQ=255)
                                                                                               tATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                                 <  1:742543/150‑1 (MQ=255)
                                                                                               tATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                                 <  1:779556/150‑1 (MQ=255)
                                                                                                    actactACCTACAGACAGAGGTCAACCCCGGCCTGGACAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGAAGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaag                                            <  1:782537/150‑1 (MQ=255)
                                                                                                     ctactaCCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGaaggaagg                                           <  2:1052460/150‑1 (MQ=255)
                                                                                                       actacCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCa                                         <  1:201768/150‑1 (MQ=255)
                                                                                                            ctacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAg                                    >  1:779362/1‑150 (MQ=255)
                                                                                                             tacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCGAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                                   <  2:972102/150‑1 (MQ=255)
                                                                                                                 gACCGTGAGCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAc                               >  2:1064766/1‑150 (MQ=255)
                                                                                                                  aCCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtggt                                                                                       <  1:1009466/93‑1 (MQ=255)
                                                                                                                  aCCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtggt                                                                                       >  2:1009466/1‑93 (MQ=255)
                                                                                                                    cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCa                                                        <  2:946564/122‑1 (MQ=255)
                                                                                                                    cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCa                                                        >  1:946564/1‑122 (MQ=255)
                                                                                                                     gagaTCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTg                           >  2:1314365/1‑150 (MQ=255)
                                                                                                                        aTCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCGAGGGTTGCGCGCAGGCCGATGTCACCCTGAAGGACGTCATCGAGCGCACCCTGATc                        >  2:347203/1‑150 (MQ=255)
                                                                                                                        aTCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGAACCCTGCTc                        >  2:1170080/1‑150 (MQ=255)
                                                                                                                          cAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGa                      <  2:1043222/150‑1 (MQ=255)
                                                                                                                           aaCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAg                     <  2:2189/150‑1 (MQ=255)
                                                                                                                           aaCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAg                     >  2:678830/1‑150 (MQ=255)
                                                                                                                            aCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGa                                     <  2:267655/133‑1 (MQ=255)
                                                                                                                            aCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGa                                     >  1:267655/1‑133 (MQ=255)
                                                                                                                            aCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGc                    >  2:1312504/1‑150 (MQ=255)
                                                                                                                            aCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGc                    <  2:1333806/150‑1 (MQ=255)
                                                                                                                            aCCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGc                    <  2:779362/150‑1 (MQ=255)
                                                                                                                              cccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGt                  >  2:769175/1‑150 (MQ=255)
                                                                                                                                tcaccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                                            >  2:1240363/4‑106 (MQ=255)
                                                                                                                                tcaccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                                            <  1:1240363/103‑1 (MQ=255)
                                                                                                                                tcaccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgacggaa                                             >  2:83936/4‑121 (MQ=255)
                                                                                                                                tcaccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaa                                             <  1:83936/118‑1 (MQ=255)
                                                                                                                                tcaccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTAt                >  1:977864/4‑150 (MQ=255)
                                                                                                                                    cTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGa                                                             >  2:812225/1‑101 (MQ=255)
                                                                                                                                    cTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGa                                                             <  1:812225/101‑1 (MQ=255)
                                                                                                                                      ggCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa          <  1:972843/150‑1 (MQ=255)
                                                                                                                                      ggCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa          <  2:621512/150‑1 (MQ=255)
                                                                                                                                       gccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAg         <  1:515880/150‑1 (MQ=255)
                                                                                                                                       gccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAg         <  2:328801/150‑1 (MQ=255)
                                                                                                                                            ccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGaa    >  1:488019/1‑150 (MQ=255)
                                                                                                                                            ccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGaa    <  1:272065/150‑1 (MQ=255)
                                                                                                                                             caTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAg   <  1:16913/150‑1 (MQ=255)
                                                                                                                                              aTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAgg  >  1:928236/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGG  >  NC_002947/2719160‑2719451

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: