Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,652,362 A→G 5.5% pseudogene (812/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,3620AG 5.5% 218.5 / 2.8 74pseudogene (812/1158 nt)intQpseudogene, Qin prophage; phage integrase family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (2/2);  total (35/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 = 4.32e-01

CAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAC  >  NC_000913/1652215‑1652509
                                                                                                                                                   |                                                                                                                                                   
cAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCTAAGAGCATATCTTTGATGTTCACTTGAGAGAGTATGGCAGAACAGAAAACCaaaaa                                                                                                                                                    >  1:699677/1‑149 (MQ=255)
cAAAAAAACAGGCTGTAACTAATCTGGTCATTTATATTATTAAGCCAGCAATCGACTCTCTCGGTAGTCAGATGACATTATCGTTACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACaaaaa                                                                                                                                                    <  2:351675/149‑1 (MQ=255)
 aaaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                                                                                                                   >  2:426650/1‑149 (MQ=255)
  aaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAATTg                                                                                                                                                  <  2:37690/149‑1 (MQ=255)
  aaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATg                                                                                                                                                  >  1:1150881/1‑149 (MQ=255)
            cGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                                                                                                                                        <  2:422186/149‑1 (MQ=255)
             gggAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                                                                                                                                       <  2:954904/149‑1 (MQ=255)
               gAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttttt                                                                                                                                     <  2:516602/149‑1 (MQ=255)
                 aCCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAGGTGCATTTTTGTTTTTct                                                                                                                                   <  2:300003/149‑2 (MQ=255)
                 aCCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAGGTGCATTTTTGTTTTTct                                                                                                                                   <  2:207832/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                                                                                                                               >  1:621565/1‑149 (MQ=255)
                     aTCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                                                                                                               <  2:894208/149‑1 (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                                                   <  2:549614/134‑1 (MQ=255)
                                                                                                                tCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGa                                                   >  1:549614/1‑134 (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)
                                                                                                                                                   |                                                                                                                                                   
CAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAC  >  NC_000913/1652215‑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: