Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,652,368 C→T 5.6% pseudogene (818/1158 nt) intQ → pseudogene, Qin prophage; phage integrase family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,652,3680CT 5.6% 211.5 / 2.4 73pseudogene (818/1158 nt)intQpseudogene, Qin prophage; phage integrase family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base T (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

CGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAC  >  NC_000913/1652227‑1652509
                                                                                                                                             |                                                                                                                                             
cGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                                                                                                                                        <  2:422186/149‑1 (MQ=255)
 gggAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                                                                                                                                       <  2:954904/149‑1 (MQ=255)
   gAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttttt                                                                                                                                     <  2:516602/149‑1 (MQ=255)
     aCCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAGGTGCATTTTTGTTTTTct                                                                                                                                   <  2:207832/149‑2 (MQ=255)
     aCCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAGGTGCATTTTTGTTTTTct                                                                                                                                   <  2:300003/149‑2 (MQ=255)
       cAATCGGGTCATTCCTCTTATTAAGCCAGAAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAc                                                                                                                                 <  2:1182170/149‑1 (MQ=255)
       cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAc                                                                                                                                 <  1:630034/149‑1 (MQ=255)
       cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAc                                                                                                                                 >  1:541035/1‑149 (MQ=255)
         aTCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                                                                                                               <  2:894208/149‑1 (MQ=255)
         aTCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                                                                                                               >  1:621565/1‑149 (MQ=255)
          tCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                                                                                                              >  2:752345/1‑149 (MQ=255)
             ggTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGAATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAg                                                                                                                           >  2:310742/1‑149 (MQ=255)
                 attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                                                                                                                       <  1:759852/149‑1 (MQ=255)
                 attaatCTTAATAAGCAAGCAATCGACGCACACCGGAGTCAGATCACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                                                                                                                       >  1:217935/1‑149 (MQ=255)
                    catCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc                                                                                                                    <  1:407228/149‑1 (MQ=255)
                     atCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc                                                                                                                    <  1:247484/148‑1 (MQ=255)
                     atCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc                                                                                                                    >  2:247484/1‑148 (MQ=255)
                      acTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCga                                                                                                                  <  1:704801/148‑1 (MQ=255)
                       cTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTATTCAACCTGAAGTGTCAGCgag                                                                                                                 >  1:95046/1‑149 (MQ=255)
                              aGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                                                                                                          <  1:92485/149‑1 (MQ=255)
                                  aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAtt                                                                                                      >  2:333730/1‑149 (MQ=255)
                                   gCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATTGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTa                                                                                                     <  1:846213/149‑1 (MQ=255)
                                   gCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTa                                                                                                     <  2:541035/149‑1 (MQ=255)
                                              tctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtt                                                                                          >  2:365361/1‑149 (MQ=255)
                                                  cGGAGTAAGATGAAATTAACGTGACTGATAAAAGAGCATATCATTTATGTTCACTTCAGAGAGTATTGCAGAACAGAAAATCAAAAATGCACCTTTGTTTTTCAACCTGAAGTTTCAGCGAGAGTAAAAAATTATGGTGACCATTCTAc                                                                                      <  2:543508/149‑1 (MQ=255)
                                                   gcaGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAcc                                                                                     <  2:141510/147‑1 (MQ=255)
                                                     aGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGt                                                                                   <  1:752345/149‑1 (MQ=255)
                                                      gTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGtt                                                                                  <  2:1991/149‑1 (MQ=255)
                                                       tCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGttt                                                                                 >  1:637465/1‑148 (MQ=255)
                                                         aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGAc                                                                               >  1:1166510/1‑149 (MQ=255)
                                                          gATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACt                                                                              <  2:18097/149‑1 (MQ=255)
                                                          gATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACATGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAACATTTTACCGTTGACt                                                                              >  1:1052057/1‑149 (MQ=255)
                                                                aTTAACGAGACTGAGCAAAGAGCATATCATTGATGGTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAATATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGTTGACCATTTTAACCTTGACACAaaaa                                                                        >  1:1089960/1‑149 (MQ=255)
                                                                 ttAACGATACTGAGCAAAGAGCATTTCATTGATTTTAACTTCAGATAGTATGGCAGAACAGAAAAACAAAAATTCAGCTTTGTTTTTCAACGTGAAGTGTCAGCGTGAGTAAAAAATTATGTTGACCATTTTAACGTTGACTCAATAAg                                                                       <  1:344134/149‑1 (MQ=255)
                                                                        gaCTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGAtg                                                                >  1:348495/1‑149 (MQ=255)
                                                                         aCTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGAtgt                                                               <  1:426650/149‑1 (MQ=255)
                                                                             atcAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGa                                                           <  2:188844/147‑1 (MQ=255)
                                                                             aGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACATTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTggta                                                           >  1:128543/1‑147 (MQ=255)
                                                                             aGCAAAGAGCAGATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAATCAAAAATGCACCTATGTTTTTCCACCTGAAGTGTGCGCGAGAGTAAAAAATTATGGTGACCAGTTTACCGTTGACTCAATAAGGCACATGTgggc                                                           >  2:500639/1‑148 (MQ=255)
                                                                               cAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTAACCGTTGACTCAATAAGGCAGATGTGGGACg                                                         <  2:217935/149‑1 (MQ=255)
                                                                                aaaGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGc                                                        >  1:55758/1‑149 (MQ=255)
                                                                                  agcgCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAg                                                      <  2:381532/146‑1 (MQ=255)
                                                                                     gCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCCGAAGGGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGa                                                   <  2:73681/149‑1 (MQ=255)
                                                                                         atCATTGATGATCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTTTTTTTCACCCTAAAGTTTCAGCGAGAGTAATAAATTATGGTGAAAATTTAAAAATTAACTAAATAATGCATAATTGGGACACAGAGATaaa                                               >  2:1059743/1‑149 (MQ=255)
                                                                                          tCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaaaa                                              <  1:310742/149‑1 (MQ=255)
                                                                                          tCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaaaa                                              <  2:1056602/149‑1 (MQ=255)
                                                                                          tCATTGATGTTCACTTCAGACAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGGAAGCGAGAGTAAAAAATTTTTGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaaaa                                              <  1:994371/149‑1 (MQ=255)
                                                                                            aTTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAAct                                            >  1:531241/1‑148 (MQ=255)
                                                                                            aTTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACg                                            >  2:236454/1‑149 (MQ=255)
                                                                                                    tCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTTTGGTGACCATTTTACCGTTGAATCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGAc                                    >  2:90092/1‑149 (MQ=255)
                                                                                                    tCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGa                                                   >  1:549614/1‑134 (MQ=255)
                                                                                                    tCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGa                                                   <  2:549614/134‑1 (MQ=255)
                                                                                                           agagagTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCACCCTGAAGTGTCAGCGAGAGTAAAAAATTATTGTGACCATTTTACCGTTGCCTCAATAAGGCAGATGTGGGACGAAGCGATAAAACGTGCCGGACTCCGCCa                             >  2:1045848/1‑149 (MQ=255)
                                                                                                             agagTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                           >  1:1113161/1‑149 (MQ=255)
                                                                                                                gTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGCGTAAAAAATTATGGTGACCATTTTACCGTTGCCTCAATAATGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaa                        <  1:701613/149‑1 (MQ=255)
                                                                                                                 tATGGCAGAACAGAAAAACAAAACTGCACCTTTGTTTTTCAACCTGAAGCGGCAGAGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATTTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaaa                       <  2:135420/149‑1 (MQ=255)
                                                                                                                 tATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAc                                                                                      >  1:1095936/1‑86 (MQ=255)
                                                                                                                 tATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAc                                                                                      <  2:1095936/86‑1 (MQ=255)
                                                                                                                   tGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATACGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAAt                     <  1:751791/149‑1 (MQ=255)
                                                                                                                   tGGCAGAACAGAAAAACAAAAAAGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTTACAATTTTACCGTTGACTCAATAAGTCAGATTTGGGACGCAGCGATAAAAAGTGCCGGACTACGCCATCGAAAAt                     >  2:784345/1‑149 (MQ=255)
                                                                                                                    ggCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATc                    <  2:621565/149‑1 (MQ=255)
                                                                                                                    ggCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAAAGTGACGGACTGCGCCATCGAAAAta                    >  1:556203/1‑148 (MQ=255)
                                                                                                                      cagaacagaaAAACAAAAATGCACCTTTGTTTTTCATCCTGAAGTGACAGCGACAGTAAAAAAGAATGGTGACCATTTTACCGTTGACTCAATAACGCAGATGTGCGACGCAGCGATAACACGTGCCGGACTCCGCCATCGAAaataat                  <  2:443085/149‑1 (MQ=255)
                                                                                                                      cagaacagaaAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGTACGCAGCGTTAAAACGTGCCGTACTCCGCCATCGAAAATCat                  >  1:766037/1‑149 (MQ=255)
                                                                                                                      aagaacagaaAAACAAAAATGCACATTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCat                  <  2:1023146/148‑1 (MQ=255)
                                                                                                                        gaacagaaAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCataa                >  1:597237/1‑148 (MQ=255)
                                                                                                                           ttgaaAAACAAAGGTGCATTTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAg             >  1:279636/3‑149 (MQ=255)
                                                                                                                           ttgaaAAACAAAGGTGCATTTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAg             >  1:932772/3‑149 (MQ=255)
                                                                                                                              aaaaaCAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCg          >  1:915429/1‑149 (MQ=255)
                                                                                                                              aaaaaCAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACAATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCg          >  2:674299/1‑149 (MQ=255)
                                                                                                                               aaaaCAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCga         <  1:333730/149‑1 (MQ=255)
                                                                                                                                   aaaaaaTTCACCTTTTGTTTTCACACTGAAGTGTCAGCGAGAGTAAAAAATTATGCTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACa     <  2:852588/148‑1 (MQ=255)
                                                                                                                                      aaaTGCACATTTGTTTTTCAAACTGAAGTGTCAGCGAGAGTAAAAAAATATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGAACTACTCCATAGAAAATCATATCAGTAGAGACATAc  >  1:760434/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAC  >  NC_000913/1652227‑1652509

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