Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,458,110 C→A 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,1100CA99.3% 500.9 / NA 147T578T (ACC→ACA) gspDgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (101/45);  minor base T (1/0);  total (102/45)
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

Only 100 of 148 total aligned reads displayed.

ATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457967‑3458256
                                                                                                                                               |                                                                                                                                                  
aTTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgc                                                                                                                                              >  1:1736375/1‑150 (MQ=255)
aTTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgc                                                                                                                                              >  1:1605777/1‑150 (MQ=255)
 ttCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgcgcg                                                                                                                                             >  1:828463/1‑150 (MQ=255)
             tGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAAcc                                                                                                                                 >  1:169031/1‑150 (MQ=255)
              gCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCt                                                                                                                                >  1:629721/1‑150 (MQ=255)
               cTGGTCAAAACCGGTTAATGGGTGGCCCTGGGCGGTTTGGTGTATTAATTTTCTAAAGACCAATGGTCAAAGGTACCTCTGCTTGGCTTTATTCCTTTAGGGGGGCAACTCTTCCGCTATACCTCCACAGCGCGCGCAAACCGCAACCtc                                                                                                                               <  2:1709418/150‑2 (MQ=255)
               cTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTGCTGTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGCGCGCGCTAAACGCAACCTg                                                                                                                               >  2:1336345/1‑150 (MQ=255)
               cTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTg                                                                                                                               >  2:1674015/1‑150 (MQ=255)
                 ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGAGTGCTGGATGATGTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCGTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGAGAGCTAAACGCAACCTGAt                                                                                                                             >  1:1392269/1‑150 (MQ=255)
                  gTCAATACCGGTGAACCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGGCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            <  2:827907/150‑1 (MQ=255)
                   tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATgg                                                                                                                           <  1:33726/150‑1 (MQ=255)
                   tCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATgg                                                                                                                           >  2:1991502/1‑150 (MQ=255)
                     aaaaCCGGGGAAACGGTGGTCCTGGGCGGTTTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCATTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTa                                                                                                                         <  1:1124412/150‑1 (MQ=255)
                       aaCCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTATGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTAtt                                                                                                                       <  1:314262/150‑1 (MQ=255)
                           ggTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtc                                                                                                                   <  2:1141677/150‑1 (MQ=255)
                           cgggatccGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTGCTTTAGTGGGGCAACTCTACCGGTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtc                                                                                                                   <  2:179161/143‑1 (MQ=255)
                               aaaCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccgtc                                                                                                               >  1:550463/1‑150 (MQ=255)
                                 aCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccgtccg                                                                                                             <  2:1098054/150‑1 (MQ=255)
                                 aCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccgtccg                                                                                                             >  2:1847662/1‑150 (MQ=255)
                                    gtggtCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTAAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCcgacg                                                                                                          <  1:1229777/150‑1 (MQ=255)
                                        tCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGTCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGCCGATTa                                                                                                      >  2:989286/1‑150 (MQ=255)
                                        tCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCGAAACCCAACCTGATGGTATTTATCCGTCCGCCGCTTa                                                                                                      >  2:1247013/1‑150 (MQ=255)
                                          cTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCAGCTATACCTCCACAGAGCGCGCTAAACGCTACCTGATGGTATTTATCCGTCCGACGATTATc                                                                                                    >  2:1231577/1‑150 (MQ=255)
                                           tGGGCGTATTGCTGTATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATcc                                                                                                   <  2:851203/150‑1 (MQ=255)
                                           tGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATcc                                                                                                   >  1:510248/1‑150 (MQ=255)
                                           tGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCAGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATcc                                                                                                   >  1:357386/1‑150 (MQ=255)
                                            gggCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCg                                                                                                  <  1:1073696/150‑1 (MQ=255)
                                            gggCGGAGTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCg                                                                                                  <  1:1841952/150‑1 (MQ=255)
                                              gCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGCCGATTATCCGTg                                                                                                >  2:1926308/1‑150 (MQ=255)
                                                 gATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACg                                                                                             >  2:1874175/1‑150 (MQ=255)
                                                 gATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACg                                                                                             >  2:1149361/1‑150 (MQ=255)
                                                     gCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACgatga                                                                                         <  1:240195/150‑1 (MQ=255)
                                                      cTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACgatgat                                                                                        >  2:1264934/1‑150 (MQ=255)
                                                          atgatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCTCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGttt                                                                                    >  1:152619/1‑150 (MQ=255)
                                                          atgatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGttt                                                                                    >  2:1118207/1‑150 (MQ=255)
                                                           tgatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTa                                                                                   >  2:152631/1‑150 (MQ=255)
                                                             atTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGTCGATTATCCGTGACGATGATGTTTATc                                                                                 <  1:152631/150‑1 (MQ=255)
                                                              tttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCg                                                                                >  1:752515/1‑150 (MQ=255)
                                                              tttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGACCGCGCTAAACGCACCATGATGCTATTTAACAGACCGACGAGAAACCGTGACCATGATGTTTAGCg                                                                                >  2:1875597/1‑150 (MQ=255)
                                                               ttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGc                                                                               <  2:732345/150‑1 (MQ=255)
                                                                tttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCt                                                                              >  1:907328/1‑150 (MQ=255)
                                                                    tAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCATCCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACt                                                                          >  1:837007/1‑150 (MQ=255)
                                                                    tAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACt                                                                          >  2:1396225/1‑150 (MQ=255)
                                                                     aaaGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTg                                                                         >  2:1651429/1‑150 (MQ=255)
                                                                      aaGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAACCGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGTTGATGTTTATCGCTCACTGt                                                                        >  2:977121/1‑150 (MQ=255)
                                                                      aaGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGt                                                                        >  1:862771/1‑150 (MQ=255)
                                                                          gCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCaaa                                                                    >  1:1298483/1‑150 (MQ=255)
                                                                            aaGTGTCAAGGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAag                                                                  >  2:1482836/1‑150 (MQ=255)
                                                                            aaGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAag                                                                  <  1:917094/150‑1 (MQ=255)
                                                                             aGTGTCAAAGGTTCCTCTGCTTGTGGATATTCCTTTAGTGGGTCAACTCTTCCGCTATACCTCCAGAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAaga                                                                 <  2:67555/150‑1 (MQ=255)
                                                                             aGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAaga                                                                 >  2:1086584/1‑150 (MQ=255)
                                                                             aGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTATAGTGGGGAAACTCTACCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGCTATTTATCCGTCCGAAGATTATACGTGACGATGAAGTTTATCGCTCACTGTCAAAaga                                                                 >  1:1355212/1‑150 (MQ=255)
                                                                               tgtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAgaga                                                               >  2:861213/1‑150 (MQ=255)
                                                                                gtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaa                                                              >  2:904376/1‑150 (MQ=255)
                                                                                gtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaa                                                              >  1:1329660/1‑150 (MQ=255)
                                                                                gtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaa                                                              >  1:1384968/1‑150 (MQ=255)
                                                                                 tCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTCATGGTATTTATCAGTCCGACGATTATCCGTGACGATGATGTTTATCGCTAACTGTCAAAAGAGaaa                                                             >  1:1266067/1‑150 (MQ=255)
                                                                                  cAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAAt                                                            >  2:1430370/1‑150 (MQ=255)
                                                                                   aaaGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATa                                                           >  1:1746427/1‑150 (MQ=255)
                                                                                   aaaGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATAATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATa                                                           >  1:689480/1‑150 (MQ=255)
                                                                                               gCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACTCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCg                                               <  2:1220364/150‑1 (MQ=255)
                                                                                               gCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCg                                               <  1:750778/150‑1 (MQ=255)
                                                                                               gCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCg                                               >  2:133546/1‑150 (MQ=255)
                                                                                               gCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAAAACACCCGtatccg                                               >  2:335514/1‑145 (MQ=255)
                                                                                                 ttGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTc                                             >  2:1058758/1‑150 (MQ=255)
                                                                                                  tGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCCCAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTAACCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCa                                            >  2:198880/1‑150 (MQ=255)
                                                                                                  tGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCa                                            <  1:1149361/150‑1 (MQ=255)
                                                                                                  tGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCa                                            >  1:1962695/1‑150 (MQ=255)
                                                                                                    gCGATATTCCTTTAGTGGGGCATCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAgg                                          <  1:1449868/150‑1 (MQ=255)
                                                                                                    gCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAgg                                          <  2:1568246/150‑1 (MQ=255)
                                                                                                     cGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGa                                         >  2:1451171/1‑150 (MQ=255)
                                                                                                           tCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTGACTGTCAAAAGAGAAATAAACCCCTTACCGTCAGGAGcaaca                                   >  1:930768/1‑150 (MQ=255)
                                                                                                           tCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaaca                                   >  2:389844/1‑150 (MQ=255)
                                                                                                           tCCTTTAGGGGGGCCACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGGCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaaca                                   <  2:1735606/150‑1 (MQ=255)
                                                                                                            ccTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGcaacaa                                  >  1:682519/1‑150 (MQ=255)
                                                                                                            ccTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGAAAATCCGTGCCGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGGTACCGTCAGGAGcatcaa                                  >  2:1273989/1‑150 (MQ=255)
                                                                                                                tAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGc                              >  1:13430/1‑150 (MQ=255)
                                                                                                                 aGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATATATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCg                             >  1:1772935/1‑150 (MQ=255)
                                                                                                                   tGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTAAAAAGAGAAATACACCCGTTACCGTCAGTATCAACAACAGCGGa                           >  2:937094/1‑150 (MQ=255)
                                                                                                                    ggggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGTCGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGAt                          <  1:1480433/150‑1 (MQ=255)
                                                                                                                    ggggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAATTTCACCCGTTACCGTCAGGAGCCACAACAGCGGAt                          >  2:1883199/1‑150 (MQ=255)
                                                                                                                     gggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATc                         <  1:1939279/150‑1 (MQ=255)
                                                                                                                     gggCAACTCTTCCGCTATACCTCCACAGAGCGCACTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATc                         >  2:770228/1‑150 (MQ=255)
                                                                                                                      ggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCg                        >  2:1028136/1‑150 (MQ=255)
                                                                                                                       gCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGca                                      <  1:1088989/135‑1 (MQ=255)
                                                                                                                       gCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGca                                      >  2:1088989/1‑135 (MQ=255)
                                                                                                                        acaCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGAc                      >  2:1036136/3‑150 (MQ=255)
                                                                                                                           ctctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACggg                   <  1:1255143/150‑1 (MQ=255)
                                                                                                                              ttCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGaaa                >  2:728062/1‑150 (MQ=255)
                                                                                                                               tCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGCCGGGAAAt               >  1:1569938/1‑150 (MQ=255)
                                                                                                                                ccGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATc              <  1:1091417/150‑1 (MQ=255)
                                                                                                                                ccGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTAAGGAGCAACAACAGCGGATCGACGGGAAATc              >  2:369765/1‑150 (MQ=255)
                                                                                                                                 cGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCa             <  1:1674015/150‑1 (MQ=255)
                                                                                                                                   cTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAATAGAAATACACCCGTTACCGTCAGGCGCAACAACAGCGGATAGACGGGAAAACaaa           >  2:140502/1‑150 (MQ=255)
                                                                                                                                    tataCCTCCACAGAGCTCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTAAATGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCaaaa          <  1:1227793/150‑1 (MQ=255)
                                                                                                                                      taCCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAgc        >  2:617043/1‑150 (MQ=255)
                                                                                                                                        ccTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAgcgc      >  2:39664/1‑150 (MQ=255)
                                                                                                                                         cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     <  1:1823401/150‑1 (MQ=255)
                                                                                                                                           ccACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTgg   >  1:432924/1‑150 (MQ=255)
                                                                                                                                            cacaGAGCGCGCTAAACGCACCCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGt  >  2:1221725/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457967‑3458256

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