Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,405,885 G→C 5.4% E150D (GAG→GAC accB → biotin carboxyl carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,405,8850GC 5.4% 351.6 / 7.1 106E150D (GAG→GACaccBbiotin carboxyl carrier protein
Reads supporting (aligned to +/- strand):  ref base G (48/51);  new base C (3/3);  total (51/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 106 total aligned reads displayed.

GTTCATCGAAGTGGGTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATC  >  NC_000913/3405738‑3406025
                                                                                                                                                   |                                                                                                                                            
gTTCATCGAAGTGGGTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGcc                                                                                                                                            >  2:1536266/1‑150 (MQ=255)
gTTCATCGAAGTGGGTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGcc                                                                                                                                            >  1:1370861/1‑150 (MQ=255)
     tCGAAGTGGGTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTgg                                                                                                                                       >  1:998596/1‑150 (MQ=255)
      cGAAGTGGGTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgt                                                                                                                                      <  2:662368/150‑1 (MQ=255)
       gAAGTGGGTCAGAAAGTCAACGTGGGCGATACCCTGTGCGTCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgtc                                                                                                                                     >  1:1551112/1‑150 (MQ=255)
           gggggCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGCACAAATCCGCTACCGTGAAAGCAATTCTGGCCGAAAGTGGACAACCGGTAGAATTTGACGACCAGCGGCTCGTCa                                                                                                                                 <  2:520423/149‑1 (MQ=255)
            gggTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCAt                                                                                                                                >  1:301216/1‑150 (MQ=255)
            gggTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCAt                                                                                                                                <  2:301216/150‑1 (MQ=255)
            gggTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGACCAGCGGCTCGTCaa                                                                                                                                <  2:1036452/150‑2 (MQ=255)
             ggTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATc                                                                                                                               >  1:784815/1‑150 (MQ=255)
              gTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCg                                                                                                                              <  2:483757/150‑1 (MQ=255)
              gTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCg                                                                                                                              <  2:454600/150‑1 (MQ=255)
                cAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAg                                                                                                                            >  1:1159741/1‑150 (MQ=255)
                cAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAg                                                                                                                            <  2:111856/150‑1 (MQ=255)
                cAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAg                                                                                                                            >  1:1268337/1‑150 (MQ=255)
                     aGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACAGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACg                                                                                                                       >  2:629118/1‑150 (MQ=255)
                      gTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGa                                                                                                                      >  2:1038756/1‑150 (MQ=255)
                       tCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAg                                                                                                                     <  2:517785/150‑1 (MQ=255)
                         aaCGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGc                                                                                                                   <  2:1129494/150‑1 (MQ=255)
                              gggCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACa                                                                                                              >  2:1383928/1‑150 (MQ=255)
                               ggagATACCCTGTTGATCGTTGAAGCCATGAATATGATGAAACAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTTACGAGCCTATGGTCGTCATCGAGTAACGAGGCGAAGAt                                                                                                             <  2:497578/147‑1 (MQ=255)
                                    tACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTgg                                                                                                        <  1:1439738/150‑1 (MQ=255)
                                     aCCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGa                                                                                                       >  1:38086/1‑150 (MQ=255)
                                      cccTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGAt                                                                                                      <  1:253860/150‑1 (MQ=255)
                                          gtgCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATaaaa                                                                                                  >  2:449664/1‑150 (MQ=255)
                                           tgCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAAt                                                                                                 <  1:899224/150‑1 (MQ=255)
                                            gCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAAtt                                                                                                >  1:1557389/1‑150 (MQ=255)
                                                cGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTa                                                                                            <  2:1068326/150‑1 (MQ=255)
                                                cGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTa                                                                                            <  2:1370861/150‑1 (MQ=255)
                                                cGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTa                                                                                            >  2:945069/1‑150 (MQ=255)
                                                cGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTa                                                                                            >  1:303973/1‑150 (MQ=255)
                                                   tGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTg                                                                                         <  1:449664/150‑1 (MQ=255)
                                                   tGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTg                                                                                         <  1:1614230/150‑1 (MQ=255)
                                                    gAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGc                                                                                        <  2:521090/150‑1 (MQ=255)
                                                    gAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGc                                                                                        >  1:926260/1‑150 (MQ=255)
                                                     aaGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGcc                                                                                       <  2:1159741/150‑1 (MQ=255)
                                                        ccATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAAc                                                                                    >  2:925077/1‑150 (MQ=255)
                                                         cATGAAAATGATGAACCAGATCTAAGCGGACAAATCCGGTACCGTGAAAGAAATTCTGGTAGAAAGTGGACAACATGTAAAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAAcc                                                                                   >  2:1565478/1‑150 (MQ=255)
                                                         cATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGGGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAAcc                                                                                   >  2:506475/1‑150 (MQ=255)
                                                         cATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAAcc                                                                                   >  2:550619/1‑150 (MQ=255)
                                                         cATCAAAATCATGAACCAGATCGAACCGGACAAATCCGGTACCGTGAAACCAATTCTTGACGAAAGTGTACAACCCCTAGAATTTGACGACCCGCTGGTCGTCATCGACTAACGAGGCGAACATGATGGAAAAAATTCTTATTGCCAAcc                                                                                   >  1:728876/1‑150 (MQ=255)
                                                            gAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACcgcg                                                                                >  1:1260794/1‑150 (MQ=255)
                                                             aaaaTGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGATCATGCTGGATAAAATTGTTATTGCAAACCgcgg                                                                               >  2:841218/1‑150 (MQ=255)
                                                             aaaaTGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCgcgg                                                                               <  1:751381/150‑1 (MQ=255)
                                                             aaaaTGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCgcgg                                                                               <  1:1268657/150‑1 (MQ=255)
                                                               aaTGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCgcggcg                                                                             <  1:484865/150‑1 (MQ=255)
                                                                   atgaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGat                                                                         >  2:545778/1‑150 (MQ=255)
                                                                     gaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGattg                                                                       <  1:15889/150‑1 (MQ=255)
                                                                     gaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGattg                                                                       <  1:1513700/150‑1 (MQ=255)
                                                                       cccAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTGATTGCCACCCGCGGCGAGagtgca                                                                     >  2:1560742/2‑150 (MQ=255)
                                                                       aCCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGattgca                                                                     >  2:1200027/1‑150 (MQ=255)
                                                                       aCCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGattgca                                                                     >  2:1304533/1‑150 (MQ=255)
                                                                        ccAGATCGAAGCGGACAAATCCGGTACCGTGAAAGAAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGattgcat                                                                    >  2:1606164/1‑150 (MQ=255)
                                                                           gATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGattgcattgc                                                                 <  2:234299/150‑1 (MQ=255)
                                                                           gATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGattgcattgc                                                                 >  2:1024395/1‑150 (MQ=255)
                                                                              cGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTa                                                              >  2:1285669/1‑150 (MQ=255)
                                                                              cGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGAAGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGt                                                                                              >  2:231072/1‑118 (MQ=255)
                                                                              cGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGAAGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGt                                                                                              <  1:231072/118‑1 (MQ=255)
                                                                               gAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTAt                                                             >  2:392697/1‑150 (MQ=255)
                                                                               gAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTAt                                                             >  1:1029630/1‑150 (MQ=255)
                                                                               gAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTAt                                                             >  1:823051/1‑150 (MQ=255)
                                                                                   cGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTAttctt                                                         <  1:1523855/150‑1 (MQ=255)
                                                                                    ggACAAATCCGGTACCGTTAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCTCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTAttcttc                                                        >  2:429656/1‑150 (MQ=255)
                                                                                    ggACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTgg                                                                                                        >  1:73823/1‑102 (MQ=255)
                                                                                    ggACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTgg                                                                                                        <  2:73823/102‑1 (MQ=255)
                                                                                    ggACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTAttcttc                                                        >  2:1311263/1‑150 (MQ=255)
                                                                                    ggACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTAttcttc                                                        >  1:613500/1‑150 (MQ=255)
                                                                                     gACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCg                                                       <  2:188081/150‑1 (MQ=255)
                                                                                        aaaTCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGc                                                    <  2:132480/150‑1 (MQ=255)
                                                                                         aaTCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCTGCGAGATTGCATTGCGTATTCTTCGTGcc                                                   <  2:784815/150‑1 (MQ=255)
                                                                                         aaTCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGcc                                                   <  1:264735/150‑1 (MQ=255)
                                                                                           tCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGTGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTg                                                 <  2:342736/150‑1 (MQ=255)
                                                                                           tCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTg                                                 <  1:1200027/150‑1 (MQ=255)
                                                                                           tCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGGCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTg                                                 <  2:927691/150‑1 (MQ=255)
                                                                                             cGGTAGCGTGAAAGCAATTCAGGACGAAAGTGCACAACCGGTAGAATTAGACGAGCCGCTGGTAGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATAGTTATAGCCAACCGCGCCGAGAATGCAATGCGGATTCTTCGTGCCTGTa                                               >  2:828251/1‑150 (MQ=255)
                                                                                             cGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTa                                               <  2:303973/150‑1 (MQ=255)
                                                                                               gTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTaaa                                             <  1:392697/150‑1 (MQ=255)
                                                                                                   cGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAAc                                         >  1:157569/1‑150 (MQ=255)
                                                                                                       aaaGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTggg                                     <  2:913564/150‑1 (MQ=255)
                                                                                                            aaTTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCa                                >  2:507750/1‑150 (MQ=255)
                                                                                                            aaTTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCa                                <  1:1285669/150‑1 (MQ=255)
                                                                                                                cTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGAc                            <  2:926260/150‑1 (MQ=255)
                                                                                                                cTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGAc                            <  1:1383928/150‑1 (MQ=255)
                                                                                                                     cGATAGTGGACTACCGCTATAATTTTACGATCCTCTGGTCTTTTTCTAGTAACGAGTCGAAGATTCTGGAAAAAATTGTTATTGCCAACCCGGGCGTGATTGCATTGCGTATTCTTCGTGCCTGTACAGTACTGGGCATGAAGACTGTCg                       <  2:728876/150‑1 (MQ=255)
                                                                                                                     cGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCg                       >  2:831308/1‑150 (MQ=255)
                                                                                                                      gAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGc                      <  1:429656/150‑1 (MQ=255)
                                                                                                                           tGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGc                 <  1:545778/150‑1 (MQ=255)
                                                                                                                           tGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGc                 <  2:179236/150‑1 (MQ=255)
                                                                                                                           gcgACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGc                 <  1:323337/148‑1 (MQ=255)
                                                                                                                                  ccGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAg          >  2:1587253/1‑150 (MQ=255)
                                                                                                                                     gTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAgcgc       >  2:272839/1‑150 (MQ=255)
                                                                                                                                      tAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGcgcg      >  1:847384/1‑150 (MQ=255)
                                                                                                                                       aGAATTTGTCGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCgg     <  2:807842/150‑1 (MQ=255)
                                                                                                                                       aGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCgg     <  1:272839/150‑1 (MQ=255)
                                                                                                                                       aGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCgg     <  1:662992/150‑1 (MQ=255)
                                                                                                                                        gAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGa    <  1:1038756/150‑1 (MQ=255)
                                                                                                                                         tcgatGACGACCAGCGGCTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGCCTGTCGCTGTGCACTCCAGCGCGGAt   >  1:399419/5‑150 (MQ=255)
                                                                                                                                         tcgatGACGACCAGCGGCTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGAt   >  1:1273823/5‑150 (MQ=255)
                                                                                                                                         aaTTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGAt   >  1:513246/1‑150 (MQ=255)
                                                                                                                                          aTTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATc  <  2:1029630/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GTTCATCGAAGTGGGTCAGAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTATTGCCAACCGCGGCGAGATTGCATTGCGTATTCTTCGTGCCTGTAAAGAACTGGGCATCAAGACTGTCGCTGTGCACTCCAGCGCGGATC  >  NC_000913/3405738‑3406025

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