Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,198,875 A→T D29V (GAT→GTT)  nfi → endonuclease V

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,198,8750AT98.5% 435.3 / NA 135D29V (GAT→GTT) nfiendonuclease V
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (58/74);  total (60/75)
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 136 total aligned reads displayed.

GGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACAT  >  NC_000913/4198741‑4199010
                                                                                                                                      |                                                                                                                                       
     cAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAGAGATCCGCCGGTTCTGATCGcc                                                                                                                               >  1:473630/1‑140 (MQ=255)
          cATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGc                                                                                                                          <  2:1244595/140‑1 (MQ=255)
                cTGTCTGAACAGTATCATCGCTAAGGAGTGATTATTGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGt                                                                                                                    <  2:520376/140‑1 (MQ=255)
                cTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCCCCGGTTCTGATCGCCGGAGCCGATGt                                                                                                                    <  2:1535863/140‑1 (MQ=255)
                  gTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAATATCCACCGGTTCTGATCGCCGGAGCCGATGTCg                                                                                                                  >  1:30588/1‑140 (MQ=255)
                  gTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAg                                                                                                                           >  2:1133156/1‑131 (MQ=255)
                  gTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGt                                                                                                                    <  2:18767/138‑1 (MQ=255)
                   tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCgg                                                                                                                 >  1:1445234/1‑140 (MQ=255)
                     tGAACAGTATCACCGATAAGGATTGTTTATGGATCTCCCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGt                                                                                                               <  1:734474/140‑1 (MQ=255)
                          aGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAg                                                                                                          <  1:951161/140‑1 (MQ=255)
                          aGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAg                                                                                                          <  1:1342024/140‑1 (MQ=255)
                          aGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAg                                                                                                          <  2:1415902/140‑1 (MQ=255)
                          aGTATCACCGCTAAGGACTGATCATGGATCAGGTGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAg                                                                                                          <  2:748500/140‑1 (MQ=255)
                            tATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCa                                                                                                        <  1:947430/140‑1 (MQ=255)
                            tATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCa                                                                                                        <  1:880763/140‑1 (MQ=255)
                            tATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCa                                                                                                        <  1:1186252/140‑1 (MQ=255)
                            tATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCa                                                                                                        <  2:243390/140‑1 (MQ=255)
                            tATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCa                                                                                                        <  2:37395/140‑1 (MQ=255)
                            tATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATGTGATCGCCGGAGCTAATGTCGGGTTTGAGCa                                                                                                        >  2:697413/1‑140 (MQ=255)
                            tATCACCGCTAAGGAATGATTATGGCTCTTGCGTCATGACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCa                                                                                                        <  1:614128/140‑1 (MQ=255)
                              tCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAgg                                                                                                      >  2:216881/1‑140 (MQ=255)
                              tCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGCTCTAATCGCCGGAGCCGGTGTCGGGTTTGAGCAgg                                                                                                      >  2:614128/1‑140 (MQ=255)
                               cACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAggg                                                                                                     <  2:1445234/140‑1 (MQ=255)
                               cACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAggg                                                                                                     <  2:905426/140‑1 (MQ=255)
                                 ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGt                                                                                                                    >  2:66422/1‑123 (MQ=255)
                                 ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGt                                                                                                                    <  1:66422/123‑1 (MQ=255)
                                 ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCg                                                                                                   <  2:869834/140‑1 (MQ=255)
                                 ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCg                                                                                                   >  2:887396/1‑140 (MQ=255)
                                 ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCg                                                                                                   >  1:308708/1‑140 (MQ=255)
                                 ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTACTGATCGCCGGTGCCGATGTCGGGTTTGAGCAGGGCg                                                                                                   >  1:744581/1‑140 (MQ=255)
                                 ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAAATCCACCGTTTCTGATAGCCGGAGCAGATGTCGGGTTTGAGCAGGGCg                                                                                                   >  2:566497/1‑140 (MQ=255)
                                    cTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTTATCGCCGGAGCAGATGTCGGGGTTGAGCAGGGCGgag                                                                                                >  1:582794/1‑140 (MQ=255)
                                     tAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAg                                                                                                          <  2:186845/129‑1 (MQ=255)
                                     tAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAg                                                                                                          >  1:186845/1‑129 (MQ=255)
                                      aaGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCa                                                                                                        <  1:1530329/130‑1 (MQ=255)
                                      aaGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                              >  1:1068090/1‑140 (MQ=255)
                                      aaGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                              >  2:224766/1‑140 (MQ=255)
                                      aaGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                              >  2:1492014/1‑140 (MQ=255)
                                      aaGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                              >  1:70181/1‑140 (MQ=255)
                                      aaGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                              >  1:48930/1‑140 (MQ=255)
                                      aaGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGAACGCCGGAACCGATGGCGGGTTTGAGCAGGGCGGAGaa                                                                                              >  1:589749/1‑140 (MQ=255)
                                            tGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGCAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGAcg                                                                                        >  1:661178/1‑140 (MQ=255)
                                              aTTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGATCCGATGTAGGGGTTGAGCAGGACGGAGAAGTGAcgcg                                                                                      >  1:664146/1‑140 (MQ=255)
                                              aTTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGAcgcg                                                                                      <  1:297528/140‑1 (MQ=255)
                                              aTTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCGGATCGCCGGAGCGGATCTCGGGTTTAAGCAGGGCGGAGAAGTGAcgcg                                                                                      >  2:402692/1‑140 (MQ=255)
                                                tATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAg                                                                                    <  1:224766/140‑1 (MQ=255)
                                                   ggATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAgcgg                                                                                 >  2:1274155/1‑140 (MQ=255)
                                                   ggATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAgcgg                                                                                 >  1:1486606/1‑140 (MQ=255)
                                                    gATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAgcggc                                                                                >  1:613213/1‑140 (MQ=255)
                                                      tctcGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGa                                                                              >  1:1305728/1‑140 (MQ=255)
                                                        tcGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATg                                                                            >  1:837352/1‑140 (MQ=255)
                                                        tcGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGATTGGGCAGGGCGTAGAAATGACGCGAGCGGCGATg                                                                            >  2:410031/1‑140 (MQ=255)
                                                         cgcgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGa                                                                              >  2:1193847/1‑137 (MQ=255)
                                                         cgcgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGa                                                                              <  2:1223069/137‑1 (MQ=255)
                                                         cgcgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGa                                                                              >  1:1223069/1‑137 (MQ=255)
                                                         cgcgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGa                                                                              <  1:1193847/137‑1 (MQ=255)
                                                         cgcgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATgg                                                                           <  1:1492014/140‑1 (MQ=255)
                                                         cgcgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATgg                                                                           >  2:271290/1‑140 (MQ=255)
                                                           cgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtg                                                                         <  2:1486606/140‑1 (MQ=255)
                                                            gTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgc                                                                        <  2:371579/140‑1 (MQ=255)
                                                             tCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgct                                                                       <  2:1106304/140‑1 (MQ=255)
                                                              cATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgctg                                                                      <  1:887396/140‑1 (MQ=255)
                                                                ttACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTgctgct                                                                    <  2:1517945/140‑1 (MQ=255)
                                                                        tCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgctg                                                                      <  1:22862/130‑1 (MQ=255)
                                                                        tCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgctg                                                                      >  2:22862/1‑130 (MQ=255)
                                                                         caacaaATCGAACTGTCTTCTTCTGCGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATcc                                                           <  1:627084/140‑1 (MQ=255)
                                                                         caacaaATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATcc                                                           <  2:48930/140‑1 (MQ=255)
                                                                              aaTCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGc                                                      <  2:163317/140‑1 (MQ=255)
                                                                              aaTCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGc                                                      <  1:216881/140‑1 (MQ=255)
                                                                                        gCTTCTTCTGTGATCCGAGAGGATCGACTCGCTAAGGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGt                                            <  2:582794/140‑1 (MQ=255)
                                                                                         cttcttCTGTGATCCTCGAGGATCGCCTCCATAAAGACCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTc                                           <  1:402692/140‑1 (MQ=255)
                                                                                         cttcttCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTc                                           <  2:1346437/140‑1 (MQ=255)
                                                                                         cttcttCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTc                                           <  2:193637/140‑1 (MQ=255)
                                                                                           tcttctGTGATCCGCGAGGATCGACTCGATAAATATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGa                                         <  2:30588/140‑1 (MQ=255)
                                                                                                  tgATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTc                                           <  2:916007/131‑1 (MQ=255)
                                                                                                  tgATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTCCTGGTGAAATATCCCTAGCTTGAGCTGGTCGAGTATaaa                                  >  2:482818/1‑140 (MQ=255)
                                                                                                       cgcgAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGc                             >  1:1007977/1‑140 (MQ=255)
                                                                                                        gcgAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGcc                            <  2:70181/140‑1 (MQ=255)
                                                                                                        gcgAGGATCGACTCGATAAAGAACCACCGGTTCTGTTCGTCGGAGCCGCTGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGcc                            <  2:589749/140‑1 (MQ=255)
                                                                                                         cgAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGACGTGACGCGAGCGGCGATGTTGCTGCTGAAATTTCCGTCGCTTGGGCTGGCCGAGTATAAAGTTGccc                           >  2:551299/1‑140 (MQ=255)
                                                                                                         cgAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGccc                           >  2:1016555/1‑140 (MQ=255)
                                                                                                           aGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGc                         <  2:1007977/140‑1 (MQ=255)
                                                                                                            ggATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCa                        <  1:995789/140‑1 (MQ=255)
                                                                                                             gATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCAt                       <  1:1274155/140‑1 (MQ=255)
                                                                                                                cGACTCGATAAAGATCACCCGGTTCTGATCCCCTGCGCCTATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGc                    <  2:744581/140‑1 (MQ=255)
                                                                                                                  aCTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGcca                  <  2:1305728/140‑1 (MQ=255)
                                                                                                                     cGATAAAGATCCACCGGTTCTGTTCGCCGGAGCCGATCTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccacca               <  1:602545/140‑1 (MQ=255)
                                                                                                                      gATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccaccac              >  1:807684/1‑140 (MQ=255)
                                                                                                                       aTAAAGATCCACCGGTTCTGATAGCCGGACCCGATGTTGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccaccacc             <  1:410031/140‑1 (MQ=255)
                                                                                                                       aTAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccaccacc             <  1:536478/140‑1 (MQ=255)
                                                                                                                        tAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccaccacca            >  2:782281/1‑140 (MQ=255)
                                                                                                                        tAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTAAGCTGGTGGAGTATAAAGTTGCCCGCATCGccaccacca            >  2:601963/1‑140 (MQ=255)
                                                                                                                             aTCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGc                         <  1:872719/122‑1 (MQ=255)
                                                                                                                             aTCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGc                         >  2:872719/1‑122 (MQ=255)
                                                                                                                              tCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCtt      >  1:927610/1‑140 (MQ=255)
                                                                                                                               ccACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGTAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTa     <  2:140222/140‑1 (MQ=255)
                                                                                                                               ccACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTa     <  1:1016555/140‑1 (MQ=255)
                                                                                                                                cACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTAc    >  2:56144/1‑140 (MQ=255)
                                                                                                                                  ccGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGTGCTGGTCGGGTAGAAATTTGCCCGCCTCGCGACCACCATGCCTTACAt  >  2:416025/1‑140 (MQ=255)
                                                                                                                                  ccGGTTCTGATCGCCGAAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACAt  <  2:661178/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACAT  >  NC_000913/4198741‑4199010

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