Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,465,130 T→C 6.4% V17V (GTT→GTC gspM → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,465,1300TC 6.4% 263.9 / 5.2 95V17V (GTT→GTCgspMgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (39/48);  new base C (3/3);  total (43/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

TCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGT  >  NC_000913/3464983‑3465277
                                                                                                                                                   |                                                                                                                                                   
tcgatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAACGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTg                                                                                                                                                    <  1:1076598/149‑1 (MQ=255)
tcgatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTg                                                                                                                                                    <  1:464024/149‑1 (MQ=255)
  gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGTAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCg                                                                                                                                                  >  2:1433732/1‑149 (MQ=255)
  gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCg                                                                                                                                                  >  1:1456229/1‑149 (MQ=255)
   atcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCgg                                                                                                                                                 >  1:225444/1‑149 (MQ=255)
      gTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGccc                                                                                                                                              >  2:1203849/1‑147 (MQ=255)
         tCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTgg                                                                                                                                           <  1:622430/149‑1 (MQ=255)
           tgtgAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCg                                                                                                                                         <  1:1437128/149‑1 (MQ=255)
           tgtgAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCg                                                                                                                                         <  2:1145693/149‑1 (MQ=255)
            gtgAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTTGGCAGAAAAATCCACATCGGAAAAACAGGTCGTTGCGGCACTGGCgg                                                                                                                                        >  2:215498/1‑149 (MQ=255)
             tgAGTTAACACAGTCCTGGCTTCCGATGGAAAAAAAAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTTCGGCACTGGCGGt                                                                                                                                       <  2:95546/149‑1 (MQ=255)
                gTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAGAAAGATCGTTACGTCACTGGCGGGtct                                                                                                                                    >  1:537619/1‑149 (MQ=255)
                gTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTtct                                                                                                                                    >  2:637422/1‑149 (MQ=255)
                 ttAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTtctc                                                                                                                                   >  2:1367385/1‑149 (MQ=255)
                 ttAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTtctc                                                                                                                                   >  1:448652/1‑149 (MQ=255)
                     cacaGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTctcgc                                                                                                                               <  1:291120/149‑1 (MQ=255)
                     cacaGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTctcgc                                                                                                                               >  2:1412862/1‑149 (MQ=255)
                        aGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTctcgctcg                                                                                                                            <  1:1262570/149‑1 (MQ=255)
                            cTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTCTGGCCAGTAAGGAACTCAGGGACATGATTACATCCTGGTTGCCACGATCACCCACACTGGACAAACAAACCGGGCCGGCACTGGCGGCTCCCCCGCCCGGCgg                                                                                                                        >  2:849488/1‑148 (MQ=255)
                            cTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGt                                                                                                                        >  1:514658/1‑149 (MQ=255)
                              ggCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGttt                                                                                                                      <  1:525129/149‑1 (MQ=255)
                              ggCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGttt                                                                                                                      >  2:1179393/1‑149 (MQ=255)
                                 ttCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGtttttt                                                                                                                   <  1:1365634/149‑1 (MQ=255)
                                  tCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGATAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTttg                                                                                                                  <  1:50711/149‑1 (MQ=255)
                                     gATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTGCCCTCGCTCGGCGTTTTttgttg                                                                                                               >  1:746680/1‑149 (MQ=255)
                                      aTGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTCGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTgg                                                                                                              >  2:776626/1‑149 (MQ=255)
                                      aTGGAAAAAACAGAAAAAGATCCGGTCAGCGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGTTCGTCGCGGCACTGGCGGTTCTCTGGCTCGGCGTTTTTTGTTgg                                                                                                              >  1:1025667/1‑149 (MQ=255)
                                      aTGGAAAAAACAGAAAAAGATCCGGTCAGAGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTCGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTgg                                                                                                              >  1:1025730/1‑149 (MQ=255)
                                         gAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCtt                                                                                                           <  2:201259/149‑1 (MQ=255)
                                           aaaaaCAGAAAAAGATCCGGTCAGTGGCGGATGGACAGTAAGGAACTCAGGGAAAAGATCAAACCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGCCGGTTCGCTCGCCCGGCGTTTTTCGTTGGCCTgg                                                                                                         >  2:1017191/1‑149 (MQ=255)
                                            aaaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTggg                                                                                                        <  1:617241/149‑1 (MQ=255)
                                             caaCAGAAAAAGATCCTGTCAGTGGTGTATGGACAGCAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTgggg                                                                                                       <  1:116851/148‑1 (MQ=255)
                                             aaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACAAGTAGAGGTGTGTCTCGGCGTTTTTTGTTGGCTTgggg                                                                                                       <  1:1073943/149‑1 (MQ=255)
                                              aaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGGGGGCGGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTATCTCGCTCGGCGTTTTTGGTTGGCTTGGGGt                                                                                                      >  1:1430751/1‑149 (MQ=255)
                                                cAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGa                                                                                                    <  2:225444/149‑1 (MQ=255)
                                                 aGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGAt                                                                                                   >  1:1203270/1‑149 (MQ=255)
                                                   aaaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATaa                                                                                                 >  1:1405296/1‑149 (MQ=255)
                                                   aaaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATaa                                                                                                 >  1:161542/1‑149 (MQ=255)
                                                   aaaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATaa                                                                                                 >  1:560721/1‑149 (MQ=255)
                                                   aaaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGGTCTCCCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATaa                                                                                                 >  1:603124/1‑149 (MQ=255)
                                                         acccGGTCAGTGGGGGATGGACGGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCga                                                                                           <  2:179128/147‑1 (MQ=255)
                                                         aTCCGGTCACTGGTGGATGGACAGTAAGGAACTCAGGGAAATGAGTAAATCATGGTGGGCAGAAAACTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCga                                                                                           <  2:1421476/149‑1 (MQ=255)
                                                          tCCGGGCAGGGGTGTATGGACAGTAAGGAACTCAGGGAATGATTAAATCATGGTGGGCATAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCgata                                                                                         <  2:1047335/149‑1 (MQ=255)
                                                           cgggTCAGGGGTGTATGGACAGGAAGGAACCCAGGGAAAGGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCgata                                                                                         <  2:361274/147‑1 (MQ=255)
                                                           ccGGTCAGTGTTGTAGGGACTGTAAGGAACTCAGGGAAATTATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTGGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCgata                                                                                         <  2:1128220/149‑1 (MQ=255)
                                                           ccGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCgata                                                                                         <  2:28956/149‑1 (MQ=255)
                                                           ccGGTCAGTGGTGGATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCgata                                                                                         <  1:148740/149‑1 (MQ=255)
                                                             ggTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCgataga                                                                                       <  1:1351491/149‑1 (MQ=255)
                                                                cAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATAc                                                                                    <  1:1367385/149‑1 (MQ=255)
                                                                     gtgtATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACa                                                                               <  1:1005023/149‑1 (MQ=255)
                                                                       gtATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGACACCAACAtt                                                                             >  1:29974/1‑149 (MQ=255)
                                                                          tGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATCAAACCGATAGAGACCTACATTgga                                                                          >  2:853065/1‑147 (MQ=255)
                                                                          tGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCa                                                                          <  2:514658/149‑1 (MQ=255)
                                                                          aggACCGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCa                                                                          <  2:759058/148‑1 (MQ=255)
                                                                                gTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCAt                                                                    <  1:1351554/149‑1 (MQ=255)
                                                                                  aaGGTACTGAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAACCAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGTCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCa                                                                  <  2:197668/149‑1 (MQ=255)
                                                                                  aaGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCGCGCTCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCa                                                                  <  1:408863/149‑1 (MQ=255)
                                                                                      aaCTCAGGGAAATGATTAAAGCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTcc                                                              <  2:572314/149‑1 (MQ=255)
                                                                                         tCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTTTTTTGTGGAGATAAAACCGATAGACACCTACCCTGCCGAGGAACAGTCCCAc                                                           >  2:1059482/1‑149 (MQ=255)
                                                                                         tCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCAc                                                           >  1:1226735/1‑149 (MQ=255)
                                                                                         tCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCAc                                                           >  2:628528/1‑149 (MQ=255)
                                                                                           agagAAATGATTAAATCATGGTGGGGAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCAcgc                                                         <  2:431107/146‑1 (MQ=255)
                                                                                           aGGGAGATGATTAAAGCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCAcgc                                                         <  2:1282986/149‑1 (MQ=255)
                                                                                             ggAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACgcgc                                                       >  1:1004623/1‑149 (MQ=255)
                                                                                                 aTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCaaaaa                                                  <  1:1186798/149‑1 (MQ=255)
                                                                                                   gATTAAATCAGGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCaaaaaa                                                 <  2:168974/149‑1 (MQ=255)
                                                                                                    aTTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGGGCATCAGTCCCACGCGCgaaaaaa                                                >  1:55195/1‑148 (MQ=255)
                                                                                                    aTTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAAt                                                >  1:580066/1‑149 (MQ=255)
                                                                                                     ttAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCCCGCGCAAAAAAtt                                               >  1:221459/1‑149 (MQ=255)
                                                                                                      tAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTa                                              >  1:425875/1‑149 (MQ=255)
                                                                                                          tCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTaagaa                                          <  2:1456229/149‑1 (MQ=255)
                                                                                                            aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGTTCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAgg                                        <  2:459446/149‑1 (MQ=255)
                                                                                                               gtggGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACa                                     >  2:554271/1‑149 (MQ=255)
                                                                                                               gtggGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAAACGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAACAAACTAAGAAGGACa                                     >  2:843760/1‑149 (MQ=255)
                                                                                                                    cAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCaaa                                <  1:116672/149‑1 (MQ=255)
                                                                                                                    cAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTGTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCaaa                                <  2:560721/149‑1 (MQ=255)
                                                                                                                    aaGAAAAATCCATAACGAAAAAACAGAACGTTGCGGGACTGTCGGTACTCTCGCTCGGCGTTTTTTGTTGACTTAGAGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGACACACAAACAAAAAATTAAGAAGGACATCaaa                                <  2:216520/148‑1 (MQ=255)
                                                                                                                       aaaaaTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCGGTCCCACGCGCAAAAAATTAAGAAGGACACCAACTgg                             >  1:544624/1‑149 (MQ=255)
                                                                                                                       aaaaaTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCACATgt                             >  2:868726/1‑148 (MQ=255)
                                                                                                                       aaaaaTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATgg                             >  1:221069/1‑149 (MQ=255)
                                                                                                                            tCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCa                        <  1:628528/149‑1 (MQ=255)
                                                                                                                             ccACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGAATCAGCCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAg                       <  1:926458/149‑1 (MQ=255)
                                                                                                                              cacaTGGGAAAAACGGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAgg                      <  2:73450/149‑1 (MQ=255)
                                                                                                                              cacaTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAgg                      >  2:1393362/1‑149 (MQ=255)
                                                                                                                               acaTCGGAAAAACAGATCGTTGCGGCACTGGCGGATCTCTCGCTCGGCGTTTTTTGTTGACTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCAGCAGTCCCACGCGCGCAAACTTAAGAAGGACATCAACTGGAGGCAGGa                     >  2:375393/1‑149 (MQ=255)
                                                                                                                                 aTCGGAAAAGCAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATc                   <  2:333013/149‑1 (MQ=255)
                                                                                                                                   cGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGAGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCaa                 <  1:1145784/149‑1 (MQ=255)
                                                                                                                                       aaaaCAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGa             <  1:1179393/149‑1 (MQ=255)
                                                                                                                                         agtaGATCGTCGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGc           <  2:1025667/146‑1 (MQ=255)
                                                                                                                                         aaCAGATCGTCGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGc           <  1:776626/149‑1 (MQ=255)
                                                                                                                                         aaCAGATCGTCGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGc           <  2:1025730/149‑1 (MQ=255)
                                                                                                                                           cAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCAc         >  1:288378/1‑149 (MQ=255)
                                                                                                                                               tCGTTGCGGCCCTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGCGTGATAACACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCAt     <  2:1373572/149‑1 (MQ=255)
                                                                                                                                               tCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTTGCTTGGGGTGATAAAACCGATAGCTACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCCGGATCAAGCGAGCACCCAt     >  2:693400/1‑149 (MQ=255)
                                                                                                                                                cGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATg    <  1:1219287/149‑1 (MQ=255)
                                                                                                                                                  ttGGGGCACTGGCGGTTCTCTCGCTCTGCGTTTTTTGTTGGCTTGGGGTGAGAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGt  <  1:1227506/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGT  >  NC_000913/3464983‑3465277

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