Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,432,994 (C)6→5 coding (567/1125 nt) smf ← protein Smf

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,432,9890C.99.1% 442.7 / ‑2.5 107coding (572/1125 nt)smfprotein Smf
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (53/53);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 115 total aligned reads displayed.

TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCA  >  NC_000913/3432850‑3433123
                                                                                                                                           |                                                                                                                                      
tGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcgc                                                                                                                                        >  2:510632/1‑140 (MQ=255)
  cTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccc                                                                                                                                      <  2:1126568/140‑1 (MQ=255)
  cTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccc                                                                                                                                      <  1:519198/140‑1 (MQ=255)
   tGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGcccc                                                                                                                                     <  2:583654/140‑1 (MQ=255)
      cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccccctg                                                                                                                                  >  1:1181826/1‑138 (MQ=255)
         cAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTc                                                                                                                              <  2:288625/140‑1 (MQ=255)
            cGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAg                                                                                                                           >  1:1542198/1‑140 (MQ=255)
                   cTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTg                                                                                                                    >  2:335352/1‑140 (MQ=255)
                   cTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTg                                                                                                                    >  2:572921/1‑140 (MQ=255)
                         aaaGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCACGGGGTGGAACATCGAGGGGAAATCCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGcc                                                                                                              <  2:1139742/140‑1 (MQ=255)
                            gCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGt                                                                                                           <  2:191282/140‑1 (MQ=255)
                            gCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGt                                                                                                           <  2:1399422/140‑1 (MQ=255)
                             ccGCTTCCACCACCAGTACACCTTTACTTAAACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTc                                                                                                          >  2:1660302/1‑140 (MQ=255)
                              cGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                                                         >  2:507760/1‑140 (MQ=255)
                                 ttCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGc                                                                                                      >  1:1073670/1‑140 (MQ=255)
                                 ttCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGc                                                                                                      >  1:1643199/1‑140 (MQ=255)
                                       caccaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcg                                                                                                <  1:1454931/140‑1 (MQ=255)
                                        accaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                                               <  1:507760/140‑1 (MQ=255)
                                            gTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCC‑‑GTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                                          <  2:495459/140‑1 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                                          <  1:335352/140‑1 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                                          >  2:1566364/1‑140 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                                          <  1:572921/140‑1 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                                          <  2:1542198/140‑1 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                                          <  2:170411/140‑1 (MQ=255)
                                               cacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcg                                                                                                >  2:1574872/1‑132 (MQ=255)
                                               cacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcg                                                                                                <  1:692922/132‑1 (MQ=255)
                                               cacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcg                                                                                                >  2:692922/1‑132 (MQ=255)
                                                   tttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                                    <  2:853857/140‑1 (MQ=255)
                                                   tttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                                    <  2:1008163/140‑1 (MQ=255)
                                                   tttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                                    <  1:1652851/140‑1 (MQ=255)
                                                    ttACTTAGACCACTTATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaagg                                                                                   >  2:210973/1‑138 (MQ=255)
                                                       cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTa                                                                                >  1:657567/1‑140 (MQ=255)
                                                       cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTa                                                                                >  1:451431/1‑140 (MQ=255)
                                                       cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTa                                                                                >  1:487534/1‑140 (MQ=255)
                                                          aGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTa                                                                             >  2:882751/1‑140 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                                      >  1:1584148/1‑129 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAg                                                                           >  1:1283119/1‑140 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAg                                                                           >  1:308696/1‑140 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAg                                                                           >  2:1299434/1‑140 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAg                                                                           >  1:1486951/1‑140 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAg                                                                           >  1:693451/1‑140 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGGATTAAg                                                                           >  1:392056/1‑140 (MQ=255)
                                                              cACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                                         <  2:392056/140‑1 (MQ=255)
                                                              cACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                                         <  2:541634/140‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       >  1:1528495/1‑140 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       <  1:1660302/140‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       <  2:451431/140‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       <  1:1474198/140‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       <  2:487534/140‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       >  1:236820/1‑140 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       <  2:657567/140‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       >  1:162493/1‑140 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       <  2:162493/140‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       >  1:949419/1‑140 (MQ=255)
                                                                  gATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCAtt                                                                     >  1:442993/1‑140 (MQ=255)
                                                                   aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                                                                    >  2:805151/1‑140 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGt                                                                                 <  1:98187/122‑1 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                                               <  2:1528495/140‑1 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGAAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                                               >  2:209285/1‑140 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAATTTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                                  >  1:543930/1‑134 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTaa                                                                            <  2:290654/124‑1 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTaa                                                                            >  1:290654/1‑124 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       <  2:492990/129‑1 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                       >  1:492990/1‑129 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATAc                                                            >  1:1440861/1‑140 (MQ=255)
                                                                             ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAg                                                          >  2:1334035/1‑140 (MQ=255)
                                                                             ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAg                                                          <  2:1283119/140‑1 (MQ=255)
                                                                              tcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGc                                                         <  1:805151/140‑1 (MQ=255)
                                                                                     gggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                                 <  1:711193/140‑1 (MQ=255)
                                                                                      ggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                                 <  2:308696/140‑1 (MQ=255)
                                                                                      ggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                                 <  1:1459123/140‑1 (MQ=255)
                                                                                      ggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                                 <  2:942224/140‑1 (MQ=255)
                                                                                      ggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                                 <  2:949419/140‑1 (MQ=255)
                                                                                       gAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGAc                                                <  2:618695/140‑1 (MQ=255)
                                                                                        aaaTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                                               <  2:1643199/140‑1 (MQ=255)
                                                                                             gTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAtt                                          <  2:1073670/140‑1 (MQ=255)
                                                                                             gTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAtt                                          <  1:273414/140‑1 (MQ=255)
                                                                                                    aGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGt                                   >  2:771685/1‑140 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTa                                  >  2:537919/1‑140 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTa                                  >  1:60790/1‑140 (MQ=255)
                                                                                                      gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGGATGGGTTTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTaa                                 >  2:683876/1‑140 (MQ=255)
                                                                                                      gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                                 >  2:1582307/1‑124 (MQ=255)
                                                                                                      gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                                 <  1:1582307/124‑1 (MQ=255)
                                                                                                      gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                                 <  1:468884/124‑1 (MQ=255)
                                                                                                        gTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCa                                            >  1:1341818/1‑127 (MQ=255)
                                                                                                        gTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCa                                            <  2:1341818/127‑1 (MQ=255)
                                                                                                          ggAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgct                             <  1:1566364/140‑1 (MQ=255)
                                                                                                           gAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgctg                            >  1:127666/1‑140 (MQ=255)
                                                                                                            aaCATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgctgc                           >  1:609902/1‑140 (MQ=255)
                                                                                                               aTCGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt                        <  1:537919/140‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgct                             <  1:1140039/133‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgct                             >  2:1140039/1‑133 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATg                     <  1:210973/140‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAAGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTAt                      >  2:1487611/1‑140 (MQ=255)
                                                                                                                      ggAAATTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCt                 <  1:766539/140‑1 (MQ=255)
                                                                                                                        aaaTTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTac               <  1:771685/140‑1 (MQ=255)
                                                                                                                         aaTTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTac               <  2:1445741/139‑1 (MQ=255)
                                                                                                                          aTTCCGAGACGAGAGCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTacac             >  2:729194/1‑140 (MQ=255)
                                                                                                                                     agagCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCa  <  2:1486951/140‑1 (MQ=255)
                                                                                                                                     agagCGCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAAGCTGCTTTATGCGCTACACCATCGATTccc  >  1:251405/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                      
TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCA  >  NC_000913/3432850‑3433123

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