Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,328,969 T→G 5.7% T277P (ACG→CCG)  yfaQ ← tandem DUF2300 domain‑containing protein YfaQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,328,9690TG 5.7% 212.3 / 4.9 70T277P (ACG→CCG) yfaQtandem DUF2300 domain‑containing protein YfaQ
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base G (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCAC  >  NC_000913/2328823‑2329118
                                                                                                                                                  |                                                                                                                                                     
ctctTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCcgtcgt                                                                                                                                                    >  1:768268/1‑150 (MQ=255)
 tctTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTa                                                                                                                                                   <  1:779889/150‑1 (MQ=255)
 tctTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTa                                                                                                                                                   <  2:1005188/150‑1 (MQ=255)
        gATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGACGGCAGTAAATCCCGGCGTActtaatg                                                                                                                                            <  1:156544/150‑6 (MQ=255)
         aTCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTACACGTCACACGATAGCGCCCATTAACTACGCCGGGATTTACTGCCCTCGTACTCTTTTc                                                                                                                                           >  2:845314/1‑150 (MQ=255)
            gagaTGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGc                                                                                                                                        >  2:1388336/1‑150 (MQ=255)
              gaTGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCtg                                                                                                                                      >  2:277786/1‑150 (MQ=255)
                tGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCtgtg                                                                                                                                    <  1:901166/150‑1 (MQ=255)
                    aCTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGAttt                                                                                                                                <  1:642876/150‑1 (MQ=255)
                    aCTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGAttt                                                                                                                                <  2:1401470/150‑1 (MQ=255)
                     cTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGGCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGAtttt                                                                                                                               <  1:690308/150‑1 (MQ=255)
                      tGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTttc                                                                                                                              <  1:529379/150‑1 (MQ=255)
                       gCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTttct                                                                                                                             <  2:1207962/150‑1 (MQ=255)
                       gCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTttct                                                                                                                             <  1:979642/150‑1 (MQ=255)
                            ttCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCTGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGc                                                                                                                        <  1:106002/150‑1 (MQ=255)
                             tCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAgag                                                                                                                       <  1:275677/150‑3 (MQ=255)
                                   tGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCcgtcg                                                                                                                                                     >  2:1095992/1‑114 (MQ=255)
                                   tGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCcgtcg                                                                                                                                                     <  1:1095992/114‑1 (MQ=255)
                                    gCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCATCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCg                                                                                                                >  2:790370/1‑150 (MQ=255)
                                      ccTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCggg                                                                                                              >  1:990446/1‑150 (MQ=255)
                                          ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCaaa                                                                                                          <  1:1331923/150‑1 (MQ=255)
                                                cAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTc                                                                                                    >  1:1407820/1‑150 (MQ=255)
                                                cAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTc                                                                                                    >  2:1264057/1‑150 (MQ=255)
                                                cAGTAAGGTCGTCTCAACGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTc                                                                                                    >  2:1477431/1‑150 (MQ=255)
                                                           tctcACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTacac                                                                                         <  1:902200/150‑1 (MQ=255)
                                                            ctcACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTAcaca                                                                                        >  1:704670/1‑150 (MQ=255)
                                                             tcACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACAt                                                                                       >  2:905708/1‑150 (MQ=255)
                                                             gcACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACAt                                                                                       <  1:633792/149‑1 (MQ=255)
                                                                 cGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTTATTGTATAGATTGCAATAAGATCGAGATCTACACATTGcc                                                                                   <  2:10630/150‑1 (MQ=255)
                                                                      cTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTg                                                                              <  1:1405591/150‑1 (MQ=255)
                                                                       tGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTgg                                                                             <  2:751929/150‑1 (MQ=255)
                                                                       tGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTgg                                                                             >  1:390157/1‑150 (MQ=255)
                                                                        gACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGc                                                                            >  1:1342628/1‑150 (MQ=255)
                                                                          cAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCaa                                                                          >  2:109763/1‑150 (MQ=255)
                                                                               gTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCCTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGc                                                                     <  1:247942/150‑1 (MQ=255)
                                                                                   tgtgGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCg                                                                 <  1:443592/150‑1 (MQ=255)
                                                                                    gtgGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCgg                                                                <  2:1628261/150‑1 (MQ=255)
                                                                                       gTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCAc                                                             >  1:1389082/1‑150 (MQ=255)
                                                                                         tACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTg                                                           >  1:1550327/1‑150 (MQ=255)
                                                                                                aGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAAc                                                    >  1:1425604/1‑150 (MQ=255)
                                                                                                aGTAAACGTCACAAGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAAc                                                    >  2:600972/1‑150 (MQ=255)
                                                                                                 gTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCATCCTGCGTCGGCACTGGCAAAAcc                                                   <  2:1257116/150‑1 (MQ=255)
                                                                                                   aaaCGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGc                                                               <  1:360618/136‑1 (MQ=255)
                                                                                                   aaaCGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGc                                                               >  2:360618/1‑136 (MQ=255)
                                                                                                   aaaCGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCg                                                 >  1:961747/1‑150 (MQ=255)
                                                                                                     aCGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGAt                                               >  1:592466/1‑150 (MQ=255)
                                                                                                     aCGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGAt                                               <  1:805354/150‑1 (MQ=255)
                                                                                                       gTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTa                                             >  2:630811/1‑150 (MQ=255)
                                                                                                          acacGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGTGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGc                                          >  1:939518/1‑150 (MQ=255)
                                                                                                          acacGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGTGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGc                                          >  1:933723/1‑150 (MQ=255)
                                                                                                          acacGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGc                                          >  2:387869/1‑150 (MQ=255)
                                                                                                          acacGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGc                                          <  1:244357/150‑1 (MQ=255)
                                                                                                           cacGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCa                                         >  2:564684/1‑150 (MQ=255)
                                                                                                             cGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCAcc                                       >  1:414318/1‑150 (MQ=255)
                                                                                                              gATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCg                                      <  2:1161306/150‑1 (MQ=255)
                                                                                                                   cgcCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCa                                 <  1:478139/150‑1 (MQ=255)
                                                                                                                     ggaaaagaGTACGACGGCAGTAAATCCCGGCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCATCGCATa                               >  2:711332/8‑150 (MQ=38)
                                                                                                                        aTTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACg                            <  2:768268/150‑1 (MQ=255)
                                                                                                                         ttAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACgg                           <  2:390157/150‑1 (MQ=255)
                                                                                                                            aGTACGACGGCAGTAAATCCCGGCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGtt                        >  1:931685/1‑150 (MQ=255)
                                                                                                                             gTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTa                       >  2:1018394/1‑150 (MQ=255)
                                                                                                                                cGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAAccc                                                  <  2:1605070/120‑1 (MQ=255)
                                                                                                                                cGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAAccc                                                  >  1:1605070/1‑120 (MQ=255)
                                                                                                                                cGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGc                                  >  2:441945/1‑136 (MQ=255)
                                                                                                                                cGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGc                                  <  1:441945/136‑1 (MQ=255)
                                                                                                                                cGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCAAAACggtactaac                    >  1:980738/1‑143 (MQ=255)
                                                                                                                                  ccGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTAcc                  >  2:753236/1‑150 (MQ=255)
                                                                                                                                      gATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGatt              >  1:1556354/1‑150 (MQ=255)
                                                                                                                                         ttACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGc           >  1:807752/1‑150 (MQ=255)
                                                                                                                                          tACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCt          <  1:1264057/150‑1 (MQ=255)
                                                                                                                                          tACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCt          <  1:1018394/150‑1 (MQ=255)
                                                                                                                                             tGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGt       >  1:702314/1‑150 (MQ=255)
                                                                                                                                              cccGGCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCATCGCATAACGGGTTAGTACCGTTTTGCTGGTg      <  1:711332/149‑1 (MQ=255)
                                                                                                                                                cgtcgtACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGcc    <  2:592466/150‑1 (MQ=255)
                                                                                                                                                cgtcgtACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGcc    <  1:1016568/150‑1 (MQ=255)
                                                                                                                                                 gtcgtACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCa   >  2:1313266/1‑150 (MQ=255)
                                                                                                                                                  tcgtACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCAc  <  1:1388336/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCAC  >  NC_000913/2328823‑2329118

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