Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 393,962 C→A 41.3% pseudogene (993/1449 nt) yaiT → pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913393,9620CA41.3% 228.9 / 269.9 208pseudogene (993/1449 nt)yaiTpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base C (71/51);  new base A (44/42);  total (115/93)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

TCAGGTGGATCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGAC  >  NC_000913/393815‑394109
                                                                                                                                                   |                                                                                                                                                   
  aGGTGGATCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAAcc                                                                                                                                                  >  1:484442/1‑149 (MQ=255)
   ggTGGATCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCa                                                                                                                                                 >  4:362906/1‑149 (MQ=255)
    gTGGATCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCaa                                                                                                                                                >  4:1437240/1‑149 (MQ=255)
    gTGGATCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTACAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCaa                                                                                                                                                >  4:964391/1‑149 (MQ=255)
       gATCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAATTc                                                                                                                                             >  3:1309680/1‑149 (MQ=255)
         tCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCg                                                                                                                                           >  1:1392530/1‑149 (MQ=255)
         tCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCg                                                                                                                                           >  2:1065899/1‑149 (MQ=255)
         tCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCg                                                                                                                                           >  3:84468/1‑149 (MQ=255)
           aaGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGAc                                                                                                                                         <  1:674468/149‑1 (MQ=255)
           aaGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGAc                                                                                                                                         >  1:341252/1‑149 (MQ=255)
            aGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACg                                                                                                                                        >  1:165260/1‑149 (MQ=255)
              aCAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCt                                                                                                                                      >  4:485406/1‑149 (MQ=255)
              aCAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCt                                                                                                                                      >  2:1199202/1‑149 (MQ=255)
                 gCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTgg                                                                                                                                   <  2:5679/149‑1 (MQ=255)
                  ccAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTggg                                                                                                                                  <  4:1450682/149‑1 (MQ=255)
                        tGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCg                                                                                                                            >  2:82272/1‑149 (MQ=255)
                          ccTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGt                                                                                                                          >  3:43576/1‑149 (MQ=255)
                                tCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGaaa                                                                                                                    <  1:833487/149‑1 (MQ=255)
                                 cTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAg                                                                                                                   >  3:412155/1‑149 (MQ=255)
                                  tACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGc                                                                                                                  >  4:559275/1‑149 (MQ=255)
                                  tACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGc                                                                                                                  >  1:184677/1‑149 (MQ=255)
                                  tACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGc                                                                                                                  >  1:1461083/1‑149 (MQ=255)
                                  tACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGc                                                                                                                  >  4:941326/1‑149 (MQ=255)
                                    ctctTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCg                                                                                                                >  3:1329942/1‑149 (MQ=255)
                                    ctctTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCg                                                                                                                <  3:1194828/149‑1 (MQ=255)
                                    ctctTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCg                                                                                                                <  3:1341766/149‑1 (MQ=255)
                                    ctctTCTGCTCACTTCGCGAACAACGTCTTTGTTGATAGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCg                                                                                                                <  3:964391/149‑1 (MQ=255)
                                     tcttctGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCgg                                                                                                               >  1:391244/1‑149 (MQ=255)
                                      cttctGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGt                                                                                                              <  2:207745/149‑1 (MQ=255)
                                       ttctGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGtt                                                                                                             >  2:1114809/1‑149 (MQ=255)
                                       ttctGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGtc                                                                                                             <  2:1281016/149‑2 (MQ=255)
                                          tGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACg                                                                                                          >  2:707935/1‑149 (MQ=255)
                                          tGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTCGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGTTATCACCAACTCCGACGCTTGGGGCTTCGGTATGAAAGCCGGTTGCg                                                                                                          >  2:963107/1‑149 (MQ=255)
                                           gCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGa                                                                                                         <  1:1110155/149‑1 (MQ=255)
                                               aCTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTgg                                                                                                                                   <  1:741693/119‑1 (MQ=255)
                                               aCTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTc                                                                                                     >  3:212546/1‑149 (MQ=255)
                                                 ttCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCaa                                                                                                   >  3:395019/1‑149 (MQ=255)
                                                      gAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTgg                                                                                              <  2:256791/149‑1 (MQ=255)
                                                      gAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTgg                                                                                              <  1:336509/149‑1 (MQ=255)
                                                      gAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTgg                                                                                              <  2:1218303/149‑1 (MQ=255)
                                                      gAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTgg                                                                                              <  4:80492/149‑1 (MQ=255)
                                                      gAACAACGTCTTTGTTGATGGTAGCTGGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTgg                                                                                              <  4:792549/149‑1 (MQ=255)
                                                          aacGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGa                                                                                          >  2:169498/1‑149 (MQ=255)
                                                             gTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTgg                                                                                              <  1:1215669/142‑1 (MQ=255)
                                                                tttGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTgg                                                                                    >  4:866332/1‑149 (MQ=255)
                                                                tttGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTgg                                                                                    >  2:731986/1‑149 (MQ=255)
                                                                    ttgATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTAc                                                                                >  1:1144683/1‑149 (MQ=255)
                                                                    ttgATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTAc                                                                                >  1:224298/1‑149 (MQ=255)
                                                                     tgATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACg                                                                               <  1:819670/149‑1 (MQ=255)
                                                                      gATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGt                                                                              <  3:60983/149‑1 (MQ=255)
                                                                        tGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCg                                                                                                                >  2:1206697/1‑113 (MQ=255)
                                                                         agTAGCTCAAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGAc                                                                           >  3:1059600/2‑149 (MQ=255)
                                                                             gCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCt                                                                       >  3:1174933/1‑149 (MQ=255)
                                                                                tGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTAc                                                                    >  3:712023/1‑149 (MQ=255)
                                                                                    cTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCa                                                                >  1:104369/1‑149 (MQ=255)
                                                                                        tctcACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGt                                                            >  1:853866/1‑149 (MQ=255)
                                                                                            aCTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTgg                                                                                              <  1:1411476/111‑1 (MQ=255)
                                                                                            aCTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTgg                                                                                              >  2:1411476/1‑111 (MQ=255)
                                                                                             cTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGa                                                                                                         >  2:603349/1‑99 (MQ=255)
                                                                                             cTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGa                                                                                                         <  1:603349/99‑1 (MQ=255)
                                                                                             cTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTg                                                       <  2:1323600/149‑1 (MQ=255)
                                                                                             cTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTg                                                       <  2:853866/149‑1 (MQ=255)
                                                                                             cTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAACCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTg                                                       >  1:1079781/1‑149 (MQ=255)
                                                                                              ttCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTgg                                                      >  4:1169905/1‑149 (MQ=255)
                                                                                              ttCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTgg                                                      >  3:282432/1‑149 (MQ=255)
                                                                                                caacaaCGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTc                                                    <  1:1003082/149‑1 (MQ=255)
                                                                                                 aacaacGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCt                                                   >  1:524378/1‑149 (MQ=255)
                                                                                                  acaacGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTTCGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTg                                                  >  3:1405481/1‑149 (MQ=255)
                                                                                                    aacGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAAc                                                                                                 >  2:33899/1‑100 (MQ=255)
                                                                                                    aacGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGtt                                                <  3:1470399/149‑1 (MQ=255)
                                                                                                    aacGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGtt                                                >  3:345159/1‑149 (MQ=255)
                                                                                                    aacGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGtt                                                >  2:54413/1‑149 (MQ=255)
                                                                                                       gACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCa                                             <  4:712023/149‑1 (MQ=255)
                                                                                                            gTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTc                                               >  2:746593/1‑142 (MQ=255)
                                                                                                             tCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTgg                                       <  1:940509/149‑1 (MQ=255)
                                                                                                              cTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGCGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGt                                      <  1:974999/149‑1 (MQ=255)
                                                                                                              cTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGt                                      <  1:446387/149‑1 (MQ=255)
                                                                                                               tGCAACCATGAGCAACGGTACTTACGATGACGGTAGCAGCAACTCCGACGCTTGGGGCTTCGGTTTAAAAGCCGGTTACGACTTCAAACTGGGTGATGGTGGTTACGTAACTCCTTACGGCAGCGTTTCTGGTCTGTTCCCGTCTGGtg                                     >  3:921023/1‑149 (MQ=255)
                                                                                                                gCAACCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGtga                                    >  2:478900/1‑149 (MQ=255)
                                                                                                                gCAACCATGAGCAAAGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGtga                                    >  3:112908/1‑149 (MQ=255)
                                                                                                                   aCCATGAGCAACGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGtgatga                                 <  4:1029084/149‑1 (MQ=255)
                                                                                                                    ccATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGAc                                >  4:60246/1‑149 (MQ=255)
                                                                                                                     cATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACt                               >  4:101414/1‑149 (MQ=255)
                                                                                                                     cATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACt                               >  1:877979/1‑149 (MQ=255)
                                                                                                                         aGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCa                           <  4:1059600/149‑1 (MQ=255)
                                                                                                                            aaCGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCt                        <  2:1427047/149‑1 (MQ=255)
                                                                                                                             aCGGTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTg                       >  3:527766/1‑149 (MQ=255)
                                                                                                                              cGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGa                      <  1:1426419/149‑1 (MQ=255)
                                                                                                                                gTACTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCCGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGc                    <  4:1405481/149‑1 (MQ=255)
                                                                                                                                   cTTACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAAc                 <  4:43576/149‑1 (MQ=255)
                                                                                                                                    ttACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCg                                                                                                                <  4:1321508/53‑1 (MQ=255)
                                                                                                                                    ttACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTAc                                                                                                           >  2:941069/1‑58 (MQ=255)
                                                                                                                                    ttACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGtgat                                   <  3:775285/130‑1 (MQ=255)
                                                                                                                                    ttACGTTGACGGTAGAACCAACTCCGACGCTTGGGGCTTCGGt                                                                                                                          >  3:428166/1‑43 (MQ=255)
                                                                                                                                      aCGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGAc              >  4:1309582/1‑149 (MQ=255)
                                                                                                                                           gACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGaa         <  3:101414/149‑1 (MQ=255)
                                                                                                                                             cGGTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGAGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGAAGACTACCAGCTGAGCAACGACATGAAAg       <  2:1079781/149‑1 (MQ=255)
                                                                                                                                               gTAGAACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGtt     >  2:456498/1‑149 (MQ=255)
                                                                                                                                                 aGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGa   <  4:412155/149‑1 (MQ=255)
                                                                                                                                                  gCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGAc  <  2:1392530/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCAGGTGGATCAAGACAGCCAGACTGCCTACATCTACTCTTCTGCTCACTTCGCGAACAACGTCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGAC  >  NC_000913/393815‑394109

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