Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,545,698 C→A 8.6% I25I (ATC→ATA nfuA → iron‑sulfur cluster carrier protein NfuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,545,6980CA 8.6% 341.5 / 13.7 103I25I (ATC→ATAnfuAiron‑sulfur cluster carrier protein NfuA
Reads supporting (aligned to +/- strand):  ref base C (46/48);  new base A (3/6);  total (49/54)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

GTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATACCTG  >  NC_000913/3545552‑3545842
                                                                                                                                                  |                                                                                                                                                
gTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCcgc                                                                                                                                               <  2:1182519/150‑1 (MQ=255)
 tCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCcgcg                                                                                                                                              <  1:1374573/150‑1 (MQ=255)
  cGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGt                                                                                                                                             >  2:1395358/1‑150 (MQ=255)
  cGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGAACAAGCGCACTTTGCAAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCcgcga                                                                                                                                             >  1:1415743/1‑149 (MQ=255)
   gAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTa                                                                                                                                            <  1:1374970/150‑1 (MQ=255)
       cTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTtg                                                                                                                                        <  1:647023/150‑1 (MQ=255)
        ttGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTtgt                                                                                                                                       <  2:16326/150‑1 (MQ=255)
        ttGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTtgt                                                                                                                                       <  1:296618/150‑1 (MQ=255)
         tGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTtgtg                                                                                                                                      <  1:1395358/150‑1 (MQ=255)
         tGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTtgtg                                                                                                                                      <  1:262177/150‑1 (MQ=255)
         tGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTtgtg                                                                                                                                      >  1:1432988/1‑150 (MQ=255)
            agagCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATACGCGGATTtgtgtcc                                                                                                                                   <  2:645012/150‑4 (MQ=255)
            agagCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATACGCGGATTtgtgtcc                                                                                                                                   <  2:1407832/150‑4 (MQ=255)
            agagCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATACGCGGATTtgtgtcc                                                                                                                                   <  1:748835/150‑4 (MQ=255)
            agagCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATACGCGGATTtgtgtcc                                                                                                                                   <  2:423188/150‑4 (MQ=255)
            agagCCTCGATGATGGGCGTATTATAACCAACTAAAATAGGCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGAtt                                                                                                                                   <  1:489374/150‑1 (MQ=255)
             gagCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATACGCGGATTtgtgtccc                                                                                                                                  <  1:312563/150‑5 (MQ=255)
               gCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAAc                                                                                                                                >  1:769017/1‑150 (MQ=255)
               gCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAAc                                                                                                                                >  1:524279/1‑150 (MQ=255)
                   cGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTg                                                                                                                            <  1:162927/150‑1 (MQ=255)
                          gggCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGcc                                                                                                                     <  2:1063/150‑1 (MQ=255)
                           ggCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCt                                                                                                                    >  2:284743/1‑150 (MQ=255)
                             cGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTaa                                                                                                                  >  1:1465006/1‑150 (MQ=255)
                                attatAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGc                                                                                                               <  1:696104/150‑1 (MQ=255)
                                attatAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGc                                                                                                               >  2:67166/1‑150 (MQ=255)
                                   ataACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGa                                                                                                            >  2:697761/1‑150 (MQ=255)
                                     aaCCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAAt                                                                                                          >  2:517686/1‑150 (MQ=255)
                                      aCCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAAtg                                                                                                         >  1:653428/1‑150 (MQ=255)
                                            tAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGt                                                                                                   >  1:526349/1‑150 (MQ=255)
                                               aaTAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTc                                                                                                >  2:763200/1‑150 (MQ=255)
                                                 tAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCtt                                                                                              >  2:888021/1‑150 (MQ=255)
                                                  aGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTa                                                                                             >  1:1190316/1‑150 (MQ=255)
                                                   gTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTAt                                                                                            >  1:1428967/1‑150 (MQ=255)
                                                   gTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTAt                                                                                            <  1:636305/150‑1 (MQ=255)
                                                   gTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTAt                                                                                            >  1:1279758/1‑150 (MQ=255)
                                                     cAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCGAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTg                                                                                          >  1:538395/1‑150 (MQ=255)
                                                     cAAATATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGATGAATGTGGCGTTTCTTAttt                                                                                          >  1:501502/1‑149 (MQ=255)
                                                      aaCTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGt                                                                                         <  1:284743/150‑1 (MQ=255)
                                                      aaCTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGt                                                                                         <  2:32701/150‑1 (MQ=255)
                                                        cTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCt                                                                                               >  2:742805/1‑142 (MQ=255)
                                                        cTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCt                                                                                               <  1:742805/142‑1 (MQ=255)
                                                        cTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTcc                                                                                       >  1:1378834/1‑150 (MQ=255)
                                                           ttAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTccgcc                                                                                    <  2:1428967/150‑1 (MQ=255)
                                                            tAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTccgccg                                                                                   <  1:697761/150‑1 (MQ=255)
                                                               gCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGAc                                                                                >  1:1227786/1‑150 (MQ=255)
                                                                 cATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGc                                                                              >  1:850200/1‑150 (MQ=255)
                                                                     aCTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCtgtg                                                                          <  2:1271004/150‑1 (MQ=255)
                                                                      cTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTgg                                                                         <  1:229594/150‑1 (MQ=255)
                                                                      cTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTgg                                                                         <  1:510083/150‑1 (MQ=255)
                                                                      cTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTgg                                                                         <  2:312463/150‑1 (MQ=255)
                                                                           aTCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGcc                                                                    >  1:1453465/1‑150 (MQ=255)
                                                                             ccGTATTTCCTATGCTGCACTAGCGAACTTTGCCAATCTGCTGGGAAATCAGTAAGAATGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATTTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGacac                                                                  <  2:1415743/150‑1 (MQ=255)
                                                                              cGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCAcc                                                                 <  1:692379/150‑1 (MQ=255)
                                                                               gTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCATGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCg                                                                <  1:747381/150‑1 (MQ=255)
                                                                               gTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCg                                                                <  1:386911/150‑1 (MQ=255)
                                                                                  tttCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGATGTGGAAGCCACCGaca                                                             <  2:93658/150‑1 (MQ=255)
                                                                                       gATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGcc                                                        >  1:671840/1‑150 (MQ=255)
                                                                                       gATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGccc                                                       <  1:777534/150‑1 (MQ=255)
                                                                                           ctgcACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGa                                                    >  1:744811/1‑150 (MQ=255)
                                                                                             gcACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTTGAAGCCACCGACACAGCCCTGaaa                                                  <  1:1320582/150‑1 (MQ=255)
                                                                                             gcACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGaaa                                                  <  2:536641/150‑1 (MQ=255)
                                                                                               acaAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAAtt                                                >  1:486908/1‑150 (MQ=255)
                                                                                                 aaGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTg                                              <  1:1043353/150‑1 (MQ=255)
                                                                                                 aaGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTg                                              <  1:1235673/150‑1 (MQ=255)
                                                                                                  aGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGa                                             >  2:925501/1‑150 (MQ=255)
                                                                                                  aGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGa                                             <  2:1432988/150‑1 (MQ=255)
                                                                                                   gcgcACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGAc                                            <  1:888021/150‑1 (MQ=255)
                                                                                                   gcgcACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGAc                                            >  1:322363/1‑150 (MQ=255)
                                                                                                    cgcACTTTGCCAAACTTCTGGCAAATCAGGAAGAAGGGACACAAATACGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGTCGTTTCTTATTGTCCGCCGGACGCTGTGTAAGCCACCGACACAGCCCTGAAATTTGAcc                                           <  1:545536/150‑1 (MQ=255)
                                                                                                    cgcACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGAcc                                           <  1:979866/150‑1 (MQ=255)
                                                                                                    cgcACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGAcc                                           <  1:1171704/150‑1 (MQ=255)
                                                                                                       aCTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACAGAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACctgc                                        <  2:534439/150‑1 (MQ=255)
                                                                                                            gCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGAcc                                   >  1:253323/1‑150 (MQ=255)
                                                                                                             ccAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACcg                                  >  2:1046825/1‑150 (MQ=255)
                                                                                                             ccAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACcg                                  >  1:141114/1‑150 (MQ=255)
                                                                                                              cAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACcgc                                 >  1:1349487/1‑150 (MQ=255)
                                                                                                              cAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCATGAAATTTGACCTGCTGACcgc                                 >  1:1329226/1‑150 (MQ=255)
                                                                                                               aaaCTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACcgcg                                >  2:626999/1‑150 (MQ=255)
                                                                                                                     ctgGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGtt                          <  1:925501/150‑1 (MQ=255)
                                                                                                                       ggCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTtga                        >  1:733605/1‑150 (MQ=255)
                                                                                                                       ggCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTtga                        <  2:526349/150‑1 (MQ=255)
                                                                                                                        gCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTtgat                       <  2:524279/150‑1 (MQ=255)
                                                                                                                         cAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTtgatg                      <  2:653428/150‑1 (MQ=255)
                                                                                                                         cAAATAAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTtgatg                      >  1:415424/1‑150 (MQ=255)
                                                                                                                           aaTCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAg                    >  2:342573/1‑150 (MQ=255)
                                                                                                                           aaTCAGGAAGAAGGGAAACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAg                    >  1:232786/1‑150 (MQ=255)
                                                                                                                            aTCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACctgctg                                      >  1:352401/1‑131 (MQ=255)
                                                                                                                            aTCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACctgctg                                      <  2:352401/131‑1 (MQ=255)
                                                                                                                              cAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTa                 <  2:1378834/150‑1 (MQ=255)
                                                                                                                              cAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTa                 >  2:510052/1‑150 (MQ=255)
                                                                                                                                ggAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAg               <  2:1199043/150‑1 (MQ=255)
                                                                                                                                   agaagGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAgcgc            >  2:72091/1‑150 (MQ=255)
                                                                                                                                    gaagGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTAAGTTGATGAGTTAAGCGCa           <  2:573232/150‑1 (MQ=255)
                                                                                                                                     aagGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCAc          <  2:1071853/150‑1 (MQ=255)
                                                                                                                                     aagGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCAc          <  1:1046825/150‑1 (MQ=255)
                                                                                                                                      ttaatcacaAATACGCGGATTTGTGATTATCCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCAcc         >  1:135118/6‑150 (MQ=255)
                                                                                                                                      agGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCAcc         >  2:49410/1‑150 (MQ=255)
                                                                                                                                         atcacaAATACGCGGATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATa      >  2:1227613/3‑150 (MQ=255)
                                                                                                                                          tcacaAATACGCGGATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATAc     >  2:406063/2‑150 (MQ=255)
                                                                                                                                             caAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATACCTg  >  1:1210853/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATACCTG  >  NC_000913/3545552‑3545842

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