Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,179,964 T→A 5.5% intergenic (+381/‑32) rplA → / → rplJ 50S ribosomal subunit protein L1/50S ribosomal subunit protein L10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,179,9640TA 5.5% 646.1 / 16.3 197intergenic (+381/‑32)rplA/rplJ50S ribosomal subunit protein L1/50S ribosomal subunit protein L10
Reads supporting (aligned to +/- strand):  ref base T (83/103);  new base A (3/8);  total (86/111)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 200 total aligned reads displayed.

CGAAGACCGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATGACT  >  NC_000913/4179816‑4180112
                                                                                                                                                    |                                                                                                                                                    
   aGACCGCAGGAGTTTCGCAAGAAACTTATTCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTtc                                                                                                                                                  <  2:476108/150‑1 (MQ=255)
      ccGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTctct                                                                                                                                               <  2:584335/150‑1 (MQ=255)
      ccGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTctct                                                                                                                                               <  1:272731/150‑1 (MQ=255)
      ccGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTctct                                                                                                                                               >  1:515665/1‑150 (MQ=255)
      ccGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTctct                                                                                                                                               >  1:622363/1‑150 (MQ=255)
       cGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTTTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTg                                                                                                                                              <  2:600272/150‑1 (MQ=255)
       cGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTg                                                                                                                                              <  1:1167688/150‑1 (MQ=255)
         cAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGc                                                                                                                                            <  1:968855/150‑1 (MQ=255)
         cAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAAAATCTCTgga                                                                                                                                            <  1:692625/150‑2 (MQ=255)
           ggAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCaa                                                                                                                                          >  1:841775/1‑150 (MQ=255)
           ggAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAAAATCTCTggaac                                                                                                                                          <  1:614893/150‑4 (MQ=255)
           ggAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAAAATCTCTggaac                                                                                                                                          <  1:339911/150‑4 (MQ=255)
             aGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAAc                                                                                                                                        <  1:813516/150‑1 (MQ=255)
              gTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACa                                                                                                                                       <  2:635553/150‑1 (MQ=255)
              gTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACa                                                                                                                                       <  1:145724/150‑1 (MQ=255)
                   gCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAg                                                                                                                                  >  2:825338/1‑150 (MQ=255)
                    cAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAgg                                                                                                                                 <  2:550241/150‑1 (MQ=255)
                     aaGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGa                                                                                                                                >  1:833573/1‑150 (MQ=255)
                      aGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAg                                                                                                                               <  1:711495/150‑1 (MQ=255)
                         aaCTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCaa                                                                                                                            >  2:994431/1‑150 (MQ=255)
                          aCTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCaaa                                                                                                                           <  1:162296/150‑1 (MQ=255)
                            ttAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGc                                                                                                                         <  1:1257458/150‑1 (MQ=255)
                            ttAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGc                                                                                                                         >  2:494678/1‑150 (MQ=255)
                            ttAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGc                                                                                                                         <  1:57149/150‑1 (MQ=255)
                                 ccccTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATg                                                                                                                    <  2:892713/150‑1 (MQ=255)
                                   ccTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGc                                                                                                                  <  1:221363/150‑1 (MQ=255)
                                     tGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAttt                                                                                                                                                    >  2:871709/1‑114 (MQ=255)
                                      gCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCttt                                                                                                               <  2:682968/150‑1 (MQ=255)
                                        gTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTaa                                                                                                             >  2:1353766/1‑150 (MQ=255)
                                         tAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTaaa                                                                                                            <  1:63115/150‑1 (MQ=255)
                                           gACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCttt                                                                                                               >  1:30707/1‑145 (MQ=255)
                                           gACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATc                                                                                                          >  2:1029950/1‑150 (MQ=255)
                                           gACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATc                                                                                                          >  2:1168717/1‑150 (MQ=255)
                                             cGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCtt                                                                                                        >  2:1086946/1‑150 (MQ=255)
                                              ggTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTc                                                                                                       <  2:1291063/150‑1 (MQ=255)
                                               gTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCa                                                                                                      <  2:1238061/150‑1 (MQ=255)
                                               gTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCa                                                                                                      <  2:321664/150‑1 (MQ=255)
                                               gTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCa                                                                                                      <  1:25468/150‑1 (MQ=255)
                                                tGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCaa                                                                                                     <  1:458674/150‑1 (MQ=255)
                                                 gACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAg                                                                                                    <  1:757316/150‑1 (MQ=255)
                                                 gACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAg                                                                                                    <  2:255522/150‑1 (MQ=255)
                                                 gACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAg                                                                                                    <  1:865699/150‑1 (MQ=255)
                                                 gACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAg                                                                                                    <  1:884784/150‑1 (MQ=255)
                                                   cAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGac                                                                                                  >  2:1098546/1‑150 (MQ=255)
                                                     gAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGc                                                                                                                                            >  1:1398834/1‑106 (MQ=255)
                                                     gAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGc                                                                                                                                            <  2:1398834/106‑1 (MQ=255)
                                                        cGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGacaaaca                                                                                             <  2:590056/150‑1 (MQ=255)
                                                          cTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAg                                                                                           >  2:1390158/1‑150 (MQ=255)
                                                               attattCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGc                                                                                          <  2:163016/146‑1 (MQ=255)
                                                                ttattCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGAttg                                                                                     >  2:1342024/1‑150 (MQ=255)
                                                                    tCTTTTATATTCTGGCTTGTTTCTGCTCACCTTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGc                                                                                 <  2:1244724/150‑1 (MQ=255)
                                                                     cTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCt                                                                                >  2:62439/1‑150 (MQ=255)
                                                                      ttttATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTg                                                                               <  2:694930/150‑1 (MQ=255)
                                                                      ttttATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTg                                                                               <  1:1426093/150‑1 (MQ=255)
                                                                         tataTTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAg                                                                            <  2:715556/150‑1 (MQ=255)
                                                                            atTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCa                                                                         >  2:205721/1‑150 (MQ=255)
                                                                               cTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAg                                                                                                                          >  2:228173/1‑98 (MQ=255)
                                                                               cTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCg                                                                      >  2:703687/1‑150 (MQ=255)
                                                                                tGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTTCTGAAGTCAGCGa                                                                     <  1:896254/150‑1 (MQ=255)
                                                                                  gCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAg                                                                   >  2:480330/1‑150 (MQ=255)
                                                                                  gCTTGTTTCTGCTCACCATAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAg                                                                   >  1:394758/1‑150 (MQ=255)
                                                                                   cTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGt                                                                  >  1:913511/1‑150 (MQ=255)
                                                                                      gTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGAt                                                                                       >  2:384152/1‑126 (MQ=255)
                                                                                      gTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGc                                                               >  1:70184/1‑150 (MQ=255)
                                                                                       tttCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGacaaaca                                                                                             >  1:1359504/1‑119 (MQ=255)
                                                                                        ttCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCa                                                             <  2:76739/150‑1 (MQ=255)
                                                                                            gCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAgg                                                         >  1:167388/1‑150 (MQ=255)
                                                                                               cACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCACACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgc                                                      >  2:993431/1‑150 (MQ=255)
                                                                                               cACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgc                                                      >  2:1428616/1‑150 (MQ=255)
                                                                                               cACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgc                                                      >  2:808910/1‑150 (MQ=255)
                                                                                                 ccGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcgc                                                    <  2:834014/150‑1 (MQ=255)
                                                                                                 ccGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcgc                                                    >  1:1279841/1‑150 (MQ=255)
                                                                                                 ccGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcgc                                                    >  2:35256/1‑150 (MQ=255)
                                                                                                 ccGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcgc                                                    >  2:519024/1‑150 (MQ=255)
                                                                                                 ccGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcgc                                                    >  1:108696/1‑150 (MQ=255)
                                                                                                   gTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTg                                                  <  1:546571/150‑1 (MQ=255)
                                                                                                        tAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGc                                             >  2:1387636/1‑150 (MQ=255)
                                                                                                         aaGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCa                                            >  2:830530/1‑150 (MQ=255)
                                                                                                           gACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCagt                                          <  1:6430/150‑1 (MQ=255)
                                                                                                                tctctcCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTg                                     <  1:1387636/150‑1 (MQ=255)
                                                                                                                  tctcCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCg                                   <  1:1025797/150‑1 (MQ=255)
                                                                                                                   ctcCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCgg                                  >  2:417989/1‑150 (MQ=255)
                                                                                                                     ccGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGAt                                >  2:101818/1‑150 (MQ=255)
                                                                                                                       gTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTc                              >  1:111291/1‑150 (MQ=255)
                                                                                                                       gTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTc                              >  1:116473/1‑150 (MQ=255)
                                                                                                                          tGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCg                           <  1:205721/150‑1 (MQ=255)
                                                                                                                              ggagTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGc                       <  1:994131/150‑1 (MQ=255)
                                                                                                                               gagTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCg                      >  1:1227497/1‑150 (MQ=255)
                                                                                                                                         tttgccAGAGAAAATCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGa            >  1:180524/5‑150 (MQ=255)
                                                                                                                                          gTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGAt           >  2:1149348/1‑150 (MQ=255)
                                                                                                                                          gTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGAt           >  1:647000/1‑150 (MQ=255)
                                                                                                                                             ccAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATaaa        <  2:856011/150‑1 (MQ=255)
                                                                                                                                              cAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATaaaa       >  1:486425/1‑150 (MQ=255)
                                                                                                                                               agagATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAAt      <  1:101818/150‑1 (MQ=255)
                                                                                                                                               agagATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAAt      >  2:905294/1‑150 (MQ=255)
                                                                                                                                               agagATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAAt      >  1:1392258/1‑150 (MQ=255)
                                                                                                                                                gagaTTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATg     >  1:230847/1‑150 (MQ=255)
                                                                                                                                                gagaTTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATg     >  2:183486/1‑150 (MQ=255)
                                                                                                                                                gagaTTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATg     <  2:40712/150‑1 (MQ=255)
                                                                                                                                                  gaTTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATGAc   <  1:1149348/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGAAGACCGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATGACT  >  NC_000913/4179816‑4180112

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