Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,458,110 C→A 100% T578T (ACC→ACA gspD → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,1100CA97.8% 2137.9 / 22.2 566T578T (ACC→ACAgspDgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base A (269/283);  total (275/291)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 574 total aligned reads displayed.

TACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457963‑3458256
                                                                                                                                                   |                                                                                                                                                  
tACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGa                                                                                                                                                  <  2:5408292/150‑1 (MQ=255)
tACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGa                                                                                                                                                  <  1:6248137/150‑1 (MQ=255)
tACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGa                                                                                                                                                  <  2:349643/150‑1 (MQ=255)
tACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGa                                                                                                                                                  <  2:1834818/150‑1 (MQ=255)
  cgagtCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGGTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgc                                                                                                                                                <  1:1076845/146‑1 (MQ=255)
     ttCAAAACGCCGTGCTGGTCAACACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcg                                                                                                                                             <  1:4323793/150‑1 (MQ=255)
     ttCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcg                                                                                                                                             >  1:4083433/1‑150 (MQ=255)
     ttCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcg                                                                                                                                             >  1:7713802/1‑150 (MQ=255)
      tCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcgc                                                                                                                                            >  1:2576842/1‑150 (MQ=255)
      tCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcgc                                                                                                                                            <  1:321876/150‑1 (MQ=255)
       cAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCt                                                                                                                                           >  2:6383332/1‑150 (MQ=255)
        aaaaCGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTa                                                                                                                                          <  2:4964707/150‑1 (MQ=255)
           aCGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAAc                                                                                                                                       <  1:4996940/150‑1 (MQ=255)
           aCGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAAc                                                                                                                                       <  1:71273/150‑1 (MQ=255)
             gCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGc                                                                                                                                     <  1:1498048/150‑1 (MQ=255)
             gCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGc                                                                                                                                     <  1:8015833/150‑1 (MQ=255)
               cGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCaa                                                                                                                                   <  2:3717648/150‑1 (MQ=255)
               cGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCaa                                                                                                                                   >  1:9424968/1‑150 (MQ=255)
                 tGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAAcc                                                                                                                                 >  2:9108843/1‑150 (MQ=255)
                    tGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGa                                                                                                                              >  1:3169194/1‑150 (MQ=255)
                    tGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGa                                                                                                                              >  2:7663554/1‑150 (MQ=255)
                    tGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGa                                                                                                                              >  1:8895202/1‑150 (MQ=255)
                    tGGTCAAAACCGGTGAAACGGTGGTCATGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCATCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGa                                                                                                                              >  1:7439856/1‑150 (MQ=255)
                     ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGAt                                                                                                                             <  1:4009781/150‑1 (MQ=255)
                     ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGAt                                                                                                                             <  2:5538447/150‑1 (MQ=255)
                     ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGAt                                                                                                                             <  1:6502672/150‑1 (MQ=255)
                     ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGAt                                                                                                                             <  1:160598/150‑1 (MQ=255)
                     ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGAt                                                                                                                             <  1:9037656/150‑1 (MQ=255)
                      gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            >  2:7367962/1‑150 (MQ=255)
                      gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            <  1:1308677/150‑1 (MQ=255)
                      gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGAGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            >  1:5558009/1‑150 (MQ=255)
                       tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgc                                                                                                                                              >  1:2990772/1‑131 (MQ=255)
                       tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATgg                                                                                                                           >  2:5669749/1‑150 (MQ=255)
                       tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATgg                                                                                                                           >  1:5078910/1‑150 (MQ=255)
                        cAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGt                                                                                                                          <  2:1257988/150‑1 (MQ=255)
                              cGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAt                                                                                                                    <  1:6115549/150‑1 (MQ=255)
                                 tGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccg                                                                                                                 >  1:83007/1‑150 (MQ=255)
                                      cGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCcga                                                                                                            <  2:2604646/150‑1 (MQ=255)
                                              cTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATc                                                                                                    <  2:7006264/150‑1 (MQ=255)
                                                gggCGGATTGCTGGATGATTTTTCTAAAGAGCACGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCg                                                                                                  >  2:8740842/1‑150 (MQ=255)
                                                gggCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCg                                                                                                  >  2:1446736/1‑150 (MQ=255)
                                                  gCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTg                                                                                                >  1:968566/1‑150 (MQ=255)
                                                      aTTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACga                                                                                            >  1:7033255/1‑150 (MQ=255)
                                                          cTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACgatgat                                                                                        >  1:5355835/1‑150 (MQ=255)
                                                          cTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCATTTAGTGGGGAAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACgatgat                                                                                        >  1:7802290/1‑150 (MQ=255)
                                                           tGGATTATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGatgatg                                                                                       <  2:6213093/150‑1 (MQ=255)
                                                                 atTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATc                                                                                 <  1:3909117/150‑1 (MQ=255)
                                                                  tttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCg                                                                                >  1:8155547/1‑150 (MQ=255)
                                                                     ttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTc                                                                             <  2:2066765/150‑1 (MQ=255)
                                                                     ttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACACAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTAATCGCTc                                                                             >  1:8960673/1‑150 (MQ=255)
                                                                      tCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCa                                                                            <  2:2685965/150‑1 (MQ=255)
                                                                      tCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCa                                                                            <  2:8237570/150‑1 (MQ=255)
                                                                        tAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACt                                                                          >  1:5168575/1‑150 (MQ=255)
                                                                         aaaGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTg                                                                         >  1:2943998/1‑150 (MQ=255)
                                                                         aaaGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTg                                                                         >  1:1832160/1‑150 (MQ=255)
                                                                           agagCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTc                                                                       <  1:1172478/150‑1 (MQ=255)
                                                                           agagCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTc                                                                       <  1:3633490/150‑1 (MQ=255)
                                                                            gagCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCa                                                                      <  1:4436431/150‑1 (MQ=255)
                                                                               cAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCaaaa                                                                   >  1:7568495/1‑150 (MQ=255)
                                                                                 aGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAaga                                                                 <  2:3112156/150‑1 (MQ=255)
                                                                                 aGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAaga                                                                 <  2:808378/150‑1 (MQ=255)
                                                                                     tCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaaa                                                             <  1:2279456/150‑1 (MQ=255)
                                                                                       aaaGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATa                                                           >  1:8756175/1‑150 (MQ=255)
                                                                                       aaaGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATa                                                           >  1:246380/1‑150 (MQ=255)
                                                                                        aaGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATac                                                          >  1:5955005/1‑150 (MQ=255)
                                                                                         aGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATaca                                                         <  1:1966850/150‑1 (MQ=255)
                                                                                                tctGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTa                                                  >  2:9569997/1‑150 (MQ=255)
                                                                                                tctGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTa                                                  >  1:8531861/1‑150 (MQ=255)
                                                                                                      tGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCa                                            <  2:8531861/150‑1 (MQ=255)
                                                                                                       ggCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAg                                           <  2:7951655/150‑1 (MQ=255)
                                                                                                         cGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGa                                         <  2:4524337/150‑1 (MQ=255)
                                                                                                              ttCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaac                                    <  1:6771836/150‑1 (MQ=255)
                                                                                                                 cTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTCAACGCACCCTGATGGTCTTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGACCTACACCCGTTACCGTCAGGCGcaccacc                                 >  2:5234920/1‑148 (MQ=255)
                                                                                                                 cTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCaacaac                                 >  2:3293454/1‑150 (MQ=255)
                                                                                                                 cTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCaacaac                                 >  2:5919185/1‑150 (MQ=255)
                                                                                                                  tttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAacaaca                                >  1:6381474/1‑150 (MQ=255)
                                                                                                                  tttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAacaaca                                >  2:5048148/1‑150 (MQ=255)
                                                                                                                  tttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAacaaca                                >  2:426149/1‑150 (MQ=255)
                                                                                                                  tttAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAacaaca                                >  1:3469931/1‑150 (MQ=255)
                                                                                                                     aGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCg                             <  2:4697871/150‑1 (MQ=255)
                                                                                                                       tGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGa                           <  2:1667865/150‑1 (MQ=255)
                                                                                                                            cAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGAc                      >  1:7235857/1‑150 (MQ=255)
                                                                                                                            cAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGAc                      >  1:957999/1‑150 (MQ=255)
                                                                                                                             aaCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACg                     <  1:4611588/150‑1 (MQ=255)
                                                                                                                             aaCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACg                     >  2:1776802/1‑150 (MQ=255)
                                                                                                                              aCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACgg                    >  1:6410552/1‑150 (MQ=255)
                                                                                                                              aCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACgg                    <  2:9278798/150‑1 (MQ=255)
                                                                                                                                tctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGa                  >  1:2783729/1‑150 (MQ=255)
                                                                                                                                 ctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGaa                 <  2:7077221/150‑1 (MQ=255)
                                                                                                                                   tCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAAt               >  2:8943134/1‑150 (MQ=255)
                                                                                                                                      gCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCaa            >  2:4101051/1‑150 (MQ=255)
                                                                                                                                            ccTCCCCAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAgcgc      <  2:6603146/150‑1 (MQ=255)
                                                                                                                                            ccTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAgcgc      >  2:2454965/1‑150 (MQ=255)
                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGACATCAAAAGCGCt     >  1:921917/1‑150 (MQ=255)
                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     >  2:2937070/1‑150 (MQ=255)
                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     >  2:75387/1‑150 (MQ=255)
                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     >  1:5708934/1‑150 (MQ=255)
                                                                                                                                              tCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTg    >  1:5817259/1‑150 (MQ=255)
                                                                                                                                                cacaGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGt  <  2:3398214/150‑1 (MQ=255)
                                                                                                                                                cacaGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGt  <  1:5416904/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457963‑3458256

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